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CXXXITI. Microscopic Studies on Beryl Crystals.—I1. Dislocations and 
the Growth of {1010} Prism Faces. ) 


By L. J. GRIFFIN, 
H. H. Wills Physical Laboratory, Bristol*. 


[Received July 20, 1951.] 
[Plates XLV.-LII.] 


SUMMARY. 


It is shown that the type of growth predicted by the dislocation theory 
of growth takes place on {1010} faces of beryl. The growth layers, in 
general, produce a complex topography. The behaviour of layers at 
surface obstacles is discussed. Examples are given of the propagation 
of a dominant layer group past lines of dislocations, and inactive 
dislocation groups, and the mode of propagation by cross-linking is 
discussed. Typical examples of lineage boundaries and limited slip 
zones are shown; these features are observed as lines of dislocations. The 
dependence of growth rate of the layers on crystallographic orientation 
is considered. The types of layer structure observed on natural and 
synthetic crystals are compared. It is shown that the dislocation theory 
of growth can be applied to give such information as the supersaturation 
and temperature at which growth took place. 


§1. GENERAL. 


THE areas between the uni-molecular steps observed on prism faces of 
beryl (Griffin 1950, 1951) are in general plane. In some cases, however, 
the uni-molecular steps stop abruptly, and it is possible to describe a 
path from the low to the high side of a step without at any point crossing 
a discontinuity in height. The point of termination of the step must 
therefore mark a point of singularity of the surface corresponding to the 
point of emergence on the face of a dislocation with a screw component 
normal to the face. Any dislocation of this kind must necessarily mark 
one end of a step, which can only run off to the edge of the face or end 
at another screw dislocation. 


§2. SPIRALS AND Loops. 


A case is shown in fig. 1 (Pl. XLV.) of a step, attached to a single 
dislocation, which takes a spiral form. Such behaviour has been predicted 
on purely theoretical grounds for the case of the advance of a layer 
attached to a single dislocation when the rate of advance of the layer is 
dependent on the crystallographic orientation (Burton, Cabrera and 
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Frank 1951). The agreement between theoretical prediction and 
observation is so close that there can be no doubt but that the layer 
in fig. 1 (Pl. XLV.) has been winding itself up during growth. 

The case of an edge joining two dislocations of opposite hand is shown 
in fig. 2 (Pl. XLVI), this edge lying in the centre of the figure. The 
behaviour predicted by Frank was that the inner unclosed loop would 
expand during growth so as to generate successive closed loops, typified 
by the outer loop of fig. 2 (Pl. XLVI.). (The triangular area cutting 
into the closed loop of fig. 2 (Pl. XLVI.), probably consists of an 
overgrowth of foreign crystalline material.) 

The dependence of growth rate upon crystallographic orientation of 
a layer can be represented by means of a two-dimensional polar diagram 
(Frank 1950). A close packed direction (Burton and Cabrera 1949) 
corresponds to a minimum growth rate and will manifest itself as an 


Fig. 3. 


almost straight edge. The dislocation pair in fig. 2 (PI. XLVI.) has been — 
operating unhampered by the presence of other dislocations so that this 
case is almost ideal for consideration of growth rates. In this case the 
longitudinal edges must correspond to a close-packed direction. (These 
almost straight longitudinal edges are parallel to the c-axis and are most 
useful in defining the orientation of any area.) : 

The two transverse edges of the outer closed loop are very nearly 
ares of the same circle, so that the rate of growth must be independent 
of orientation for the range covered by these edges (--30°). It is therefore 
possible to construct a polar diagram as shown in fig. 3. This diagram 
coincides with the “ razor-blade ” outline over the circular ares at the 
ends and at the mid points of the longitudinal edges. The dotted sections 
of the diagram are uncertain, and it is only possible to say, from the 
present information, that the polar diagram must lie outside these 
sections. 
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When two dislocations of the same hand are closer together than a 
certain distance (equal to 47 times the size of the critical nucleus) the 
theory predicts that the dislocations will work together cooperatively 
and produce steps which form two arms of a. spiral. Examples of this 
behaviour have also been observed. 

The examples so far quoted are the ideally simple cases and are 
comparatively seldom observed. More usually a complex layer system 
is seen on any area of a face, and we shall now consider some characteristic 
examples. 

§3. BEHAVIOUR AT OBSTACLES. 


A typicai example of the behaviour of a layer system at a surface 
obstacle is shown in fig. 4 (Pl. XLV.). This layer system is of 
considerable extent and a small part of it has met with an obstacle to 
its steady growth, this obstacle being probably a crystal inclusion. The 
layers are forced to divide into two parts by the obstacle and as they 
advance towards the right of the figure the obstacle decreases in extent. 
It is therefore necessary for the layers to “ fill in”? by inward growth of 
the straight edges. Due to the fact that growth is more rapid parallel 
to the c-axis than perpendicular to it, the two sections of the system only 
unite at some distance beyond the obstacle. The points to be noted are :— 

(a) The rapid elimination of the re-entrant angle at the junction of the 
two sections. As the re-entrant is parallel to the direction of rapid 
growth this would be expected on almost any theory of growth. 

(6) The relative rate of growth parallel and perpendicular to the 
c-axis. The ratios of these rates are dependent upon the local conditions, 
being 1-7 and 2-1 for the two parts of this system, and in other similar 
cases 2-9 and 2-76, 2-36 and 2:37, 2-1 and 1-9. 

(c) Each edge is a contour, with a vertical interval of 7-9A. The 
obstacle thus results in a triangular V-shaped valley with a maximum 
depth of 38 molecules (300 AN 

An obstacle of this kind appears to result. in an increased curvature. 
of the edges close to it, but as the effect on the two sections of the system 
is asymmetrical it is difficult to determine whether the effect is regular. 

This difficulty can be overcome in the case of fig. 5 (Pl. XLVI.) in 
which the layers advancing from the right are split into three sections by 
a succession of crystal inclusions, which again act as obstacles. The 
central system formed in this way is affected equally on both sides and 
it can be immediately seen that it behaves as if a “ virtual source ” was 
present. 

The approximate position of the virtual source can be found by 
producing back lines drawn through the points of junction of the curved 
and straight edges ; it lies on the upper edge of the narrow neck between 
the two obstacles. The ratio of the radius of curvature of the transverse: 
edges of this small system to their distance from this ‘‘ source ”’ is constant 
and equal to unity. Moreover, this relation is still maintained, after 
junction with the lower section, for those parts of the first four layers of 
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the joint system which originated from the central section. eile 
expanding from the * virtual source : thus behaves in precisely the same 
way as the closed loop in fig. 2 (Pl. XLVL.). . 

The behaviour of layer systems at obstacles provides valuable 
confirmatory evidence that the observed steps are normally uni-molecular. 
Thus in fig. 5 (Pl. XLVI.) the central and lower sections first combine 
again to form a single system. Counting from the last edge before 
division commences on the right, there are 37 edges in both the lower 
and central sections before a junction takes place. Similarly, considering 
the junction beyond the third obstacle there are 57 edges in both the upper 
and lower sections between the first division of the edge and the first 
joint edge. Similarly in fig. 4 (Pl. XLV.) in both the upper and lower 
systems there are 70 edges between the first division of the incident layer 
and the first joint layer. Such behaviour would be very difficult to 
comprehend for steps of more than one molecule high, for one would 
expect the steps in such circumstances either to split into smaller units 
or to combine to form multiple steps. 


Fig. 6. 


$4. Propagation Past Inactive Distocatrons. 


The theory predicted that certain dislocations and dislocation groups 
on a face could be more active than others, resulting in those active 
dislocations “dominating” less active dislocations (i.e. the layers 
emanating from active dislocations propagate past ‘* dominated ” groups 
or single dislocations with only a slight delay and a slight change in 
outline). The mode of propagation in such cases and the type of cross- 
linking which occurs is shown in fig. 6. The layers are advancing from 
right to left and in (i) one layer has just passed the dislocation. In (ii) the 
layer originating from the inactive dislocation has grown in the normal 
way, while the layer advancing from the right is quite close to it. The 
kink in the layer on the left is being eliminated as the layer extends. 
In (iii) the incident layer is being held up at the dislocation, as this layer 
is now close to the region of deformed lattice around the dislocation, 
and an effective small repulsion exists between the dislocation and the 
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incident layer. By stage (iv) sufficient molecules have built on to the- 
requisite section of the incident layer, probably because of a statistical 
fluctuation, to enable the layer to link across with the layer running 
from the dislocation. The dislocation is consequently attached to the 
shorter edge, which originally formed part of the incident layer. 
Meanwhile the kink in the left-hand layer has been almost eliminated. 
This type of behaviour can be repeated indefinitely. 

Fig. 7 (Pl. XLVIL.) and fig. 8 (Pl. XLVII.) show two examples of this. 
kind of behaviour from the ve number observed. In fig. 7 (Pl. XLVII.) 
a dominant layer system is advancing from the left and meets with 
. a line of dislocations which are all of the same hand and all the stages 
of cross linking drawn in fig. 6 can be seen. 

Fig. 8 (Pl. XLVII.) illustrates more clearly the nine of the layers. 
after the dislocations have been left behind, the dominant layer system 
moving from right to left in this case. In both figures it can be seen that 
the layers pass on with only very slightly modified outline, and, that 
some five layers after the dislocations are encountered, the layer outlines 
show no trace of the action of the dislocation groups. 

Fig. 8 (Pl. XLVI.) is also of interest because of the evidence it provides 
of varying activities of layer systems, a subject which will be dealt with 
fully in a later paper. On the left of this figure there is a sudden con- 
siderable increase in layer spacing.. Now layer spacing can most simply 
be equated against supersaturation, but in a case like this it is not likely 
that the supersaturation varied over this area. It would rather appear 
that the widely spaced layers on the left were produced by a dislocation 
group of low activity, this low degree of activity being due to a low 
supersaturation or other causes. At a later date in the history of the 
growth of this face another system has become active (possibly because 
of an increase in supersaturation, or the creation of new dislocations), 
producing more closely spaced layers which have gradually obliterated 
the old layer system, 7. e. the new system has become the dominant one 
in this area. Thus the area in fig. 8 (Pl. XLVII.) shows the farthest 
of advance of the new system before growth ceased. 


§5. LINEAGE STRUCTURES. 


It is generally accepted that most crystals consist of a lineage structure, 
i.e. of -erystals blocks, sometimes macroscopic, which are slightly 
disoriented with respect to one another. The atomic structure of the 
boundaries between such blocks can be explained most satisfactorily in 
terms of dislocation systems lying in the boundaries, for dislocations 
necessarily produce a rotation of the lattice at distances far from the 
dislocation. Such a picture has proved valuable in investigations on 
grain boundaries in metals, but the evidence in favour of the presence 
of dislocations which results from such work, is necessarily indirect. 
On faces of beryl crystals, however, the structure of such boundaries 
can be observed directly. An example of a well marked lineage boundary 
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is shown in fig. 9 (Pl. XLVIII.). There are some 85 dislocations in the 
boundary, the left-hand part of which is partially obscured by foreign 
material, all the dislocations being of the same hand. The distance 
between dislocations, measured in atomic spacings, is equal to the angle 
of rotation between the two crystal blocks. In this case the angle of 
rotation across the boundary is approximately 2-8 minutes of arc. 

It is certain that these features consist of a line of dislocations. The 
layers would not curve round in this way if this feature was merely a 
surface step, for a considerable body of evidence shows that in all such 
cases the layer configuration would be quite different. For example, the 
behaviour at obstacles shown in fig. 4 (Pl. XLV.) and fig. 5 (Pl. XLVI.) 
can be contrasted with fig. 9 (Pl. XLVI.). Moreover, the layers running 
from lineage boundaries often terminate on neighbouring dislocations. 

The step marking this boundary is very obvious under narrow pencil 
illumination and this is understandable for, because of the dislocation 
structure, the higher end of the boundary must project 670 A. above the 
central region of fig. 9 (Pl. XLVIII.). - 

There are many tilt boundaries which are much smaller than the 
above example. One such case is shown in fig. 10 (Pl. XLY.). From 
each dislocation of the lower line a spiral edge starts, the straight sections 
being so closely packed that they are barely resolved. An edge also starts 
from each dislocation of the upper line, which again is marked as a series 
-of dots, the two systems combining to form a large spiral. (The area 
blocking out part of the spiral is probably composed of foreign material.) 
The second, lower, set of straight edges forming the spiral is not parallel 
to the original set, a very unusual feature. These edges end in an 
irregular boundary which has been formed by the main growth region. 
From the position and behaviour of the series of edges above the dislocation 
‘system it appears probable that these may originally have formed part 
of the spiral. The angle of rotation at the two parallel boundaries is 
approximately 8 seconds of are in each case. 

A very heavily lineaged area is shown in fig. 11 (Pl. XLVIII.). Almost 
every layer in this figure runs from a lineage boundary, these boundaries 
lying close together in the lower part of the figure. (The obvious black 
feature which runs diagonally across the figure is a crack.) The layers 
are so closely spaced that the longitudinal edges are not resolved. On 
one side of this lmeaged area there are some 140 dislocations of one hand, 
with some 180 dislocations of opposite hand on the other side. The 
result of these two groups being almost equal is that the difference in 
orientation of the lattice on either side of the area will be small, and will 
in fact amount only to approximately 20 seconds of arc. The actual 
difference in orientation produced by either group is much greater, 
being about 2-4 minutes of arc. 

The position and nature of these lineages make it possible that they 
originated in the form of limited slip-zones, a type of feature discussed 
immediately below. The constituent halves of these zones would grow 


Microscopic Studies on Beryl Crystals 1343 


apart if the dislocation lines were set at a slight angle to each other, so 
that each half of a zone could eventually act independently as a lineage 
boundary. 

It will be noticed that the examples of lineage boundaries which have 
been shown are all parallel to the c-axis of the crystal. This appears to 
be true for the majority of such features which have been observed. 


§6. Lonrep SLiP-zonEs. 


Another type of dislocation distribution very often observed is that 
corresponding to a limited slip-zone. Such features will be composed 
of a line of dislocations which is divided about the centre into two equal 
groups of opposite hand, the line lying along the trace of the 
appropriate glide plane in the face. A feature of this kind which lies 
parallel to the c-axis is shown in fig. 12 (Pl. XL{X.). This is unusually 
large and consists of two groups of approximately 90 dislocations each. 
(Evidence supporting the interpretation of these features as limited 
slip-zones will be given in a later paper. We are only concerned here 
with these slip-zones as typical dislocation groups.) 

An example of a limited slip-zone perpendicular to the c-axis is shown 
in fig. 13 (Pl. XLIX.). Isolated slip-zones of this orientation are much 
less common than cases parallel to the c-axis. However, a not uncommon 
feature consists of a number of contiguous small slip-zones which form a 
long, closely packed line of dislocations perpendicular to the c-axis. 
In fig. 13 (PL. XLIX.) the dislocation group has been sufficiently active 
to cause a local variation in the shape of the dominant layer system 
which is advancing from the left. 


§7. Propacation or Layers past 4 DisLocaTION GROUP. 


The case of the propagation of layers past dominated dislocations of 
the same hand has already been discussed. The more general case consists 
of the propagation of a dominant layer system through an inactive group 
of dislocations, and there are some significant differences in behaviour in 
this case. 

In fig. 14 (PL. L.) the dominant system is advancing from the right 
and meets a group of dislocations which are only sufficiently active 
to cause local variations in the outline of the layers. Almost all the steps 
in the bottom part of this figure run between pairs of dislocations and 
several of the layers tend toward the characteristic razor-blade shape. 

Several steps in this group have a reversed S-shape. Such behaviour 
is not contrary to the theory but, as we shall see, is merely due to the 
dominant characteristics of the incident system. For example, a very 
acute elbow runs from a single dislocation in the upper part of the figure. 
It is obvious that this edge, which originally curved downwards, must form 
a re-entrant of this type on junction with the incident layers. Moreover, 
this re-entrant cannot be rapidly eliminated as it lies in the direction 


corresponding to a minimum growth rate. 
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In order to illustrate some of the points which arise a part of the 
dislocation group is isolated and drawn in full line in fig. 15 (i). It is 
necessary to consider a limited area only in order to simplify the actual 


Fig. 15 (i). 


Fig. 15 (iii). 


Sy 


z= (eS 


rr 
J/ 


Microscopic Studies on Beryl Crystals 1345 


behaviour. Thus, the curved edges running down on the left of 
fig. 14 (Pl. L.) terminate at dislocations which are all of one 
hand, These terminations have not been shown in order to keep the 
diagrams down to a reasonable size, and also because only a simple type 
of cross linking takes place in this area. Hence, over the area shown 
in fig. 15, these curved edges can be considered to propagate regularly 
without any considerable deviation from the true conditions. 

In these diagrams the propagation of the layers has not been drawn 
exactly to scale but only approximately so. This means that deviations 
from the order of events shown in these diagrams might have occurred 
had the crystal continued to grow, but the aim is to show the general 
way in which propagation by cross-linking takes place. 

In (i), (ii) and (iii) five successive positions of the layers are shown. 
The advance of the layers from their original position in any of these 
diagrams is shown by a dotted line, and the position of layer I, for example, 
after an advance will be referred to as I’. 

Let us consider the position of the layers in fig. 15 (i) after a slight 
advance has taken place. A’ is held up at the parent dislocation of B’, and 
cross-linking will soon take place. Similarly, C’ is now very close to the 
dislocation of D’. E’ and F’ show no important change (G’ is growing 
normally, but G was initially of an anomolous shape due to' causes to 
be discussed in a later paper). The advance of H has caused cross-linking 
to take place with I, so that the joint layer I’H’ and a new double 
dislocation loop are formed. J’ and K’ show little change, but L and 
M have cross-linked to form the layer L’M’ and a double dislocation 
loop L’ has been formed. A further stage of advance of these layers is 
shown in full line in (ii). The lowest four layers have now cross-linked 
to give O and P and Q and R, R continuing with a considerable kink 
because of the cross-linking. E’ has cross-linked with F’ to give § and 
the layer T which is held up at the lower dislocation of the loop V. 
J’ has linked with K’ to form the double dislocation loop X, and cross- 
linking of this loop with W is imminent. No marked change has taken 
place in the other layers. 

Another advance immediately causes T to cross-link with the double 
loop V, so that a dislocation is now attached to the right-hand section 
of T, and the layer T’U’ carries on with an obvious re-entrant. W and 
X have linked so that the two dislocations of the loop X now function 
separately. Y and a have linked so that the left-hand part of Y’ is 
attached to a dislocation. Meanwhile the joint layer a’Y’ almost 
immediately links with 6’ so that a’ is again attached to a dislocation, 
while the joint layer b’Y’ passes on with a large re-entrant due to the 
rapid cross-linkings. The edge b’Y’ being nowhere attached to a dis- 
location (in this figure at least), is thus the first layer to pass through 
this group. 

The position after a final advance is shown in fig. 15 (iii). O’, P’, Q’ have 
advanced undisturbed to the positions c, d and e. (In fact, of course, 
another layer would by now have advanced from below to cross-link 
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with c.) R’ has linked with the double loop T’, leaving the two dislocations 
attached to the layers e and f. T’ has linked with U’ to give the layers 
Landi. Layer L is held up by one of the dislocations of the loop m and 
linking is about to take place. Otherwise there is no considerable change. 

This group of dislocations which has been studied is both small and 
comparatively widely spaced. In the case of many of the larger closely 
packed groups it is obvious that the general outline of the layers will be 
complicated and varying very rapidly, so that any one layer in its passage 
through the group must very soon lose any semblance of composition 
from its original layer and will undergo a great many changes in outline. 

It can be seen from fig. 15 (i) and (iii) that the original and final positions 
of the layers are similar. This would be more marked if the original 
position of the layers did not show certain slight anomalies explicable 
in terms other that growth. A certain periodicity is bound to occur if 
the conditions remain uniform, and if the incident edges advance in a 
regular way, a cycle of quite short period should result. Normally the 
conditions may fluctuate sufficiently to prevent the establishment of such 
a cycle, although a semi-periodic behaviour should result. 

Fig. 15 shows that when dislocations, whatever their hand, are spaced 
with one or more layers between them, only a considerable kinking of the 
layer edges is likely to result. However, when successive layers run from 
dislocations of opposite hand, a succession of double-dislocation loops and 
single dislocation layers must be formed very rapidly. If a considerable 
number of such dislocations of opposite hand lie together, then a 
complicated pattern of looping and linking must ensue as any one layer 
progresses through the system. 


§8. AppITIONAL EviDENCE OF THE DEPENDENCE OF GROWTH 
RaTE ON ORIENTATION. 


In fig. 16 (Pl. L.) a dominant system is advancing from the left 
when it meets a barrier formed by a closely spaced line of dislocations. 
The main system is forced to spread downwards from the upper limit of 
this barrier, until it meets the independent layer system running from 
the dislocations. At the junction between the two systems a slightly 
curving edge results which runs at a slight angle to the c-axis and 
continually encroaches on the curved edges of the upper system. This 
encroachment continues until the slightly curved edges spread across to 
make a junction with the straight longitudinal edges (see fig. 17 
(Pl. LI.)) and so completely eliminate the curved edges. At this point | 
the edges are slightly concave in curvature and make an angle of 
approximately 25° with the c-axis. ; 

This junction continues for some 5 mms. along the face (i. €. some 
1,000 layers) and the edges gradually become straighter and progressively 
less inclined to the c-axis. A locally dominant system then altered the 
conditions, and these edges were replaced by the normal type of curved 
edges. The general conclusion to be drawn from this behaviour appears 
to be that edges inclined at small angles to the c-axis can grow more 
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rapidly than the normal curved edges. This conclusion is supported by 
observations of other similar events which gave a total range of 
inchnations of these rapidly growing edges of 4°-25°.. It is therefore 
possible to make additions to the polar diagram shown in fig. 3 to produce 
the more complete diagram, fig. 18. The dotted lines still indicate the 
uncertain sections. It should be noted that these polar diagrams are 
only strictly applicable to other beryl crystals if they have grown in 
identical conditions. Considerable variations are observed in the 
proportions of these diagrams, but the shape is usually of the same kind. 
However, examples have been observed when even the shape is quite 
different. 


Fig. 18. 


§9. SyNTHETIC EMERALD. 


In addition to naturally occurring beryl crystals originating from 
‘various sources, two synthetic emeralds were examined. These crystals 
had been grown by Nacken by a hydrothermal process (Van Praagh 
1948), the details of which have not been revealed. The faces of both 
-crystals were very rough on a microscopic scale and a regular structure 
could only be detected on some of the faces. In all cases the surfaces 
-were too irregular to allow very fine structure to be observed. 

All the faces suitable for observation showed a number of very similar 
hills standing out above the surface. A typical feature of this kind as 
-seen under oblique illumination is shown in fig. 19 (PI. LI). All 
the steps appear to be multi-molecular and the hill builds up to a definite 
point, the inference being that one central group of dislocations has been 
responsible for the growth of the hill. 

This type of behaviour can be contrasted with the typical plateaus seen 
.on natural crystals, as shown for example, in figs. 20 and 21 (PILL): 
These plateaus stand some wavelengths above the general surface level, 
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but they still retain the typical shape characteristic of a simple uni- 
molecular layer as shown in fig. 2 (PI. XLVI.). The main plateaus of 
these figures are very flat, very few uni-molecular steps being visible. 
The single spiral shown in fig. 1 (PI. XLV.) lies on the plateau of 
fig. 20 (Pl. L.) and its position can be identified by the gross black 
feature: This appeared to be one of the very few dislocations active on 
this plateau. 

The two other main points of difference between the structures on natural 
and synthetic crystals are :— 

(i.) The differing shapes of the layers. The longitudinal edges on the 
synthetic crystal have a considerable curvature; the significance of this 
point is discussed below.. 


Fig. 22. 


(ii.) The presence of a rough structure on the transverse gradients of 
the natural crystals. This rough structure actually consists of close-packed 
“rods ’—features which will be fully discussed in a later paper. It 
need only be said that they invariably occur on steep gradients parallel 
to the c-axis on the natural crystals studied. No trace of such features. 
was observed on the synthetic crystals. 

It is possible to construct a polar diagram of growth rate for these 
synthetic crystals. The transverse edges are again almost ares of the 
same circle, so that the growth rate must be independent of orientation 
over a range of +20°. The polar diagram is shown in fig. 22 where 
the dotted lines again indicate uncertain areas. 


S10. APPLICATION OF DisLocATION THEORY OF GROWTH. 


The theory predicts that for directions not far removed from a close- 
packed direction an edge will have a radius of curvature (1+2B) times 
its distance from the point of origin, where B= exp (¢/kT) and ¢ is 
2 neighbour-neighbour interaction energy. The radius of curvature of 
the almost straight edges of the outer loop of fig. 2 (Pl. XLVI.) is some 
120 times the distance from the centre. Therefore $/kT is about 6-2. 
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The average values of these constants for the two synthetic emeralds 
examined were :— 


B pikT 
Specimen 1 3°83 3-42 
Specimen 2 3-71 3°39 


‘These values are very similar and it appears probable that. both specimens 
have grown in almost identical conditions. If one assumes that both 
the natural and synthetic crystals have grown from similar media then 
the differing values of ¢/kT can be ascribed to variations in T, and the 
natural crystal must have grown at an absolute temperature of a little 
more than half that of the synthetic crystals. The value of Nacken’s 
‘Operating temperature is not known, but it was probably close to the 
critical temperature of water. The natural crystal must therefore have 
grown at a temperature of approximately 80°C. This may appear 
‘surprisingly low, but it is of interest to note that the temperature of 
growth of a-beryl crystal found in a pegmatite vein was 175° C., as 
determined by a liquid inclusion method (Ingersoll 1947). Moreover, 
the temperature of 80° C. can probably be taken as the lower limit, for 
the values of B were determined from measurements on multi-molecular 
steps, and such steps usually give results which are inaccurate unless 
one concludes that the values of B are greater than would hold for 
uni-molecular steps. 

It is possiblé to estimate the size of the critical nucleus (p,) by 
measurement of the distance between two dislocations of opposite hand 
which are linked by an edge which has been expanding outwards. The 
distance between the dislocations must be equal to, or greater than, 
Amp, where p, is defined by p,=a¢/kTo, where a is the inter-atomic spacing 
and o is the supersaturation. Knowing the value of ¢/kT we can 
deduce from fig. 2 (Pl. XLVI.) that the supersaturation in this area was 
about 0:5 per cent. 

The accuracy with which the supersaturation in any area can be 
determined will always be dependent on the presence of two dislocations 
with a spacing of just 47p,. This means that the closest spacing should 
be taken as giving the most accurate value, but it is sometimes difficult 
to be sure whether a given edge has been expanding. As an example 
there are tabulated a number of observations on another face of the 
erystal on which fig. 2 (PI. XLVI.) was taken. In this case an asterisk 
denotes loops which may possibly not have been expanding. 

From this table it would seem that the loops corresponding to values 
of 0:29 and 0:35 microns had not been expanding. Neglecting these 
two values and the rather high final value, the average value of p, over 
this region of the face was probably quite close to 0-49 microns. Taking 
d/kT=6-2 this is equivalent to an average supersaturation of about 


_ I per cent. 
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Pp, (microns) 


Active group oe 
‘Oo 


Groups of fig. 13 0-46 


0-41* 


Centre of very large dominant group 0-55 
0-29* 


0-75 


§11. CONCLUSIONS. 


The observations show that the dislocation theory of growth is not 
only correct in principle, but also in detail. This is the more surprising 
when it is recalled that this theory was derived primarily for homopolar 
crystals, formed from spherical or cubic molecules, growing from the 
vapour phase. In this case beryl is an ionic crystal with complex 
molecules and growth has taken place from solution. 

The dislocation theory of growth provides, for the first time, a 
quantitative account of crystal growth, and it is possible to learn a good 
deal about the conditions under which any crystal has grown from simple 
measusements on the faces. This theory also marks a considerable advance 
in the unification of growth phenomena and modern theories of the solid 
state, two fields which had previously been largely divorced. 

Lineage and transition surfaces have been explained for some time in 
terms of dislocation theory as arrays and grids of dislocations. The 
observations provide direct visual evidence that this explanation is 
correct. 
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EXPLANATION OF THE PLATES. 


Puate XLV. 
Kiger, 


The dislocation lies at the origin of the spiral. The kink in the right-hand 
edge of the spiral has possibly been caused by impurity obstructing 
steady growth. . 

Fig. 10. 


.The black area at the top is a hole in the crystal surface. The two boundaries 
lie a little down from this. This is a case where resolution has to be 
sacrificed to visibility. 


Puate XLVI. 
Fig. 2. 
The inner step-line runs between two dislocations of opposite hand, the points 
of emergence of these dislocations being imaged as obvious black dots. 


Puate XLVILI. 
Fig. 7. 

The line of inactive dislocations lies close to the centre of the figure, and its 
presence is made noticeable by the kinking of the edges close to it. The 
dislocations are shown as black dots. In the top part of the figure the 
layer edges are running parallel to the c-axis, but are too closely spaced 
to be resolved. 

s Fig. 8. ‘ 
The dark band at the bottom marks the edge of the face. The line of 


dislocations starts near the centre of the edge and from there runs upward 
and slightly to the right, almost to the top. 


Puate XLVITI. 
Bis. 9: 

The lineage boundary runs from left to right just below the centre, its left-hand 
extremity being partially obscured. To the left of the mid-point of the 
boundary the step-lines become much more closely spaced and less clearly 
visible, while the corresponding longitudinal sections are too closely 
packed to be resolved. 

Rigs li: 


The position of the different boundaries can be most easily found by following 
the downward path of the curved edges. There are at least six 
boundaries, all situated in the lower central part of this area. 


Puate XLIX. 
Fig. 12. 


The slip-zone is the line running horizontally from which 
all the curved edges commence. 


Fig. 13. 
The dislocations are quite widely spaced along this slip-zone which runs 
vertically upwards to the left of centre. Most of the layers to the right , 
of this line run between pairs of dislocations of opposite hand. 
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Puate L. 
Fig. 14. ‘ 


The black band at the bottom marks the edge of the face. The inactive group 
lies centrally in the lower two-thirds of this area. 


Fig. 16. 

The closely spaced line of dislocations runs vertically upwards on the left of 
this area; individual dislocations cannot be distinguished at this 
magnification. The edges at an angle to the c-axis form at the junction 
of the two systems a little to the left and down from the centre of the 
figure. 

PuatTe LI. 
Figs Te. 

‘The orientation is, as usual, given by straight edges which are barely resolved 
in the upper third of this area. The system which is being eliminated 
lies on the left and the rapidly growing layers are seen on the right. 
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CXXXIV. The Photo-electric Disintegration of Three- and Four-Particle 
Nuclei * 


By J. C. Gunn and J. Invina, 
Department of Natural Philosophy, University of Glasgow*. 


[Received June 12, 1951.] 


ABSTRACT. 


The two-particle and complete photo-electric disintegration cross- 
sections for three- and four-particle nuclei are calculated, using Gaussian 
and new exponential wave-functions. By suitable adjustment of a 
scale constant photo-disintegration, curves can be found with the 
maximum cross-section at any chosen energy above the threshold. For 
a given position of the maximum the Gaussian wave functions require 
a bigger nucleus than the exponential, and then give a larger maximum 
cross-section. The theoretical predictions are concerned with the small 
amount of experimental data available, mainly from the inverse processes 
of p-D and p-T capture. The exponential wave functions appear to fit 
more reasonably with this evidence. 


§1. INTRODUCTION. 


Since the detailed study of the deuteron and two-nucleon scattering 
data has not yet yielded a satisfactory nuclear interaction, it is of interest 
to investigate further nuclear problems involving more than two nucleons 
so that the adequacy or inadequacy of the proposed two-body interaction 
for such problems may be demonstrated. The n-D and p-D scattering 
problems have already been investigated by Buckingham, Burhop and 
Massey (1941, 1948) with various forms of the two-body interaction, 
using the resonating group structure method, first introduced by Wheeler. 
Comparison with experiment should then indicate the best form of the 
interaction. These authors also deduced the photo-electric disintegration 
cross-section of tritium into a neutron and deuteron. Their method, 
however, while neglecting the polarization of the deuteron, involves 
a great deal of computation. The object of the present investigation 
was to see how far an account could be given of the three- and four-body 
photo-disintegration problems with the use of simple analytical wave 
functions. Verde (1950) has treated the photo-disintegration of the 
three-particle system with Gaussian wave functions in an elementary 
manner. While the present work was in progress there has also appeared 
a discussion of the He* photo-disintegration, using Gaussian wave functions, 
by Flowers and Mandl (1951). 
SER. 7, VOL. 42, NO. 335.—DEC. 1951 AY 
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In this paper the photo-electric disintegration of both three- and 
four-particle systems is considered, attention being limited to electric 
dipole transitions. The treatment is crude in the sense that the final 
states of the systems are taken to be plane waves. Equally, for the 
initial tritium and helium nuclei, only the space symmetric part of the 
wave function is considered. The calculations, however, indicate the 
great sensitivity of the cross-sections to the value of the parameters 
occurring in the wave functions of the initial states. Cross-sections for 
inverse processes are calculated by the method of detailed balance. 

Calculations were first carried out with Gaussian wave functions at 
about the same time as those of Verde. Since the Gaussian wave function 
has an incorrect asymptotic behaviour, it seemed worthwhile to investigate 
the effect of a wave function with a better behaviour in this respect. 
The wave functions worked out by Svartholm (1945), using the variation- 
iteration method, were found to be rather complicated to handle in 
photo-disintegration calculations, as also were wave functions of the ~ 
form exp {—a«(r,.+733+ .. .)} which have been used by a number of 
authors in binding energy calculations, for example, Frohlich et al. (1947). 
The new wave functions, used by Irving (1951) for calculations of the 
binding energies of three- and four-particle systems, however, were found 
to lead to fairly simple analysis in the calculation of photo-electric 
disintegration cross-sections. Detailed calculations have been carried 
out for these wave functions and the results compared with those for 
Gaussian wave functions. The theoretical results are compared with the 
experimental data from p-D capture by Fowler et al. (1949) and p-T capture 
by Argo et al. (1950)*. If the parameters appearing in the wave functions 
are determined by variational calculations of the binding energies of 
three- and four-particle systems, very low values for the cross-sections 
are obtained. One reason for this is that the variational calculations 
of the binding energies with a central two-body interaction are incorrect 
and lead to values of the variation parameter which, in the case of the 
new wave functions, concentrate the nucleus too much. Calculations of 
the binding energies, when a mixture of central and tensor forces is used, 
indicate that more reasonable values of the variation parameters are 
then obtained. 

It is to be hoped that when detailed experimental evidence of p-D and 
p-T capture cross-sections—or, alternatively, the inverse processes—is 
available, the theoretical formulae arising from the exponential wave 
functions may be found to give the correct variation of the cross-section 
with energy. In fact, if the parameter can then be chosen to give the 
correct magnitude for the cross-section, thus determining the ‘ size” 


* Experiments with the Glasgow 30 MeV. Synchrotron have demonstrated 
the y—-p process with He*, but are not yet complete (private communication 
from J. R. Atkinson). Preliminary measurements with the G.E.C. betatron 
(Gaertner and Yeater) on the same process indicate an integrated yield 
of 4x 10-°® cm.? MeV. with mean quantum energy about 27 MeV. We are 
very grateful to these authors for a report of their results. 
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of the nucleus, the value obtained should agree with the value calculated 
by a variational method which takes into account non-central forces. 
A characteristic difference between the exponential and Gaussian wave 
functions may be mentioned here. This is that the exponential wave 
functions, for a given “size” of nucleus, lead to an earlier maximum 
cross-section, than the Gaussian. For a given position of the maximum 
above the threshold the magnitude of the cross-section is SES with the 
exponential wave functions. 


§2. GENERAL FORMALISM. 


As we intend to deal in this paper with the photo-disintegration of 
both three- and four-particle systems, it seems worthwhile first to 
summarize the general formalism applicable to the photo-disintegration 
of a system of A particles of equal mass. It is in the first place convenient 
to separate the motion of the mass-centre by the change to the 
coordinates 


al 


A 
R=; “2 r,;, eee tite su ( 1) 
j= 


A; 
rj=rj;—-r,, jAl 

where r; denotes the coordinates of the A particles. In the problem 

chiefly considered in this paper one of the particles only is ejected, say, 

particle 1 ; a suitable set of coordinates is then 


, 1 = } 
Cie een ST) ah | 
; lire Se iy ee Ad mean) 
ESESSIE Ree" INC | i 
Pot, ee ao 4, Ae J 


It may easily be verified that the Jacobian of each of these 
transformations has unit modulus. The second set of coordinates 
corresponds physically with the mass-centre and relative coordinates of 
the residual nucleus of (A—1) nucleons, together with the coordinates 
relative to this mass-centre of the ejected nucleon. 

The wave function of the initial motion may be taken in the form 


Y=u(r1;) exp t(Po . R) 
and similarly for the final state 
YP ,=U,(fon) EXP Up . r,;+P’ . R’). 

These wave functions will in general have to be given the right spin, 
isotopic spin and symmetry properties. It is most simple, however, to 
consider these for individual cases. We shall, in this investigation, 
limit consideration to electric dipole transitions for which the matrix 
element reduces to the form 


A 
[eon exp {—i(p . edd X'7(r;—R) . E,| U(113) 101 om: 


AYE? 
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Here 7, is the usual isotopic variable, Ey represents the electric field vector 
at, say, the mass centre of the system, and the integral corresponding 
with overall momentum conservation has been removed. The probability 
per unit time of disintegration into a group of final states # by a quantum 
of energy k, with its polarization vector along the z-axis, is 


2 
p(Ey). 
(3) 


Here p(E,) is the density of the final states per unit energy interval, and 
the necessary averaging over the initial nucleon spin states, and sum 
over the final spin states must be carried out. Natural units =c=1 are 
employed. 


w=472ke? 


. A 
| #0 on) exp (—1p - @;) {2 7,(tj;— R).| U(1;) 121 Iam 


§3. CHoice oF WAVE FUNCTIONS. 


For analytical simplicity Gaussian wave functions immediately suggest 
themselves. Their bad asymptotic behaviour, however, seems likely to 
lead to an underestimate of the photo-disintegration cross-sections. 
These, however, have been worked out using wave functions of the space 


symmetric form 
pNP exp(—p? 27g). ee ee 
w>J 

In order to obtain an improvement on the Gaussian wave functions 
it seemed desirable to find functions, which, whilst retaining reasonable 
analytical simplicity would have a more correct asymptotic behaviour. 
The wave functions developed by Svartholm were considered for this 
purpose, but in fact greater analytical simplicity is given by the general 
form 

p=Nt exp {—p( 21y)t}/( 27g)". ee ee 
>) wy 

An extensive series of variational calculations has been carried out 
for wave functions of this form, and has been separately described by 
one of us (Irving 1951). 

With both types of wave function a choice has to be made of the 
appropriate value of » for the deuteron, tritium and helium ground 
states, denoted subsequently by pp, wy and w,. This choice is discussed 
later when the photo-disintegration cross-sections are derived. 


§4. PHoro-DISINTEGRATION OF H® anp He® (Two-ParTIcue). 


For electric dipole transitions the probability per unit time of 
disintegration into a group of final states F is 


A 2 

i=1 
where the initial quantum is regarded as polarized in the z-direction. 
The initial state O, say, for the tritium case, is the 28 ground state of 
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H?. The final state / contains one neutron moving as a plane wave with 
momentum p, the resultant deuteron being in either the °(c) a4, or the 
3(c)) state. Carrying out the elementary spin summations leads to 
the result 


2 
p(Hy) . . (6) 


or, for the differential cross-section for the absorption of a quantum 
k with the neutron ejected into a solid angle dQ 


2 
w= 37 ke® 


[ur Piz CXp {—wUp. 01)} u; de, dros 


2 


v4 


e2 ; ; 
do= 73- kpM| | uf 1. exp {—iP .@,)}u dp, dey] dQ... (7) 


Here u, , u, are the initial tritium, and final deuteron wave functions 
expressed as functions of ep, and r,;. For the case of He® the result is 
exactly similar except that allowance must be made for the Coulomb force 
between the outgoing proton and deuteron. 


(a) Photo-disintegration of Tritium. 
Using the Gaussian wave functions defined in equation (4) the differential 


cross-section for disintegration of H* with ejection of a neutron of 
momentum /p is 


: 3 ME 
—= LENE 2 LD aber 
do=|(Z) e?kp?M cos 5G Fey exp ( 4 da. (8) 


where E is the kinetic energy available for the disintegration so that: 
p?=$ME, with M the nucleon mass. The total cross-section o is then 
given by 
Men, ae ME | 

57 aa cen ae exp(—F3). Ss 1: 
Suitable values of x, are found from variational calculations with different 
Gaussian wells, 1/u,—4:2 10-18 em. corresponding with the results of 
Margenau and Warren (1937), and 1/u,.—3-65 x 10-3 cm. with the more 
recent force constants (Rosenfeld 1948) 

Le=1-9X 10-8 em., 3J=43-7 MeV., J /SJO=0-62. 

This second value gives a Coulomb energy of about 0-75 MeV. for He’. 
The cross-section is plotted in fig. 1 (a) with those values of p., and 
with pp/p_p=1-12. 

The calculations were then repeated with the new form of wave 
functions. For simplicity in the evaluation of the integral (6) these were 
taken in the form 
exp {—u,(rig+7%s+731)3} _ Npexp {=v 2p0(p? +P )ee +r) 


ta A arth) en (10) 
with e=r,—3(t2+Fs) | 
; (11) 
and i oe (r3—Yo) | 
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Figs. 1 (a) and 1 (4). 


(1) 


& 
(om 
(in 107m?) 


6 


1) 


@) 


(3) 


o 
n 


10 1S 20 E(Mev) *> 
Cross-sections for two particle photo-disintegration of H® 
(a) for Gaussian wave functions with the values 
(1) 1/pa=4-2 x 10-38 cm., (2) 1/up=3-65 x 10-13 em. 
(b) for exponential wave functions with 
(1) Ijup=3-0 x 10-18 em., (2) 1/un=2-5 x 10-38 cm., (3) 1/up=2-0 x 10-8 em. 
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The variational calculations of binding energy then show that 
1/up=2-0X 10-3 cm. is a suitable value of the parameter, corresponding 
to a binding energy of 7-1 MeV. for a standard Yukawa potential 


1/e=1-:17X10- cm., 3J=67-3 MeV., qg=0-69 . . (12) 


(of. Irving (1951) ). If, rather than using the best variational value, 
Hy is chosen to fit the Coulomb energy of He, then 1/y,=2-5 x 10-43 em., 
corresponding with a Coulomb energy of 0:8 MeV. (binding energy of 
H?=—6-3 MeV.) is suggested. In the curves of fig. 1b the rather extreme 
case of 1/u,»=3-0 x 10-8 em. is also shown. The binding energy with the 
potential (12) is then 5 MeV. 

For the deuteron we take 


nee Nas pe] (ar). een eee (13) 


with r as defined in (11) and y determined by the deuteron binding 
energy Wp, y=2(MW),/3)?. The integrals involved are long but quite 
elementary and lead to the result for the total cross-section 


9 KES? yf 
o= 25¢2 wae We {f(A) . ar) ea . : . (14) 
J 7TA—2 15 15 “asa 
wits 70 ary =I t G=TTIO (aA) 


and A=(p?+ 2u%)/y?. This cross-section is shown plotted as a function of 
photon energy in fig. 1b. 


(b) Photo-disintegration of He?. 

In this case allowance must be made for the Coulomb interaction 
between the outgoing proton and deuteron; otherwise the general 
formulae are as for H?. In view of the approximate nature of the 
calculation, it is felt that this allowance can be adequately made by 
multiplying the H? form of cross-section by a barrier penetration factor 

(Bethe 1937), at the same time taking note the of different binding 
| energy of He? (7-74 MeV.), which brings the threshold to k=5-5 MeV. 
Detailed curves are not drawn for this case but the p-wave barrier 
penetration factors for a few energies are shown in Table I. 


TaBLe I, 


E (MeV.) 0-5 1-0 1-5 2-0 2:5 
exp (—P,,) (R=4:0 x 10-¥ cm.) 0-02 0-11 0-27 0-49 0-81 
exp (—P,,) (R=5-0 x 10-* cm.) 0-05 0-25 0-64 1-0 — 
E=energy in mass centre system. 
R= barrier range. 


Reviewing the results, we have found for the two-particle disintegrations 
it must be admitted that neither the results for the Gaussian nor for the 
exponential wave functions can be very accurate. Considering for the 
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moment, for comparison, the case of photo-disintegration of the deuteron, 

-it is known that a fairly good result can be obtained with the use of an 
asymptotic form for the initial wave functions and plane waves for the 
final neutron and proton. This possibility arises from the fact that the 
neutron and proton forming a deuteron exist for the greater part of the 
time outside the range of the nuclear forces. For the three- and four- 
particle nuclei this effect becomes successively weaker, and the central 
portion of the dipole matrix element integral becomes more important. 
Allowance for distortion of the outgoing waves is difficult to make, but 
it will probably tend in the direction of increasing the cross-sections at 
low energies and bringing the maxima nearer to low energies. In favour 
of the exponential wave functions it can at least be stated that they have 
a fairly correct asymptotic form, and should give the general shape of the: 
cross-section better than the Gaussian form. (It may be recalled that 
even for the case of photo-disintegration of the deuteron Gaussian wave 
functions give a rather poor result.) . 

Unfortunately, only very little experimental information is available 
with which to compare the theoretical predictions. It is limited to the 
p-D capture cross-section reported by Fowler et al. (1949). These 
measurements indicate an approximate cross-section 


d= 0-74 05 =< 102 om 29s eee 


for the H*(p , y) He* capture process in the region 0:5—1-5 MeV. The 
capture cross-section for this process is given in terms of the corresponding 
inverse photo-disintegration by 


o(E,)=(k2/2MB)o,(k). . . . . . . (16) 


Here E,, is the proton bombarding energy in the laboratory system 
(deuteron stationary), and related to the earlier E (energy in the mass 
centre system) by E,—3K. The corresponding quantum energy / is 
given by k=E-+B, where in this case B=5-5 MeV. From (16) the 
theoretical capture cross-sections may be derived using the earlier (9) and 
(14). For the present we merely note that for small energies o,cE}. It 
has, thus, neglecting Coulomb barrier effects, a very swift rise above 
threshold, for any type of wave function or choice of parameters. The 
difference between different theoretical models is only observable when 
measurements beyond the initial rise are available. This is by no means 
true in the present case. Indeed, the experimental results given in (15) 
apply only up to E,,=1-5 MeV., when the Coulomb barrier effects are still 
appreciable, and can easily be fitted, at least in shape, with either type 
of wave function using the barrier penetration factors of Table I. The 
most significant use to make of (15) seems to be to derive the magnitude 
of the corresponding photo-disintegration cross-section. We find that 
this should have the value 4x 10-%8em.2, for E=1 MeV. The 
corresponding cross-section without barrier should be at _ least 
1610-8 cm.? This should be approximately equal to the cross-section 
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for the H*(y, n)H? process at 1 MeV. above the threshold. Comparison: 
with figs. 1 (a) and (6) shows that none of the theoretical cross-sections has 
such a rapid rise above the threshold. Even the exponential wave 
functions with I/u,=3-010- cm. gives only a cross-section of 
5x 10-** cm.” In this energy range the distortion of the outgoing wave 
may be important and give rise to the discrepancy noted. 


§5. PHoTo-DISINTEGRATION oF H? anp He? (THREE-PARTICLE). 


As most of the binding energy of H? or He?® is lost when the dis- 
integrations H?>n+D, He®+p+D are carried out, we might expect, 
in view of the greater phase space available, that, except just near the 
threshold, the three-particle disintegrations, of which H?>+n-+-n+p is 
typical, should predominate. Calculations along the lines of the preceding 
ones confirm this belief. They have been carried out both for the Gaussian 
and for the new exponential form of wave functions, once more using the 
approximation of plane waves for the outgoing particles. With the same 
notation as previously the two cross-sections are 


Gaussian 
ake? M 
c= Fp ME? exp (— 5). > ee es ere Op | 
Exponential 
1507?e? (uy\4 kERM? 
eae (4) (+38) ees ys (LS) 
2M 


For the Gaussian wave functions the maximum cross-section comes 
approximately when 


1 2 
E=500 Mev. /{ —in 10-3 em. } 
MT 


and the magnitude of the cross-section varies approximately as Lys. In 
Table II. these maximum cross-sections are shown for various values 
of 1/u;, for the disintegration of H?. 


TaBLe II. 


Maximum cross-sections for H’+n-+-n-+p 
(Gaussian wave functions). 


L/pr (units 10-13 em.) 3-65 4.2 5-48 
Emax (MeV.) 35 28-4 17 
mas (units 10-2? em?) 2:2 2:9 4-8 


For the exponential type wave functions rather similar values of 
cross-sections are obtained ; for reasonable values of 1/1;, however, the 
maxima of the cross-sections again occur at lower energies. The 
maximum cross-section is approximately proportional to 1/u2, and the 
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corresponding energy to »%. In fig. 2 the results are shown with 
1/uy=25X10-% cm. It will be noted how these values dwarf the 
two-particle disintegration cross-sections of § 4, except at the threshold 
where only the two-particle process is possible. This will add to the 
already considerable experimental difficulties in the way of observation 
of the He*(y, p)H?, H3(y, n)H? processes. 

We have earlier noted the dangers in the use of the Born approximation, ° 
i.e. plane waves in the final state, in excitation problems of this sort. 
The method is to some extent consistent in that, for the initial state 
also, wave functions are used which have only asymptotic validity. 


Fig. 2. 


in 2 
(in 10°2%cm?) 


c 5 a 1S E(revy 20 


Three particle photo-disintegration cross-section of H3, for exponential 
wave function, 1/uy=2-5 x 10-8 em. 


§ 6. T'wo-ParTicLe PHoTO-DISINTEGRATION OF «-PARTICLE. 


Here we consider only the transition He*(y, p) H3, neglecting at first 
the Coulomb interaction between the outgoing proton and triton. This 
is, in any case, negligible for quanta above about 22 MeV. From the 
general formulae given in § 2. wé have 


2 


Bek 
w=nthe? |utp,. exp (—ip . P,) 4. de, dts; dra pike eae ee 


The initial state is the 4S ground state of He?, and the final state F contains 


one proton moving as a plane wave with momentum p, the resultant 
triton being in a 2S state. 
Introducing the coordinates 


U=—F,+3(fe+r3+44) 
V=4(f3+4r4)—}(re+r3+r,) , - +. > (20 


w=r,—-r, 
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we have 


4 
aa 2 ry Bur + 2407+ 20? | 
i — 
. (21) 
4X r= 3(240%+ 2w?). 
i9=1 
Gaussian wave functions : The calculations are elementary leading to 
: : ME 
o=3 x 4473/2e2yoia__Pr__ ps2 ex {- ae Weer tS we 22 
Mal Span 4p) PL a2 a 


The parameters , and wy, may be determined by variational calculations 
of the binding energies of the «-particle and triton, assuming, for 
example, a Gaussian well. The nuclear parameters 


1fe=1:9X10-¥ cm., 3J=43-7 MeV., q=—0-62 


then give 1/u,1/u,—3-65x10-% cm. yielding binding energies of 
approximately 4 MeV. and 23 MeV. for the triton and «-particle respec- 
tively. The Coulomb energy of He* is then approximately 0-75 MeV. 
The nuclear parameters 1/k=2-25X10-! cm., 3J(=35-6 MeV., q=0-6 
determined by Margenau and Warren (1937) to fit the triton binding 
and deuteron data give slightly too small a value for the «-particle 
binding energy. For this case 1/y,1/uy»=4:2x 10-3 em. In fig. 3 (a) 
o is graphed against E for 1/u,—1/uy, several values of the parameters 
being taken to show the dependence on them. = 1/u,=1/up~=5-25 x 10-% cm, 
has already been considered by Flowers and Mandl (1951). 
The calculations were repeated with exponential wave functions of 
the form 


ty=NiP exp {Hal 2 ri) 3/2 ri} 
i>j i>j 
N12 exp {—p,(3u?-+ 24v?+ 2w*)?} (23) 
era ee (S12 4k OF) re 
and p= NA? exp {mon t13s +151)? hy (rie 13+)? 
__ NiPexp {—F/3ynl24e?+2u%)1) 
man S/n 
The cross-section is then given by 
7X 444 
PimeSea 2 


where f@= | (735-+846524)/8+ 21/(1-+2) 42/48) 


(24) 


ue 
emia Ti MEM em {fey}... (25) 


3 cot1z 


(115524+ 63022435) | 


i 
and C= Cee CC 
ra 7. T 
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The graph of o against E is shown in fig. 3 (b) for different values of the 
parameters. 1/4~=2-5x10-% cm. gives the correct Coulomb energy of 
He?. 1/u,—1-7 10-8 em. gives a binding energy of 28 MeV. for Het; 
Iu, =2-0 x 10748 em. a binding energy of 20-5 MeV.; 1/u,=2-5 x 107 cm., 
a binding energy of 12-5 MeV. The rapid rise of the cross-section is 
evident. The exponential wave functions give earlier maxima than the 


Fig. 3 (a). 
8 


Co 
(10-27 cm2) 


7, 


(e) 10 20 ‘30 40 50 €(MeV) 


Cross-sections for two particle photo-disintegration of He# 


(a) Gaussian wave functions with 1/u.=1/uy having values 
(1) 6:10 x 10-48 em., (2) 5-25 x 10-18 em., (3) 4:20 x 10-33 em. 


Gaussian, the values of the parameters being smaller as well; that is, 
a more concentrated nucleus is found for cross-sections of comparable 
magnitude, 


The cross-section for the inverse capture process can again be calculated 
by the method of detailed balance yielding 


o({E,)=(k2/3ME)a,(k). . . . . . . (26): 
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Here, as previously, E,, denotes the proton bombarding energy in the 
laboratory system, E(=}?E,,) is the corresponding energy in the centre of 
mass system, and o,(k) denotes the photo-disintegration cross-section 
for quantum energy k. To illustrate the general form of the theoretical 
capture curves, the capture cross-section is given in fig. 4 for the 
exponential form of wave function with 1/u.—=1/u,—2-5 x 10-8 em., and 
for the Gaussian form of wave function with 1/u»=1/1,=6-10X 10-13 em. 


Fig. 3 (0). 


8 
(on 
{10-?’cm?) 


7 


10) 10 20 30 40 50 (MeV) 


Cross-sections for two particle photo-disintegration of He* 
i i —13 
b) Exponential wave functions with (1) 1/u,=1/u9=2:5 x 10-® cm., 
(2) i Bae 10-23 cm., Iup=2:5 x 10-8 em., (3) I/jue=1-7 x 10-* cm., 
Luy=2-5 X 10-18 cm. 


The experimental evidence available relates mainly to the capture 
process, coming from the measurements of Argo et al. (1950). BOOS 
has only been carried out up to a bombarding energy H,=2°6 MeV., 
which is just approaching the top of the barrier for p-wave protons. | It 
is clear that the shape of the experimental curve can, as noted by Flowers 


1366 J.C. Gunn and J. Irving on the Photo-electric 


and Mandl (1951), be explained by a combination of the initial shape 
without barrier o,E,!, and suitable barrier penetration factors, more 
particularly in view of the rather arbitrary calculations involved in the 
latter. Once again, however, the magnitude of the capture cross-section 
is interesting. This is quoted as approximately 10~°” em.” at E,,=2-6 MeV. 
The corresponding photo-disintegration cross-section is about 10-6 em.* 
at k—2 MeV. None of the theoretical curves predicts so sharp a rise in 
the cross-section above the threshold, though once again the exponential 
wave functions give a better agreement. If the experimental evidence . 
regarding the photo-disintegration process mentioned in a footnote in 
§ 1 is confirmed, then support is provided for a cross-section approximately 


Fig. 4. 
(@) 
4 
(i0-28cm2) 
3 (2) 
2 
| 
md 5 10 15 


20 25 
E, (MeV) 
Cross-section for proton capture by tritium, neglecting Coulomb barrier 


(1) with exponential wave function, ]/jup=1/u,=2-5 x 10-8 em., (2) with 
Gaussian wave function, 1/up=1/u,—6-10 x 10-3 em. 


of the form given by the exponential wave functions with — 
I/pyp=1/u,=2-5x 10-8 cm. These measurements give an integrated 
yield of 4 10-°° cm.? MeV. for the cross-section, with about 40 per cent 
of the events below k=24 MeV., which could not be reconciled with any 
of the curves based on Gaussian wave functions. 

In conclusion then it appears that a reasonably good fit to the photo- 
disintegration and inverse processes for the three- and four-particle 
nuclei may be hoped for with the exponential wave functions. They are 
at least superior to the Gaussian functions in predicting the early rise 
of the cross-sections above the threshold. There is at present, however, 
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too much uncertainty regarding the choice to be made of the parameters, 
Hy and p,. The binding energy calculations with central forces (Irving 
1951) are certainly inadequate for this purpose, but the hope is that 
their extension to include non-central forces may yield a consistent 
picture. We should then, in the exponential wave functions, have an 
analytical model of reasonable simplicity for three- and four-particle 
nuclei, with the possibility of extension to other light nuclei. 


Complete disintegration of Het. 
Using the same approximations as in the case of the three-particle 


system, the cross-section for the process 


He*+y>n+n+p+p 
is easily calculated. 
For the Gaussian wave functions, we have 


1 7/2 Mur 
2 


ME 
2) SE eee AP By cee 7), 
o= 945 aa © ae kE exp me ee es OA 0 


The values of o,,,, for several values of the parameter are then as follows : 


pg (lO 2-em.:) E (MeV.) Cua k 27 cm") 
3:6 66-4 3-04 
4-2 48-0 4-58 
5-25 30-0 8-66 
6-1 22-0 13-7 


For the exponential wave function, 


22/2 , 1 pike? 
“45° MP2 (B+ 2u2/M) 


yielding the following results for o,,,. : 


Lt, (10-13 cm.) E (MeV.) Cmax (10-2? em.2) 
1:7 28-7 3-7 
2-0 20-0 6:2 


2-5 12-2 12-8 


The aan for the exponential wave function again occur at much 
lower energies than those for the Gaussian wave function. It would be 
interesting if experimental information could be obtained on this point. 
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[Plate LIIT.] 


ABSTRACT. 


The recrystallization texture of high purity (99-995 per cent) cold 
drawn aluminium wire has been found to be double, viz., [322], [100]. 
The [100] component is relatively weak. The scatter in the major [322] 
component increases, and the amount of the [100] component decreases, 
the higher the temperature of anneal. 


RECENTLY it has been suggested (Hibbard 1950) that the double [111], 
[100] deformation texture of drawn face-centred cubic metal wires 
previously described (e.g. Schmid and Wassermann 1927) is actually 
an intermediate texture which changes to a single [111] texture with 
sufficient deformation. The [111] texture has been predicted on theoretical 
grounds (Hibbard and Yen 1948, Calnan and Clews 1950). Observations 
Hibbard made on drawn wires of six metals showed that for a large 
deformation (98 per cent reduction in diameter) all the wires had 
essentially a single [111] texture. In the case of aluminium, however, 
there was a minor [100] component, which Hibbard ascribed to room 
temperature recrystallization. His reasons for this were that aluminium 
of the purity used (99-994 per cent) has been known to recrystallize at 
room temperature ; that the [100] component on the X-ray photograph 
was characteristic of diffraction from large recrystallized strain-free 
grains (2. e. “‘ spotty” rather than diffuse); that the [100] component 
increased in intensity with increasing deformation above 90 per cent ; 
and that the deformation texture of less pure drawn aluminium wire 
has been reported previously as single [111]. 

It is interesting to observe that v. Gdler and Sachs (1927) working 
with less pure aluminium (99-93 per cent) did find evidence of a [100] 
recrystallization texture. The wire they used had been reduced in 
diameter by 80-6 per cent and the deformation texture was [111] with 


* Communicated by the Authors. 
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a trace of [100]. An annealing temperature of at least 200°C. was: 
required before recrystallization was evident as faint spots on the 
Debye-Scherrer rings of the X-ray photographs. After complete 
recrystallization, requiring temperatures of at least 375° C., the [111] 
deformation texture was observed to be preserved and an appreciable 
amount of [100] texture was also present. 

Recent experiments in this department have also shown the existence 
of a minor [100] component in the recrystallization texture of high purity 
aluminium wire. The actual texture is double, viz., [322], [100]. The 
majority of grains have the [322] texture, the minor [100] component 
being relatively weak. 

The aluminium, kindly supplied by the British Aluminium Co. Ltd., 
was of 99-995 per cent purity and had received a 95-2 per cent reduction 
in diameter. Plate LIII. (a) is a flat film X-ray photograph showing the 
deformation texture, taken with partially filtered MoK« radiation with 
the beam perpendicular to the vertical wire axis. It shows the [111} 
texture on a background of complete rings from the randomly oriented 
grains. On the first ring there is a weak reflection at A (Plate LIII. (a)) 
which can be accounted for if a small percentage of the grains have the [100] 
direction along the wire axis. In this case a slight intensity maximum 
would be expected at B on the second ring, and this can be seen on the 
original negative. A weak [100] component then exists in addition to 
the strong [111] component and the small percentage of grains of random. 
orientation. The reflections at A due to the [100] component are not 
spotty, 7.e. not characteristic of large recrystallized strain-free grains. 
The wire was annealed at various temperatures and in all cases a 
new texture appeared. Plate LIII.(b) is a typical photograph. Here 
the wire had been annealed at 275° C. for two hours and a comparison 
with Plate LIII. (a) shows the change in texture. At regions A and B. 
(Plate LITI.(b)) on the first and second rings respectively the fairly weak 
reflections are at the positions predicted by a [100] texture. For these 
reflections the maximum angular deviation of the [100] direction from the 
axis of the wire is about 5°. For the remainder of the reflections an 
analysis has shown that the texture can be accounted for almost exactly 
with the [322] direction along the wire axis. The randomly oriented grains 
disappear. Anneals were given from approximately 150° C. to 550° C. 
and this new texture was always found. It is due to primary 
recrystallization ; subsequent grain growth did not affect the texture. 
The amount of the [100] component decreased, and the scatter in the 
[322] component increased, the higher the temperature of anneal. In 
Plate LIT. (6) the maximum angular deviation of the [322] direction 
from the axis of the wire is between 2° and 3°. 

v. Géler and Sachs (1927) used aluminium of three different purities 
and found the minor [100] recrystallization component strongest for the 
purest metal. 
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Their results were as follows :-— 


(i.) a preservation of the [111] deformation texture and the appearance 
of a [100] component with recrystallization of 99-93 per cent aluminium 
(previously mentioned) ; 

(u.) a partial preservation of the [111] deformation texture and a 
suggestion of a [100] component with recrystallization of 99-0 per cent 
aluminium ; and 

(ii.) a partial preservation of the [111] deformation texture with 
recrystallization below 500°C. but increasing randomness with re- 
crystallization above this temperature and no [100] recrystallization 
component for 98-7 per cent aluminium. 

Other observers (Schmid and Wassermann 1928) found no evidence of 
a [100] recrystallization component for 99-95 per cent aluminium. 

Nevertheless, the results of v. Géler and Sachs suggest that the amount 
of the [100] recrystallization component depends on purity, and it may be 
significant that with the higher purity material in the present work 
temperatures lower than the 375° C. of v. Géler and Sachs produced this: 
component. This seems to support Hibbard’s suggestion that the room 
temperature [100] component in his aluminium wire is due to re- 
crystallization. In the particular case observed here, however, the 
existence of the minor [100] recrystallization component does not 
definitely establish that the room temperature component is due to 
recrystallization. Two possibilities exist :— 

(i.) that it is due to deformation only and that recrystallization 
preserves it ; and 

(ii.) that it is already present due to recrystallization and higher 
temperature anneals merely retain it. 

The fact that the room temperature [100] component here was not 
“ spotty ” and that temperatures of at least 150° C. were required before 
the X-ray rings did become “‘ spotty ” tends to discount (ii.) while some 
of Hibbard’s previously mentioned observations would seem to discount 
the first possibility. It is evident that further investigations are necessary 
to determine the true origin of this room temperature [100] component 
in aluminium wires after large deformations. 

The major [322] recrystallization component has not hitherto been 
reported. It may be characteristic of this wire only, as purity, prior 
deformation and mode of recrystallization are determining factors for 
the resulting texture. In many cases in face-centred cubic metals a 
certain relative orientation, corresponding to a rotation of 30° to 40° 
about the [111] axis, has been found between a highly oriented cold 
worked matrix and strain-free grains growing in this matrix (Beck and 
Hu 1949). The same relative orientation was found with coarsening or 
secondary recrystallization. The retention of the [111] drawing texture 
on annealing, previously observed, would be in keeping with this 
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orientation relationship, but the production of the [322] texture on 
annealing from the [111] drawing texture, observed here, does not agree 
with this relationship. ; 

With increasing purity of available metals it is more necessary than 
ever, when stating recrystallization textures, to. quote the purity, 
temperature of anneal, mode of recrystallization and the magnitude of 


the previous deformation. 
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ABSTRACT. 


The quantum mechanics of a system of identical interacting particles 
must lead to the classical hydrodynamic equations of motion at high 
temperatures, because of the correspondence principle. On the other 
hand, the behaviour of helium IT shows that this is not always the case at 
low temperatures. In this paper it is shown that in certain cases the 
quantum description requires an extra parameter, which is the potential of 
a new velocity field superimposed on the classical motion. Expressed in 
semi-classical terms, the condition for the existence of this new parameter 
is that the probability, in the equilibrium state, of a particle having a 
very large de Broglie wavelength (that is, a negligible momentum) is 
finite. This condition is satisfied in one model of a superfluid system, a 
condensed Bose-Einstein gas, but not ina crystal. <A tentative theoretical 
interpretation of two basic equations of the empirical two-fluid theory of 
helium IT is given, in which this new parameter determines the velocity 
of the superfluid. 


§ 1. INTRODUCTION. 


BorN AND GREEN (1947, 1948) described a method of studying transport 
processes in a system of identical interacting particles obeying quantum 
mechanics by using the analogy with the corresponding classical system. 
Gurov (1948, 1950) carried out a similar calculation. In both cases the 
conclusion was that the equations of hydrodynamics and heat transfer 
were identical in form with the classical equations. These treatments 
assumed that the three quantities describing a “‘ normal” state of the 
system (Born and Green 1947, Chapman and Cowling 1939) were the same 
as for a classical system : the density, velocity and temperature fields. In 
other words, if at one instant the system is in a state described by given 
fields of these three parameters, its state at a later time is sufficiently 
described by the new fields of the parameters resulting from the thermo- 
hydrodynamical equations of motion. From the correspondence principle 
it is clear that the three parameters of classical mechanics are sufficient in 
quantum mechanics at high temperatures, but, as we shall see, this is not 


always true at low temperatures. 
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§2. PROPERTIES OF THE REDUCED Density MATRICES. 


Suppose the system consists of N identical particles of mass m inter- 
acting by central forces whose potential is U, and let the Schrodinger 
representative of the density matrix of the corresponding Gibbs ensemble 
{Dirac 1947) be 

D(qi,- ++ v3 Gi + - Ay)» 
where q;,...-@Qy are the position vectors of the N particles, and D is 
Hermitian. The reduced density matrices, which are also Hermitian, are 
defined by 


xf... [Dian .- dns Gi --- Ge Getz--- Gv) 4Gut1---dqy. - (1) 
The equation of motion for D is 
ih OD/ot=HD— DH, al) eee 


where H is the Hamiltonian operator of the system. After setting 
qo=q), -- - Gvy=qy, this may be integrated over qp, . . . qy, giving 
ihOR (qr 3 q1)/Ot= —(h?/2m)(Vi—VP)Ri(qi; a1) + JTU(qi—ae) 
—U(qi—q2)]Ro(qi, G23 7,492)4q2, «+ - (3) 

where V,;=0/0q;, etc. 

By using Wigner’s result (Wigner 1932, Moyal 1949) that a density 
matrix R,(q;; qj) corresponds to a classical probability density F(x, p) 
for the position x and momentum p of one particle such that 


Ri(qi ; 47)~JépF(3[qi+47], p) exp [ep - (qi—qi)/h], - - (A) 
combined with the classical result 
F(x, p)oc exp [—p?/2mkT],  .~.~ <= ne een 


or alternatively by approximately solving the Bloch equation 
oD/o(kT)4=—HD (Husimi 1940), we obtain a high temperature 
approximation 


Ri(qi ; qi) © exp [—(qi—q7)?mkT/2h?], . . . - (6) 
for equilibrium at the temperature T. Therefore 
lim R,(q;; q;)=0. esas ds = ey eee 


1qi—qi’ | > 
At low temperatures (7) may not be satisfied for every type of system. 
To see the physical significance of the failure of (7) let us assume for the 
moment that R,(q;; q;) depends only on |q,—q;| and approaches a 
constant value L when |q,;—qj|—0oo. Then the Fourier inverse of (4) 
yields a divergent integral which may be interpreted as 


F(x, p)=f(p)+ La(p), . - . - - . + (8) 
where /(p) is a regular function. That is to say, the probability of one 
particle having zero momentum is the finite fraction L/R,(x; x). Since (8) 
is in fact an expression for the mean occupation numbers of single-particle 
levels, L must vanish (and (7) therefore be satisfied) in a Fermi-Dirac 
system, where the exclusion principle limits the occupation numbers. 
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33. ASYMproTic BEHAVIOUR OF THE REDUCED Densrry Marrix 
FOR ONE PARTICLE. 
In the general case when (7) does not hold we shall consider the behaviour 
of Rj(q;; q7) for values of |q,—q)| large compared with both A(mkT)-¥2 


and the range of U. The integral in (3) may be written, using the 
Hermitian property of Ro, 


fU(r)[Re(qz autre; gi, qi+r)—R3(qi, qi+r; qy,qitr)|dr. . (9) 


eae that the asymptotic form of R.(q,;, q;+r; qj, q,+r) for large 
1q:—q;| is 


R.(qz, q,;+r ; qi, q,+1r)~A(qz, r)Ri(qz 3 qi), a s a (10) 
just as the neighbour distribution function of a liquid obeys 
Ri(dn 2 Gr Ge)~Ri(qis ai)Ri(qo; qe) . - + (11) 


for large |q,;—q,|. For large | q,—q; | (3) becomes 
th ORi(q1 5 q1)/0t ~ — (h?/2m)(Vi—V7?)Ri(qi 5 qi) 
+{X(q)—X*(q)IRilar; 4), - - (12) 
where X(q)=J A(q, r)U(r) dr=V(q)+7W(q), . . . . (18) 


Since the variables in (12) are separable the general solution is a linear 
combination of expressions of the form 


VAs) a) ee eer te cake sy (14) 
where thoP/dt=— (h/2m)V2F+ XY 2. 2... (15) 
At equilibrium D=e~’", and since the Schrodinger representative of H 
is real, those of D and the R’s are also real. The state of the liquid is 
uniform, and, therefore, for large enough | q,—q;| it may be assumed that 
R,(q1 3 qi) ~ const. 
This is of the form (14) with Y=const. If the density matrix of a non- 
equilibrium state is sufficiently similar to e~*" the required solution of (12) 
will still be a single expression of the form (14) : 
Ry(q13 41) ~ Yq) ¥*(q))- oe eee. (16) 
The example of a condensed Bose gas may clarify the argument of this 
section. Equation (16) is true, with Y the wave-function of the single- 
particle state into which condensation has taken place. Husimi’s formula 
(1940) 
Re(qz 423 41, 42)=Rilar; a1)Bi(qe; 42)+Ri(qi; 42)Ri(qe 5 2) 
is not valid for a condensed Bose gas; for example, at absolute zero, when all 
the particles are in the ground state, 


Ri(qr; 41)= (41) P*(qi) ; 
Ro(qu, M23 A> Iz) = H(i) ¥ (qe) P*(qi) ¥*(qz)- 
The correct formula is 
R2(qa Qe 5 qi 42 )=Ri(q1; 47) Ri(q ; qs) +Bi(qr; q2)Ri(q ; q:) 
—¥(q,)P (qz2)P*(qi)P*(qz) - - - + + (1%) 
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(compare London 1943). Assuming that this formula holds for non- 
equilibrium states, (10) and (11) may be verified, with 
A(q, r)=R,(q+r; qt+r)+R(q; q+r)¥(q+r)/P(q)—P(q+r) P*(q+4). 


Madelung’s transformation (1927) brings equations (15) and (16) into a 
form similar to the hydrodynamic equations of motion with velocity 
potential 4, density | ¥|?, external forces of potential V and a source 
distribution of strength 2| ¥/?W/h : 


Ri(qi; 41) ~Qa(ai; a1) exp {em[d(qi)—P(ai)I/A}, - + (18) 
where Qo(qr: ai)=|¥ (a1) | 1% (a7 
0|¥ P/ot-+V-(|W 2? Vd)=2 | PW/h, st Sieh ke 
0¢/dt+4(Vb)2+V/m=h2 V2 |¥ |/2m?|P [= 0. | 
At equilibrium Q,, depends only on the density and temperature, so 
that even for non-equilibrium states it is a function of the parameters of 
the classical description. On the other hand, the velocity V¢ corresponds 


to no quantity in the classical description; it differs from the classical 
velocity in its laws-of change and in being irrotational. 


$4. Tar Two-FLurp Mops. or Hetivm II. 


The two-fluid model may be interpreted tentatively as follows: Assume | 
that real functions Q,,(q,; q’;) and Q,(q;; q;) can be chosen, with 
Q,, > 0 as |q,;—q;|—, such that if at one instant 


Bi(qi; 47)=Qn(4i 5 Gi) exp {4m(qi—q:) - u(3[q-+47])/h} 
+Q.(qi; a7) exp {imm[$(qz)—$(q2)]/h} 
+O(Vu)+O(VV9) 22) 5S See ee 

(where the velocities u and V¢ do not alter appreciably in a distance 
h(mkT)~"”), then at later instants it is of the same form, with new real 
functions Q,,, Q,, u, and 4. If Q;—0 this would describe a system with 
classical velocity field u (Gurov 1948). The density is 
p(x)=mR(x; x)=mQ,,(x; x)-+mQ,(x; x) 
at Ne daa hese me Gees Sc SD 
and the current density (Born and Green 1948) 
J(x)=—(th/2)(Vi—V)Ri(qu s 7) [qi=qi=x] 

= pPp(X)U(X)-F p(X) Vo( x) on ee 
These two equations are the basis of the two-fluid theory (Landau 1941, 
Tisza 1947) ; in Landau’s version the superfluid velocity is irrotational. 


§5, EXAMPLES. 


Two examples will illustrate the argument. In a crystal the single- 
molecule wave-functions are confined to small regions, so that (7) is true 
and classical ideas should explain its mechanical properties. In a condensed 
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Bose-Einstein gas with negligible repulsive forces and periodic boundary: 
conditions (Bogoliuboy 1947) 


Ri(qi; 4;)=B-1Zn, exp il (q:—qi), | 
1 erate ees (23) 
Q..(q: ; q1)=%/B¥0, | 


where B is the volume of the container and », is the mean occupation 
numbér of the single-particle state whose wave-number is 1, given by 
m=[exp (1?h?/2mkT)—1]-}, [i40], | 
In=N. aber e. (24) 
j 
Bogoliubov claimed that such a system shows superfluidity. 

We conclude that (7) is the condition for the classical equations of heat 
transfer and hydrodynamics to apply, and that when (7) is not satisfied a 
new quantity enters the equations. This quantity is the potential of a 
new velocity field superimposed on the classical motion. 
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ABSTRACT. 


The iteration method of Svartholm (1945) is employed to improve the 
approximate wave function used by Tomonaga (1947) to describe the 
scalar meson field of an isolated nucleon. It is found that this wave- 
function of Tomonaga is surprisingly good, even in the region of 
intermediate coupling, the change in coupling parameter for given energy 
‘being at most 3 per cent after the first iteration. 


$1. 

CoMPARISON of recent experimental data concerning nucleon—nucleon 
interactions and meson production with the results of perturbation theory 
calculations, which depend on the smallness of the meson—nucleon 
interaction, shows that this interaction may not in fact be regarded as 
“weak”. For the opposite limit of ‘‘ strong ” coupling, Wentzel (1940) 
has proposed a semi-classical method, and more recently, Tomonaga 
(1947) has discussed a new approach to the problem, which covers also 
the intermediate region between these limiting cases. In this work of 
Wentzel and Tomonaga it was necessary to cut off the high momentum 
components of the meson field in an arbitrary way in order to obtain 
finite results. Recently Dyson (1951) has shown how by renormaliza- 
tion of mass and field strength a Hamiltonian can be obtained for which 
such difficulties do not occur. Accordingly it is of interest to re-examine 
‘Tomonaga’s method, and it is the purpose of the present work to test 
the validity of the trial wave-functions proposed by Tomonaga, in the 
region of intermediate coupling. 


§2. 


Tomonaga investigates the stationary state of a nucleon coupled with 
a scalar meson field, considering the nucleon as infinitely heavy, so that — 


its recoil and the creation of nucleon pairs are neglected. Here these 
assumptions are retained. 


For charged mesons the Hamiltonian is, in this case, 


H = | oferta DED Mk —L | GU k\agt bE +(at+b,)7 Mk, . (1) 


* Communicated by Professor R. E. Peierls, F.R.S. 
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where a;, a% are annihilation and creation operators for the Fourier 
component k of the positive meson field, b,, b* the corresponding 
operators for negative mesons, and w= /(u2+k?). With G(k) as used 
below, the relation of J to g, the usual coupling constant, is l=g/(47). 
The divergence difficulties associated with the Hamiltonian (1) are 
avoided by including in G(k) a cut-off function F(k), and in the work 
below we shall find it convenient to use the expression 


G(k)=F(k)//o=e?™/ Va, 

M being of the order of the nucleon mass. (In fact we have chosen 
M5.) - ; 

Following Fock (1934) we let functions x, ,,(k,...k,; K,...K,,) 
describe the presence of positive mesons of momenta k,...k,, and 
m negative mesons of momenta K,...K,,. These functions are 
symmetrical for interchanges within (k, ...k,,) and (K, ...K,,,) separately, 
and they vanish unless =m, m-+-1 because of the requirement of charge 


conservation (for the case of a proton). The Schrédinger equation for a 
stationary state of the system, 


Bie nga eee ee eee te (2) 
becomes 


ExXo0= —lfG(K1)x19(k1)d?hy, 
(ee ae 
i n 
(SLB 6K yn ually § Bae Bea Riga Ky) 
$-V(O+1)[ oy erdxnetya(le ++ Busi Kase Ba)Phn af 
(E—2,.0—2 0) Xn 41, n 


if n+1 
= — kk ———_. bee ed Kaeo 
cl Wer 2. GK) Xn n (Ki k,_, kj, n+1 1 ) 


$V(n-+1)| OKyaa)xnen, n+ (Ky eee Ke aoa - KK} 


(3) 


Tomonaga proposed as an approximate solution of (3) the product 
function 


Xn, m( Ki #162 |e K, ee K,,)=Cn, m?(K1) Hide $(k,, )xb(K,) CAFO 2(K,,); (4) 


the functions ¢, % being normalized, and the ¢,,, being chosen to 
minimize the expectation value of the energy of the system. For 
simplicity ¢ and 7 are taken to be of identical form, the most natural 
choice being 


$(k)=(k)=G(k)/(Kyw) . . . 2 ee + 0) 
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where K2=/G?(k)d%k/w2. A perturbation solution of (3) suggests the 
suitability of the form (5) for ¢ and %, and in fact for neutral mesons such 
a product form is an exact solution with (5) and suitable coefficients c,,. 
It is clear that this is an adequate approximation in the limit of ‘‘ weak ” 
coupling, where the possibility of having more than one meson present. 
is negligible. In the case of strong coupling Tomonaga shows that this 
wave-function leads to the same energy value as that given by the more 
rigorous theory of Wentzel. Physically this is clear, since the correlations. 
between positive and negative mesons are unimportant when only y,,, » 
involving many mesons contribute appreciably, and further, a wave- 
function of the form (4) is an exact solution for neutral mesons. Tomonaga 
also presents numerical calculations covering the region of intermediate 
coupling between these well-known limiting cases. It is clearly of interest 
to investigate to what extent the neglect of the correlation between the 
mesons, implied in (4), is justifiable in this intermediate region. _ 


§ 3. 
As in Svartholm’s (1945) treatment of the momentum-space integral 
equation for the deuteron, we use the fact that the Schrédinger equation 
.(3), which has the form 


(E—T)¥=—lVY, 3s... 2S 2 eee 
may be derived from the variational principle 6J=0, where 
J=(Y*(E—T) PCP AVG oe eee 


The minimum value J,,;, is simply the coupling parameter / for the 
assumed energy value. 

In the present work it has been found convenient to restrict the wave- 
function (4) by requiring c,, ,, to be of the form 


Cn, manny nt mI NE 2. | ae 


0 2n R2n 2n+2 Q-n 
es Ga a 5 an )=h228)+ phi(a8) 


where 


n=0 
I, and I,, being the Bessel functions of imaginary argument. With the 
wave-function (4) (referred to as y‘) then at once 


a 


(E—T)qp= (y'*(E—T)y) =) {224Iy(248)+(208— 4) 1,298) , (9) 


Voo= (x'O#V x) =" 2K ho: {Ty(2e8)-+1,(2e8)}, . . . . . . (96): 


where A= Jw¢?(k)d?k=2°671 wp. The ratio (E—T) 9/Vo9 is then minimized 
by choice of «, 8. This may be carried out analytically in the limiting 
cases of weak and strong coupling. With the notation «2=—B/), for 
weak coupling a=«+..., B=de+... and Kyl=«(1+x2/4+...), and 
for strong coupling «=B=«x with K l= /2«(1—1/(4«2)+...). Clearly 
the value of / does not depend critically on the choice of « and f, since 
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Z is stationary. The values given here for / agree thus far with the 
corresponding expressions obtained from the more rigorous theories in 
these limits. For the several values of E chosen in the intermediate 
region it was found that the energy value obtained using the approxi- 
mation (8) differed by less than 0-2 per cent from that obtained from the 
best possible choice for 


nm, m* 


§ 4, 
Svartholm (1945) has proposed an iterative method for the solution of 
eigenvalue problems of the type (6), (7). From an initial function Y 
a sequence of functions {¥”} is defined by the relation 


ADA OB AC ehh Attar) eo, vere re Okt) 
Under the conditions that (T—E) and V are positive definite, he has 
proved that the sequence Jp, lo, 14, .-- 5 Ins Unater Ings --- i8 Monotonic 


decreasing, where P(E) ial Vion and b,41)2= V nn| (T—E)y 41, n+1° 
In the weak coupling limit the first steps of this process may be carried 
out directly. At once 


1 0 
Xoo =K,xx Nas ? 


XW=K yy(1+-oBrfo(I—Aefo+...), $ 2 2. (LD 
X1P=Kox\Paws/(w1 +o). | 
whence 
(He), =2AK2(1 +. Fx2/2-+...); | es 
Vy =2AK3(1+ (4 —43)e?7-+. -+)s 
where ¥ is the definite integral [[w3¢?(k,)¢2(k,)d°k,d°k./(A(wy+a9)). 
(In fact ~—0°568.) The sequence {I,,} is then 
Ip=r(1+n2/4-+...)/Ko, \ any 
b= —K(L (1 — 9 3/21. . -)/ Ka. J 
It is apparent from (13) that the important modification of the wave 
function is that y‘” allows to some extent for the correlation between 
the various mesons. In the same way, (/)—I,)//) has been calculated as 
far as terms of order x*, as shown in curve (A) of the figure. - 


§ 5. 
In the general case, the iterated functions are given by 


XO (k; K)=Kox0,(k; K){(2xo)/B+Aaj}/(E-2,o—Ano), . (14a) 


XO (kK 3 K)=Kox) 5 n(; K){(2,@)/a+A8 j/(E—2,@—2qw), (14 6) 
from which (E—T),, and V,, may be calculated. For this the integral 


representation 


ee (ei Pa ee wk (LS) 
—EK+2w 0 
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is found convenient, and the expression finally obtained is 


(E—T),;=—Kop iF eb dt | (t+ 31) (D,+ 2AaBDj/) 
a 


+208 (Stott) /ut+Diep (1.45%) |. Pieri 


in which’ I,—I,(2«8D,), I,=1,(2«8Dy), I,=1,(2«8Dy), Do, D,, and D, 
being the integrals D,,(t)=Jw%e~°/"¢?(k)d?k/u”, expressible in terms of 
Bessel functions K,, and their integrals. A similar but more complex 
expression may also be obtained for V,, but this will not be required here. 

The expression (16) has been obtained numerically for three energy 
values, and the coupling constants J), 1, calculated. The results are 
expressed in the table. 


—(E—T)o9/u| —Voo/u —(E—T),;/p lo 


13-89 16-03 18-88 0-866 
6-64 11-56 20-65 0-574 
3°44 8-86 23-86 0-380 


Table of values of «, 8, (E—T)o9, Voo, (E—T),, U9, and J, for some values 
of E. 
Fig. 1. 


oox €4,-L,)/2, 


—"E Se 


The proportional change produced in J, by the first iteration, as a function 
of the energy. 


In the limit of strong coupling, asymptotic expressions for Iy, I,, I, 
may be used to obtain an analytical expression for (E—T),,. With 
a=f=« this leads finally to 

(E—T)47= 


es 


Tan MI (24"D (0)—A2)/(8Ax*) +...) 


whence 1,= y/2«(1—p?D,(0)/(4A2«2)-+...), in comparison with the value 
b= 2«(1—1/(4«°)). The proportional improvement in the value of J, is 
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sketched in the accompanying graph, together with the curves (A) and (B) 
obtained in the limiting cases. 

It should be noted that because of the arbitrariness in the cut-off used 
in (1) the physically significant quantities are not 1), l,,... but « and B, 
which by (8) give a measure of the average number of mesons present in 
the system. 


§ 6. 
The eigenvalue equation (6) will have a series of eigenvalues 1, 
1%,...U™,.. . with corresponding eigen-functions py, . . . ,,, . . . Satisfying 


the orthogonality relations (,,, Vyz,,)=0 unless m=n. If the initial wave 
function Y° for the iteration procedure is 


(z= 
i =2Anihy; 
n 


then the rth iterated wave function Y” ig 
De 2s (LM), 
via 


so that if /{,... greatly exceed the lowest eigenvalue 1‘, the rate of 
convergence to the wave function % is rapid. The sequence of 
approximations {/}} of Svartholm may then be expressed in the form 


2'| On PF nn|(Um) 

n 

Np = py SE ATER 

: 2| An ? Pre MUON ate 
n 


where 7,,,,, is the expectation value of (T—E) in state ;,. Accordingly the 
first step in the sequence {/{}} should provide the major correction to [) 
provided the ratios 1//™ are sufficiently small, which will be the case 
if there do not exist any low-lying excited nucleon states—in fact, with 
this assumption, 

(10 — 1) (HIM) —TO/UD 
and 


@ 2 
(9—41j2)[loe * | F nnlT oo» 

n=117O 
if Y is a good approximation to sy. In the present work we have found. 
that (J)—1;)/l) does not exceed 3 per cent so that the wave function of 
ines appears to be a surprisingly good approximation to , even 
in the region of intermediate coupling. 
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ABSTRACT. 


Polytypism of silicon carbide can be understood if the assumption is 
made that the crystal grows out of the primary nucleus as a single sheet 
several atomic layers in thickness. The sheet winds itself spirally 
upwards along a bunch of screw dislocations which are situated at the 
spiral axis. If there are stacking faults in the original growth layer, 
they repeat: themselves with the periodicity of the pitch of the screw 
throughout the crystal and thus determine the type and the dimensions 
of the unit cell. Similar typism occurs in graphite and possibly in 
diamond and other substances. 


Arter Frank (1951) drew our attention to the existence of spiral 
growth-steps on carborundum crystals, we examined sevéral hundred 
well-developed crystals of artificial silicon carbide. The crystal surfaces 
were not treated in any way and were examined under an optical 
microscope, using oblique illumination at low magnification. Spiral 
growth markings, similar to those described by Verma (1951) in the 
preceding paper, were observed on about 3 per cent of the crystals 
examined. Some of the steps are obviously many atomic layers in 
height, whereas some are revealed by preferential deposition or crystal 
overgrowth of foreign matter, the step acting as a nucleus for overgrowth. 

In some cases, preferential etching seemed to take place along the groove 
of the spiral. In one case a previously invisible spiral has been made 
plainly visible by smearing and wiping off grease from the crystal face. 
Some grease deposit remained caught near the spiral step, thus modifying 
the reflecting power of the crystal surface. 

Curved spirals of a type described by Menzies and Sloat (1929) and 
straight-edged triangular and hexagonal spirals were observed. The 
outline of the crystal usually bears a close relationship to the shape of 
the spiral. 

We have reached independently similar conclusions regarding the 
origin of the polytypism of carborundum to those put forward by 
Frank (1951) in the preceding paper. 


* L.C.I. Research Fellow. 
+ Communicated by the Author. 
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In carborundum crystals the [0001] face, which shows the spiral 
growth-steps, is the slowest-growing face. If we limit ourselves to the 
conditions of growth at low supersaturation, then the presence of a 
primary nucleus initiating a single screw dislocation line would give rise 
to a needle with the dislocation line as its axis. The cross-section of the 
needle would remain constant and equal to that of the primary nucleus. 
Similarly, two straight screw dislocation lines which are not parallel 
would give rise to a very thin plate of a thickness equal to that of the 
primary nucleus. 

Assuming straight dislocation lines, a crystal can grow in three 
dimensions only when three screw dislocations are present, not all in the 
same plane. This rule does not hold for curved disloeation lines, but in 
general, three-dimensional growth cannot proceed unless there are points 
of emergence of dislocation lines on at least three different (non-parallel) 
faces of the crystal. 

One can imagine that the first stages of the crystal growth are very 
thin needles or plates based on one or two straight dislocation lines. 
However, when dimensions of these exceed certain critical limits, they 
become very fragile. Depending on the stresses imposed by the environ- 
ment, deformation and breakage would ultimately occur. If the broken 
parts carry with them screw dislocations, as needles must do, they can 
serve as new nuclei (this probably accounts for the hitherto unexplained 
multiplication of ice nuclei in the high atmosphere), whereas cracks and 
slip deformations in thin plates, as suggested by Frank in the preceding 
paper, would give rise to growth-steps having large Burgers (1939) 
vectors. 

It is likely that in carborundum crystals, small steps grow faster than 
larger ones. The original short screw dislocations, which initiated the 
plate growth, would in most cases remain unaltered by cracking, so that 
the plate sides would remain fast-growing. However, if the large face 
of the plate acquires a screw dislocation with a large Burgers vector, it 
begins to grow, but remains the slowest-growing face. Even if further 
dislocations with the small Burgers vector appear on the same face, the 
large dislocation would tend to remain dominant ; the fast-growing small 
steps would merge with the large steps when overtaking them. 

The crystal can be essentially pictured as a single growth-sheet several 
atomic layers in thickness, which winds itself upwards spirally using the 
large Burgers vector as its axis. Provided there is a mistake in the 
growth-sheet and the dislocation is stable, 7.e. no new dislocations 
appear or disappear during growth, the crystal would exhibit periodicity 
determined by the length of the Burgers vector and the type of the 
mistake in stacking. 

It is most likely that silicon carbide is only one typical example of the 
occurrence of polytypism by spiral growth from screw dislocations 
possessing a large Burgers vector. Similar behaviour is to be expected 
in numerous other substances, the crystals of which grow out by ‘the 
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screw dislocation mechanism. Silicon carbide shows this type of 
behaviour to a marked extent because of the great stability of its screw 
dislocation system. Other softer substances probably do not possess 
such stability and the phenomenon would occur only on a microscopic 
scale. Graphite behaves in such a manner, as polytypism similar to that 
of carborundum has been observed by Lipson and Stokes (1942 a, 1942 b). 
by X-rays and Hoerni (1949) on individual microcrystals by electron 
diffraction technique. However, the most interesting case might prove 
to be that of diamond, which should probably behave in the same manner 
as silicon carbide. Curved faces and edges appear on diamond. Tolansky 
and Wilcock (1946) have shown that the curved faces are the result of | 
successive stepped growth-sheets. These almost certainly originate as 
spirals at growth-points located near the centre of the face. It is well 
known that at least two modifications of diamond, types I. and I1., 
exist. It is therefore to be expected that if not the whole, then at least 
substantial parts of each crystal of diamond should be all of the same 
type, as in silicon carbide crystals. Thibault (1944) also mentions the 
minerals coquimbite and paracoquimbite, both Fe,O, . 380, . 9H,O and 
parisite and synchisite, both (CeF), Ca(CO 3), as displaying polytypism 
similar to silicon carbide. 

I wish to thank Professor J. Monteath Robertson, F.R.S., for his 
interest, Dr. I. M. Dawson and Dr. F. C. Frank for valuable discussions 
on the subject of dislocations and Professor M. J. Buerger for pointing 
out the necessity of three dislocations and Dr. E. D. Currie and Dr. A. 
Taylor for the loan of carborundum samples. 
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SUMMARY. 


A technique is given for carrying out rapidly the computations involved 
in fitting a polynomial by least squares to unequally-spaced data. The 
technique is obtained by adapting the “ triangular resolution ’’ methods. 
of solution of linear equations and inversion of matrices given by Fox 
(1950) and Fox and Hayes (1951) in a way similar to that in which Guest 
(1950) has adapted a slower method of Doolittle. The method proposed 
affords the additional advantage of being able to decide, early in the 
computation, which degree of polynomial curve is appropriate to the 
data. The theory is presented to establish the validity of the technique 
and two numerical examples are carried out to deal with two different 
sets of requirements. 

References are made throughout the text to the appropriate sections 
of the above three papers, but the present paper is intended to be complete 
in itself, and does not require a prior knowledge of those papers. 


§1. INTRODUCTION. 


THE use of orthogonal polynomials in the fitting of regression polynomials 
to equally-spaced observations is now well-known, accounts of it having 
been given by several authors, e.g. Fisher and Yates (1943). In this 
case the method has considerable advantages over the straightforward 
method of solving the normal equations derived from the power terms 
of the regression polynomial, but these advantages are lost in the case 
of unequally-spaced observations and it is usual to fall back on the 
straightforward method as the quickest way to carry out the computation. 

Guest (1950) has, however, given a method of computation for this 
case which, while solving the normal equations of the power terms, 
adapts the solution to keep in view the technique of orthogonal polynomials 
and thereby preserves many of the advantages of this technique. No 
work is saved in the calculation of the regression coefficients, but Guest’s 
method does obtain all curves of degree less than p, and their associated. 
residual sums of squares, in the calculation for the curve of degree p, thus 
enabling the curve of degree best suited to the data to be selected at the 
end of the calculation. This process has the further important advantage 
that the calculation of the standard errors of the coefficients of the power 
terms and of the standard errors of points on the estimated regression 


~ curve is considerably simplified. 


* Communication from the National Physical Laboratory. 
5A2 
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The method of solving the normal equations which Guest adapts for 
this purpose is that of Doolittle, but much of the gain due to using such 
a systematic method of solution is nullified since he writes down the 
results of all the multiplications in the computation before summing them 
instead of, as is usual, accumulating the products on a desk calculating 
machine and writing down only the sums. Consequently, as Guest states, 
the arithmetic becomes extremely tiresome in the fitting of curves of 
degree higher than the third. This extension would, however, be quite 
straightforward if Doolittle’s method were followed rather more closely. 

There are now available, however, methods of solution of linear equations 
and inversion of matrices which are far quicker than this method of 
Doolittle. Some of these, depending on the technique of ‘ triangular 
resolution ’’, have been recently presented in papers by Fox (1950) and 
Fox and Hayes (1951) with particular emphasis on arrangement on the 
computing sheet. We are concerned here only with those methods for 
symmetric matrices. The present paper adopts these methods, as Guest 
has done with Doolittle’s method, to the fitting of polynomial curves to 
unequally-spaced data, preserving, at the same time, the above-mentioned 
advantages of orthogonal representation. 

The major gain in the proposed method, beyond that of the use of a 
quicker method of inversion, is due to the inclusion of a method of 
obtaining very early in the computation—in fact, in its first stage—the 
contribution to the total sum of squares of each of the orthogonal 
polynomials. With this information, we are able to determine which 
degree of polynomial is best suited to the data and thus carry out the 
rest of the calculation for a polynomial of this degree, thereby saving 
the work of computing, for example, an eighth degree polynomial when 
a fourth degree one proves to be quite adequate. 

Moreover, the elementary matrix theory given by Fox (1950) to develop 
the “triangular resolution” technique can easily be extended to fulfil 
our present purpose, and enables the theory underlying our method of 
fitting polynomials by least squares to be presented very concisely. 

The theory of our method is presented in the following sections, but 
its only purpose is to demonstrate the validity of the computational 
technique given in the examples, which can easily be carried out by 
computors having no knowledge of the matrix theory. 

Two examples are given which involve slightly different forms of the 
basic method, designed to meet two different sets of requirements. The 
data used in the examples are those given by Guest (1950), to enable a 
direct comparison to be made with his method. 


§2. THEORY. 
(a) Fitting of Power Terms. ; 
The usual algebra for the fitting of the power terms of the regression 
polynomial by Jeast squares is now outlined in matrix notation. Matrices 
will be dencted by capital letters in Clarendon type, except for those 
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with only one column, which will be denoted by small letters in Clarendon 
type. 

We have a “ dependent” variable y which we take to be related to 
the “ independent ” variable x by the equation, of degree q in 2, 

OX DalaO, 0, tant ter eit ee fe ide tel’) 

Bete m =| lie 2. 08); b'=[b, by... b,|, the row matrix of the 
coefficients of corresponding powers of x, and 5 is the error term, which, 
as usual, we take to be normally distributed with mean zero and variance 
o” independent of a. 

Our problem is to estimate the b,, the standard errors of these estimates, 
and the standard error of any point on the estimated regression curve. 

We make pairs of observations x,, y, (s=1, ..., ) of x and y 
respectively, and our n equations of the form of (1) are represented by 

Ve DEE Ones i an tay eee te iy 

where y —[y, y..--Y,], X=[x: xX....x,] and 8 is the column matrix 
of the random errors 6,, 55, ..., 5,, which are assumed to be statistically 
independent of each other. 
Regression Coefficients. 

We then minimize 25?= 6'6 = [y’—b’X][_y—X’b] with respect to b, 
obtaining the normal equations 


A 


XX'b=Xy, 
where b is our estimate of b. 
We shall write this as “ 
Ab=h, (3) 
where 
A= XX’, ae ED Bol Se i ona Rea Ce) 
h=Xy. Py ee ae Oe ea) 
Our estimated regression curve is then 
y= x tome aor EAT (6) 


where b is the solution of equation (3). 
Residual Mean Square. 
The sum of squares of residuals is 
[y’—b’X]ly—X’b] 


which reduces to 


y'y—b'h. 
The mean square of residuals, which estimates o? is then 
ee ey AAD Ce 


Standard Errors of Regression Coefficients. 
We have, from equations (2) and (3), 
b—b=A-!X8. 
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Thus, since E(6)=0, : 

E(b)=b 
and the sampling covariance matrix of the b , is therefore 
E(b—b)(b’—b’)=E(A“!X88'X’A-}). 
But, since further the 5, are independent and have variance o?, we have 
K(88’)=o7I, 
where I is the unit matrix of order n, and therefore 
E(b—b)(b’—b’)=02AXX’A-1, 


which, using equation (4), 
== 07 AT) "2 ipo © ee 
where o? is estimated by s%. 
Our estimate of the sampling variance of b, is therefore 


V(6,)=sza" Tyo oe ey acon atk ee eee 
where a”) is the element (7, 7) of A~!, r=0, ..., q. 


Standard Error of Estimated y. 


If now, for any given value of 2, we wish to estimate the corresponding 
value of y, we use our estimated regression curve 


y=x'b 
and the sampling variance of this estimate is 
E{y—E(y)}?=E{x’(b—b)(b’—b’)x}, 
which, by equation (8), 
=o? x Ati x) ) ac... sits Cn 
where s%, is again used to estimate o? 
Our estimate of the sampling variance can therefore be written as 


r=08s=0 
Ss 8 (r; 8 
Le. r=08=0 
giving 
5 
“<q 
2 
V(y)S8h ee Wags yes ss ee 
h m=0 
where 
(Fh eee r, 
Us verdes! Osea UZ, 820. 
r+s=m 


Our computation, therefore, has to provide the solution, b, of the 
equation (3), and also the inverse of the matrix A so that we can determine 


V(b,) and V(y,), our estimates of the sampling variances, from equations (9) 
and (11). 
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(b) Orthogonal Polynomials. 


We shall now give the algebra for the least squares fitting of the 
regression curve by means of orthogonal polynomials in order to 
demonstrate how this approach can be of assistance in the fitting of the 
power terms of the regression polynomial, dealt with in the previous 
section, particularly with regard to an early decision on which degree 
of polynomial is best suited to the data. The methods for the solution 
of linear equations and the inversion of matrices dealt with by Fox (1950) 
and Fox and Hayes (1951) lend themselves easily to this approach. 

The basis of these methods for the symmetric case, in which we are 
interested, is the resolution of the matrix A into the product of a lower 
triangular matrix L and its transpose L’. That is, we obtain L such that 


ASE TRS raat er dot tae 12) 


where L has zeros above the principal diagonal. 
The orthogonal polynomials, p,, which we shall consider in fitting our 
polynomial curve are given by 


TT 
Doan be) ass = Oa llme art 
s=0 


where /‘” 5) is the element (7, s) in L~}, also a lower triangular matrix. 
In matrix notation these can be written as 


ple en ee ee are oes ( 13) 
Our n observations, p,,..., p,, of p are then given b 
Pe=Laik. 


where P=[p, p2..- ppl. 
That these polynomials are, in fact, orthogonal can be seen at once, 


since 
PPA Lt XX (1) 
=i; A(L.!)” from (4), 
Det Lins), from (12), 
giving 
ie DAs es a loa folkd NN oe ess (14) 
where I is the unit matrix of order (q+ 1). 

These orthogonal polynomials are not the ones used by Guest (1950), 
but are multiples of them. One coefficient in each of the (y+ 1) polynomials 
can be chosen arbitrarily : Guest has made the coefficient of the highest 
power of x in each polynomial equal to unity, whereas we have chosen 
the polynomials which produce equation (14). 

In fitting to our observations a regression curve in the form of a sum 
of these orthogonal polynomials, we write the fundamental equation (1) in 
the form 

Sienl Hos hen eG pee a ee UP 


where 
p'c= x’b. 
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Similarly, in place of equation (2), we have 
Y= P48 a ee) ee (2’) 
and we estimate e by é, obtained by minimizing 6’$ with respect to c. 
The algebra is then exactly the same as in the previous section with 
X replaced by P and b replaced by ¢. 
Also, therefore, 
A( = XX’) is replaced by PP’, which is equal to i; 
and 
h( = Xy) is replaced by Py, which is equal to L~th. 


A 


Our estimate ¢ 
equation (3)) 


is therefore given by the equation (corresponding to 


A A 
G=L71hye ot ator ee 
and our estimated regression curve is 
yp Ck Se ee 


This curve is, of course, identical with that of the previous section 
(equation (6)), as can easily be shown algebraically using equations (3’),. 
(13), (12) and (3). 

Similarly, equation (7) becomes the equivalent equation 


1 x 
s2= ——— (y’y—c’L“"h), 


n—q—l1 
or, in a more usual form, from equation (3’), 
1 
2—= ——_— (y’y—¢’e), 
SR n—q—1 (yy ) 
or mi lach we ao ema 


2 1 x 2 x AQ 
aa ae 
Finally, equation (9) becomes, 
Vice)=s85; ok Ss taal ative Gren vein ca 
which is not usually of any interest to us, and equation (11), written in 
the form of equation (10), becomes the equivalent equation 


V(y,)=s%, “p'pe <P a ed 
It is equation (3’), together with equation (7’), which proves to be of 
value in obtaining the standard errors associated with the regression 
curve. We know that ¢? is the reduction in the total sum of squares due 
to fitting the orthogonal polynomial of degree r, and, therefore, what is 
equivalent, ¢; is the extra amount removed from the total sum of squares 
when a power series of degree r is fitted instead of a power series of 
degree (r—1). (It may be noted, incidentally, that this latter fact can 
easily be established as an algebraic identity without reference to 
orthogonal polynomials, ¢ then being defined by equation (3’). The use 
of these polynomials, however, allows a better understanding of the 
computational procedure and of the meaning of the numbers that occur: 
in the calculations.) 


Fitting of Polynomials to Unequally-Spaced Data 1393 


The coefficients ¢, are, in fact, produced as by-products in the first 
stage of the solution of linear equations for the symmetric case, as given 
by Fox (1950) in his Example 4: L is produced a column at a time, and 
as column r is produced so also is the element 7 of a row matrix z’, which 
is the solution of the equation 


Lz=h: 
(In Fox’s notation this is Ly=b.) That is, 
Z=L~th. 


z is, therefore, from equation (3’), our column of coefficients ¢,. By 
_ considering the reduction in the total sum of squares due to each ¢, in 
turn as it is produced at the end of column r of L, we can decide which 
degree of power series is most appropriate to our data, and can complete 
our computation for this degree. We thus save the waste of labour 
which is consequent on calculating a regression series which proves to 
be either too short or too long. 


§3. DerarLs or ComPUTATION. 


A decision has to be taken as to what degree of polynomial curve is 
to be considered at the start, that is, what is the highest power of x we 
are going to use in calculating the sums of squares and products (the 
elements of A and h). Since, however, we shall know very early in the 
subsequent caiculation which degree of polynomial is, in fact, appropriate 
to our data, this decision is of only minor importance: if we find that 
we have used too high a power we merely delete the unnecessary elements 
of A, L, etc. ; if we find that we have used too low a power we calculate 
the further sums of powers and products and elements of L as required, 
incurring no unnecessary calculation other than the trivial adjustment 
of the row and column sums, s, and s,, referred to below. 

To facilitate comparison between the method here presented and that 
given by Guest (1950), we shall not concern ourselves with this question, 
but shall carry out the computation exactly corresponding to that of 
his example, the fitting of a cubic curve to the sixteen pairs of observations 
he gives. These are 


8 
< 


Or 


y 
“4: 
5 
0 
ah 
8 
d 
ck 
“4 


0-40 
0-53 
1-09 
0-76 
- 0-80 
ne Wy 
ESLO 
1-45 
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Also following Guest, we shall not concern ourselves with the calculation 
-of the power sums of the x, or the sums of products of the y, with the 
powers of the x,, but shall assume these to have been obtained by the 
method he suggests. In other words, we shall start our computation 
' having already obtained the matrix A (the element (r, s) of which is given 


by a,,= 2 att’, r, s=0,...; q) and the column h (the element 7 of which 


1=1 
is given by h,= & y;2,"). 
i=! 


We shall, however, work in terms of 2/10 instead of a, since this makes 
all the diagonal elements of A of roughly the same size and order unity, 
permitting the advantage of keeping a fixed number of decimals throughout 
the calculation. 


Method I. 


In outline, this method obtains the regression coefficients, b, of the 
power terms from equation (3) by the technique given by Fox (1950) 
in his Example 4, the first stage of which produces L and ¢. The latter 
-enables the residual mean square to be calculated from equation (7’), and 
the decision to be taken as to the degree of polynomial appropriate to 
the data. The inverse of A is then computed from L by Method II. of 
the paper by Fox and Hayes (1951), and so we are able to obtain the 
standard errors of the regression coefficients by using equation (9) and 
of points on the regression curve by using equation (11). 

_ The computation is set out in Example I. below. 

A and h are placed side by side on the computing sheet. For checking 
purposes, the column sums of A and the sum of the elements of h are 
entered in the next row as s’,, and the row sums of A and h together 
are entered in the next column as s,. 

We then calculate L’ and ¢ (this is a slight variant of Fox’s technique, 
which deals with L and ¢’) a row at a time from the equations 


LUA. eidcn *+ tua oe een a 
and 
Le=h, . . (alternative form of (3’)) 
in the form 
column r of L’ x column s of L’=element (r, s) of A. mo) 


Here s takes the values r up to (¢+1) in turn for each r from 0 up to q in 
turn, where the column (q-+1) of L’ is to be read as € and the element 
(r,q¢+1) should be read as the element r of h, r=0,..., g. At each 
‘Stage, every element in equation (15) is known, except the element (r, s) 
-of L’, or the element 7 of ¢, which can therefore be obtained. 

When row r of L’ and € has been thus produced, the elements in this 


row are summed to give element r of s;. The computation of the row 
is then checked by 


column r of L’x the column s,—element r of s,, . . (16) 
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in which all the elements are known except for those in s,; which are 
multiplied by zeros in L’. 
It is at this stage that we can decide whether or not a polynomial of 
‘degree r is sufficient for fitting our data, and so whether or not we need 
‘compute any more rows of L’ and ce. We enter ¢? in position r of the 
column “sums of squares due to orthogonal polynomials ”’ and subtract 
this from the entry (r—1) in the column “ residual sum of squares ” to 
obtain the entry r in the latter column. (For r=0, we subtract ¢? from 
_ 2y?.) In the next column we enter the appropriate number of degrees 
-of freedom (=n—r—1, for r=0, 1, 2, .. .) and so obtain the residual mean 
square. These last four columns could be entered beside ¢, but to save 
computing columns have been entered as a separate table at the end of 
the main calculation. We can then test the hypothesis ¢, (and coefficients 
of all higher degree polynomials)=0 by an F test on the ratio of the sum 
of squares due to the orthogonal polynomial of degree 7 to the residual 
sum of squares in the same row, since, under that hypothesis, these are 
independently distributed as y?o7, on 1 and (n—r—1) degrees of freedom 
respectively. 
We proceed in this way until we obtain a ¢,, say ¢,, which does not 
reduce the sum of squares significantly. We then delete row k of L’ 
-¢, Aand hand column k of L’ and A, and all subsequent rows and columns 
as appropriate, adjust s’, and s, accordingly, and proceed to the next 
-stage of the calculation using the remainder of L’ and e. 

If, on the other hand, we complete the calculation of L’ and ¢ without 
obtaining a ¢, which does not reduce the sum of squares significantly, 
we must compute further rows and columns of A and h and, from them, 
further rows and columns of L’ and é until we do obtain such a ¢,, and 
then we proceed as before. It should be noted that such an increase 
in the size of A and h, or the final decrease in their size, does not affect 
-the elements of L’ and ¢ already obtained, the only alteration required 
being to the row and column sums s, and s3. 

As stated earlier, however, we shall not concern ourselves in our example 

with this:problem of deciding on the appropriate degree of polynomial 
-to use, but shall carry out the straightforward calculation for fitting the 
-cubic. We are assuming, in effect, that we have carried out our 
computation so far with five rows and columns (q=4) in A etc., found 
that ¢, does not reduce the sum of squares significantly, and consequently 
deleted row 4 and column 4 of A etc. and adjusted s3; and s;, leaving 
ourselves with the elements given in the example. 

The next stage is the computation of b. 

We have, from equations (3) and (12), 


LL’b=h, 
_and, therefore, : 
L'b=L-'h, 
-which, by equation (3’), gives 
be: See ae co (he) 


1396 J. G. Hayes and T. Vickers on the 


To obtain b we use equation (17) in the form 
row r of L’x the row b’=element r of ¢€, 
allowing + to take the values q down to zero in turn, so obtaining 


Oy ae» Ug a SULT, 


b is then checked by the equation 
the row sj x the row b’—total of column h. 


Finally, we compute A~! from L’ by Method II. in the paper by Fox 
and Hayes (1951). From equation (12) we can obtain 


L'A1=L-1, 


and use this, to obtain A~! without needing the unknown part of L~?, in 
the form 


row r of L’xrow s of A~!=element (7, s) of L~1 


==); whenr<s . (18) 
==1/l_,.) when r=s, 
where /,., is the diagonal element r of L’, r=0, ..., q. ‘ 
For this purpose the quantities 1/l,,. are recorded in a row beneath b. 
By letting s take the values q, g—1, ..., 7 in turn for each r from q 


down to zero, we obtain from equation (18) the elements of A~! row by 
row, and in each row from right to left in turn. 

When row 7 has been thus obtained, it is copied into column 7 and 
checked by the equation 


row r of A~! x the row s;=unity. 


Having obtained A~', we can calculate the sampling variances of the b, 
using equation (9), and of points on the regression curve, using equation (11), 
in the form 


V(Ye)=SR(Uo + Uy e+ Uge+. . .+Uy,0?9), 


where the uw, are the sums of the elements of A~! in the appropriate 
backward diagonals. Forexample, wv, is the sum of the elements underlined 
in Av! of example I. and uw, is the sum of the elements marked by an 
asterisk. 


EXAMPLE I. 


A h Ss; 
16-00000 —0-03500 4-78877 —0-17961 47-06000 67-63416 
—0-03500 4:78877 —0-17961 2-28074  19-74100 26:59590- 
4-78877 —0-17961 2-28074 —0-26553 16-36283 22-98720 
—0:17961 2-28074 —0-26553 1-29699 8:90970  12-04229 


$3 20-57416 6-85490 6-62437 3°13259 92-07353 


Fitting of Polynomials to Unequally-Spaced Data 1397 


L’ é Ss 
4:00000 —0-00875 1:19719 —0-04490 —11-76500 16-90854. 
2-18831 —0-07729 1-04206 906816 12-22124 
0-91733 —0-14306 3:24722 4-02149 
043430 1-04289 1-47719 


b’ 1-80352 3-13867 3-91435 2-40131 
Ip 0:25000  —-0-45697 1:09012 —-.2-30256 


A- 
0-17596 0-07893 —O0-38575 —0-19341 
0-07893 1:38489 —0O-34720 —2-49547* 

—0:38575 —0-34720 1:31731* 0-82682 

—0:-19341 —2-49547* 0-82682 530177 


V(6,) 0-0213 0-1673 0-1591 0-6405 


Sy2=233-729 


Residual sum | Degrees of | Residual mean 
of squares freedom square 


Ce—sum Of 
ry | squares due to 
poly. degree r 


0 138-415 95-314 15 

1 82-232 13-082 14 0-9344 
2 10-544 2-538 13 0-1952 
3 1-088 1-450 12 0-1208 


43 aN ae 
V(y-) =0-1208 | 0.1760--0-1579 al +0-6134 (35) —1-0812 (ia) 


5 a \6 


J} 


S) 


x 4 
—3-6736 (35) +1-6536 € 


ig 10Le— 10, V(y-) =0-3803, 
> for 7=5, V(y-) =0-0215. 
Method If. 


It often happens that we are not interested in the sampling variances 
of points on the regression curve, and are required to find only the 
regression coefficients and their sampling variances. In this case (from 
equation (9)) we need only the diagonal coefficients of A~? and it is 
quicker to compute L~1 instead of A~1 and to obtain the diagonal 
coefficients of A-1 from this. Our method of computing L~ is that of 
Example 2 in the paper by Fox (1950). 

The computation for this case is given in our Example II. It proceeds 
as in our method I. as far as the production of L’ and é and the decision 
as to the appropriate degree of polynomial except that L, ¢' and Ss, are 


1398 J. G. Hayes and T. Vickers on the 


recorded, exactly as in Fox’s Example 4, instead of their transposes, as. 
in our Example I. That is, equations (15) and (16) are replaced 
respectively by the equivalent equations 
row r of LX row s of L =element (r, s) of A, 
and 
row r of Lx the row s; =element r of s,. 
L and é@’ are thereby obtained column by column. 
Having obtained L, we compute (L~!)’ a column at a time by means 
of the equation 


LES 
in the form 
row r of Lxrow s of (L—1)’=1, when r=s, 
==(); when 7r+s, 
where we let s take the values 0, Te ..., qin turn for each 7 from 0 to q. 
In practice, s takes only the values 0, 1, ..., 7, because of the zeros. 


in L and (L~?)’, the latter having zeros below the principal diagonal. 
When column r of (L~1)’ has been obtained, its elements are summed to: 
give element r of s;, and the column is checked by méans of the equation 


row 7 of Lx the row s;=unity. 


In passing, it will be remembered from equation (13) that element (7, s) 
of (L~+)’ is the coefficient of x” in the Cree one! polynomial of degree s. 
From equation (17), we have 


b=(L-1)/'é 
which can be used, to obtain the coefficient be in the form 
b,=row rot (L=#)' x the row Gs 2:2. eee es 


When all the b, have been obtained, their computation is checked, as 
in Method I., by 


the row b’ x the row s;=sum of elements of h. 


To obtain the diagonal coefficients, a”, of A-! we have from 
equation (12) 
(U2) U4 Are eee 


which can be used for this purpose in the form 
row r of (L-1)’ x row r of (L-1)’=a™”, 
The a” yield the sampling variances of the b, by using equation (9). 


EXampete II. 


A h S; 
16-00000 —0-03500 4:78877 —0-17961 47-06000 67-63416 
—0-03500 4-78877 —O-17961 228074. 19:74100 26:59590 
4:78877 —0-17961 2:28074 —0-26553 16-36283 22-98720 
—0-17961 2:28074 —0-26553 1-29699 8-90970 12-04229 


S; 20-57416 6-85490 6-62437 3°138259 92-:07353 


Litting of Polynomials to Unequally-Spaced Data 1399» 


L 
4-00000 
—0-00875 2-18831 
1:19719 —0-07729 0-91733 
—0-04490 104206 —0-14306 0:43430 


: 11-76500 9-06816 3-24722 1 04289 


q@ 


S, 16-:90854  12-22124 4-02149 1-47719 


(L™)’ 
0-25000 0:00100 —0-32619 —0-08400 
0-45697 0-03850 —1-08377 
1-09012 0-35909 
2-30256 


si 0-25000 0-45797  0-80243 ~—«:1-49388 
‘Db’ 1-80350  3-13864 —-3-91435 ~—-2-40132 
V(b,) 00-0213 0-1673 0-1591 0-6405 


' Table of Sums of Squares, etc. as in Example I. 


As has been pointed out, Method IT. is suggested for the case where 
the sampling variances, V(y,), of points on the regression curve are not 
required. If, however, after carrying out Method II. we desire to 
calculate V(y,), the inverse of A can be obtained row by row for this 
purpose from (L~*)’ by means of equation (20) in the form 


row r of (L~1)’xrow s of (L~+)’=element (r, s) of A~}, 


the check being the same as in Method I. 

Method I. is, however, the quicker method if it is intended from the: 
start to obtain A7?. 

It may be seen further that it is possible, in using Method IT., to obtain 
the regression coefficients for a fitted power series of every degree up to q 
in the process of obtaining those for degree g. This merely requires the 
writing down of all the sub-totals of the summations involved in 


calculating the b,. from equation (19) since the coefficient, b., say, of x” in 
the regression power series of degree s is given by the equation 


aA s 
ba eee, oh ate Crepe Ren ea 
t=0 
where J” is the element (¢, 7) of (L~4)’. 
Thus, for example, the coefficient of x (that is, when r=1) in the fitted 
series of degree | is 0:-45697 x 9:06816=4-14388 ; the coefficient of x in 
the series of degree 2 is this quantity added to 0- ‘03850 x 3:24722, giving 


4:26890; and that in the series of degree 3, which is bs in Examples I. 
and II., is this latter quantity added to —1-:08377 x 1-04289, giving: 


3°13864. 
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All the coefficients obtained from equation (21) are given in the appended 
table. 


TABLE I. 
Mecec: Coefficient of 
IIe es l x/10 22/100 | 23/1000 
0 294125 
1 295032 4-14388 
2 1-89111 4:26890 3:°53986 
3 pa 1-80350 3-13864 3:91435 2-40132 
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CXL. Dislocations in Thin Plates. 


By J. D. EsHetsy and A. N. Srrou, 
H. H. Wills Physical Laboratory, University of Bristol*, 


[Received August 15, 1951.] 


ABSTRACT. 


The stress due to a screw dislocation passing normally through an 
infinite plate or a disc is largely confined to the neighbourhood of the 
dislocation line, in contrast to the case of a dislocation in an infinite 
medium. Two screw dislocations in a plate attract or repel one another 
with a short-range force in place of the inverse first power law for infinite 
parallel dislocations. The stress due to an edge dislocation is not 
essentially different in the plate and infinite body so long as the plate 
remains flat, but in some circumstances the stress may be largely relieved 
by buckling of the plate. 


§1. INTRODUCTION. 


Forty (1951) has studied the growth of tabular crystals containing one 
or more dislocations. Dawson and Vand (1950, 1951) have published 
photographs of thin crystals with growth spirals terminating on 
dislocations. It seems worth while, therefore, to contrast the properties 
of a dislocation in a thin plate and in an infinite body. 

We discuss in detail the following configurations of a straight screw 
dislocation: (i.) meeting normally the surface of a semi-infinite body, 
(ii.) running normally through an infinite plate and (iii.) along the axis of 
a disc, with or without a stress-free hole excluding the core of the 
dislocation. Case (i.) may be of interest in connection with the detailed 
topography of the surface near the point of emergence of a dislocation 
(Frank 1951). Fig. 25 of Forty’s (1951) paper corresponds almost 
exactly to case (iii.). . 

The corresponding problem for an edge dislocation is only briefly 
touched on. 

It is interesting to note that the steps in Forty’s growth spirals may 
be an arbitrary number of lattice spacings (of the order of 100), so that we 
have to deal with effectively the classical Volterra dislocation with. 
arbitrary Burgers vector, in contrast with the theory of metals where 
the Burgers vector is assumed to be limited to one of a few simple lattice 


displacements. 
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§2. A Screw DistocaTion NoRMAL TO THE SURFACE 
OF A SEMI-INFINITE Bopy. 


For a screw dislocation in an infinite body along the z-axis of cylindrical 
coordinates (7, 0, z) the non-vanishing components of displacement and 
stress are 

b b1 
Weg Ue ite owas (1) 


The elastic image field which annuls the traction due to the dislocation 
on the plane z=0 is 


Oras pb is pb 1 

lg — 5 Rags 1 ano eR Ree) es (+ ic) » (2) 
with R?=7?+2? and the other components zero. It is produced by a 
distribution of couples along the negative z-axis, twisting about the 
z-axis and with density proportional to distance from the origin. This 
may be verified by integrating twice with respect to z the expression 
given by Love (1927, p. 187) for a point-couple. The sum of the states (1) 
and (2) thus gives the elastic field about a screw dislocation perpendicular 
to the free surface of the semi-infinite solid z>0. 


$3. A Screw DIsLocATION IN A PLATE oR Disc. 


To solve the problem of a screw dislocation in an infinite plate we must 
annul the traction on the planes z=-+d. It is easy to see that this can 
be done by introducing an infinite series of images of the type (2). A 
useful expression is, however, more easily found as follows. We take 
from the image representation the fact that only w,, 7,, and 7,, do not 
vanish. Then 

T 92=POU9/ OZ, T -p=H(OU9/Or—mU4/7) 
and the equilibrium condition is 
Oty , 1 Og We , OR, 
Ort "ror rt! Ot 
A simple solution is ws=exp (+-kz)J,(kr). Multiplication by a function 
of k and integration yields a more general solution. The solution of our 
problem is easily found to be 


=0. 


b \" sinh kz 


Uomare ir eae (kr) = vijay ll Hes Sea 
for then 
(Teden a= ef" J, (kr) See 
Since 
sech 7= 22(—1)" ” exp {—(2n+1)a} 
and 
eke, (kr) A 


Jo bk eter) 
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(3) becomes, with d,=(2n+1) d, 


: 
ae ta, cecal 1A a PtP] d,te+4/[G, +2) a}: 
exhibiting wu, as the sum of a set of i images of the type (2). Returning 
to (3), uw, can be expressed as a series of modified Bessel functions Ia aby. 
contour integration or by following the analysis of Riemann (1855) in a 
similar problem. In this way the elastic field of a screw dislocation 
passing perpendicularly through an infinite plate of thickness 2d is found 


to be 
(a 0, CT cee g = + pe (sin inn) = K, (sr) sin, 
T= Es (sin $n77) K ( a) cos 3 (4) 
rate Bet ein Inn) K (Fr) sin, | 
i 0, Try—=T 96=T eT z=. 


with summation over integral n. 

The case of a screw dislocation in the annular disc bounded by the 
surfaces z= +d, r=r,, r, (r;<r,) can be treated by replacing each K, in 
the uw, of (4) by a linear combination of K, and I, and adjusting the 
constants to annul the traction on the cylindrical surfaces. The result is 


b Z nT? r isGe 
wee wm Borla) on (ale 
7 ( nar Wat NTZ 
Fee og an Ao ( yi] “) +B,K, (Sr) he od’ (5 
ues z ue nar nar\| _. naz 
Pee eai nies | 
te, Ton —— “00 — Neer or =0, 
where 
td OS Ce 
An= 3 7359 2" T ar, [2d)K(nar,/2d)—1,(nar,/2d)K (nmr, /2d) 


and B,, is obtained from A, by writing I, for K, in the Suitla4 
The energy required to form the Seaton in the annulus is 
od 


Waib[ del dere, (6) 
=—d Ti 
When the outer radius 7, is infinite 
wildy Kalra 2d) 
ba n odd n®(nzr,/2d)*K(n77r;/2d) ‘ 
Ese co anect (on a oe 
nr 3 ( : ERD 
_ pb d 3 F 
ae (soar): rise (7) 
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When 7 ;=0, 7,= © the dislocation at the origin exerts a force 


F=0' [ 1, 0e= es UKs (sr) “2 es 
vid us n odd 
on another screw dislocation with Burgers vector b’ distant r from it. 
Since when a>1 I,,(x) and K,,(a) respectively increase and decrease 
exponentially with x the elastic state is given closely by 
b bz pb z 


U.= — 6, Ug=— ey 


Bar pre Uc 
if the point 7 is a few multiples of d away from both the inner and outer 
edges of the annulus. (9) represents exactly the combined effect of (1) 
together with (2) and its reflection in the plane z=0, with couples of 
opposite hand. In (4) the terms in K annul the effect of that part of 


the couple distribution which would lie inside the plate. 


§4. Screw DisLtocaTion : Discussion. 


It will be seen that the elastic state of the plate (equations (4)) is quite 
different from that of a slab of width 2d marked out perpendicular to 
a screw dislocation in an infinite body (equations (1)). The shear stress 
7», is confined to the neighbourhood of the dislocation. That this must 
be so is clear: rt». vanishes at the surface of the plate, and far from the 
dislocation this state of affairs must persist throughout the thickness of 
the plate. The energy of a screw dislocation in an infinite plate, with 
its core excluded by a stress-free hole, is finite, while for a dislocation in 
an infinite cylinder the integral (6), (wb?/47) . 2d. In (r,/7,;), diverges as 
7, > «©. Again, in an infinite body two parallel screw dislocations attract 
like electrostatic line-charges, with a force (ubb'/27) . 2d/r per length 2d, 
whilst each term in (8) is the force due to a line-charge made up of 
particles attracting with the Yukawa potential const. 7-1 exp {—nar/2d}. 
Because 7,, does not vanish the screw dislocation will also interact with 
an edge dislocation running parallel to the z-axis. Since 7,, is an odd 
function of z there will be no net force between them, but only a couple. 
It is also clear that the image-force attracting the dislocation towards - 
the edge of a lamina will be very small unless the dislocation is only a 
few multiples of d from the edge. In the simple case where the dislocation 
is distant D from the edge of a semi-infinite plate this force is given by (8) 
with b’=b, r=2D. For comparison the image force on a length 2d of 
a screw dislocation running at a distance D parallel to the free surface of 
a semi-infinite solid is ~b%d/D. 


§5. Ep@r DIsLocaTIons. 


The stresses 7,,, To, 7,9 Produced by an edge dislocation along 
the z-axis in an infinite body are derived from the Airy function 
x=const. r In r sin 6, while 7,,=v(7,,++749) and the remaining components 
are zero. To solve the problem of an edge dislocation traversing a plate 
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we must find an image stress system giving 7,,=0 and T,,—= const. sin 6/r 
on the planes z=-+d. The necessary analysis would be similar to that of 
Sneddon (1946). We should have to evaluate integrals like those in his 
§ 5, but with Bessel functions of higher order on account of the angular 
dependence. The results would not be simple, and the case of an edge 
dislocation in an annular disc would be quite intractable. The general 
nature of the result can, however, be made out quite simply. If in the 
solution (Love, op. cit. p. 225) for an infinite hollow cylinder with an edge 
dislocation we replace A by 2Au/(A+ 2) the resulting expressions give the 
stresses and displacements averaged across the thickness of a disc cut from 
the cylinder and having surfaces free of stress (generalized plane stress : 
Love, op. cit., p. 207). The average stresses will fall off as 1/r and the 
elastic state will not be very different from what it was when the disc 
formed part of the infinite cylinder. 

As long as the plate containing the edge dislocation remains flat there 
is thus no widespread relaxation of stress, as there is for the screw 
dislocation. . However, the energy may be reduced by buckling of the 
plate. Take a sheet of paper with a hole in it and make a tuck of constant 
width 6 (small compared with the size of the hole) running from the hole 
to the edge of the sheet. The paper will form a surface given roughly 
by the equation z=(by/2 sin 6)/27. If the paper is flattened out to form 
a plane dislocated lamina a large amount of strain-energy will be 
introduced. On the other hand we should not expect buckling to occur 
in a thick dislocated disc. 

To find out the relations between the Burgers vector and the thickness 
and inner and outer radii of the annular disc for which buckling is 
energetically favourable we should have to use the theory of plates with 
strain in the middle surface (v. Karman (1910)), modified to take account 
of initial stress. These equations are non-linear and it is difficult to 
solve them with the necessary boundary conditions. It seems clear, 
however, that in certain circumstances an edge dislocation will be able 
to relieve most of its stress by slight buckling of the plate, except within 
a few multiples of d from its centre, leading to a state of affairs similar 
to that discussed for the screw dislocation. 
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CXLI. The Influence of Deviations from the Debye Spectrum on the 
Electrical Conductivity of Metals. 


By F. H. J. Cornisu and D. K. C. MacDonatxp, 
Clarendon Laboratory, Oxford*. 


[Received August 3, 1951.] 


In the modern theory of the electrical resistance of metals, as originally 
developed by Bloch (1928, 1930) and Wilson (1936), (see also Mott and 
Jones (1936)), the assumptions made about the thermal vibrations of the 
crystalline lattice are essentially those which lead to the Debye spectrum 
familiar in the theory of specific heats. Bloch himself derived rigorous 
formule valid for low (T<@) and high (T>@) temperatures, and 
Griineisen (1933) suggested a general interpolation formula for all 
temperatures : 
Ts oT 2° dz 
1) med 66 he (e—1)\(1—e-*) ; . ° . e . . (1) 


based on Bloch’s work, which agrees with the rigorous expressions at 
both limits of T. More recently Sondheimer (1950) and Rhodes (1950) 
have extended Bloch’s analysis to deal rigorously with the intermediate 
temperature region and it appears that the Griineisen formula differs 
nowhere by more than ~10 per cent. Since in our work we have been 
concerned with trying to interpret experimental deviations from the 
simple Bloch—Griineisen theory of several hundreds per cent at low 
temperatures, we shall accept (1) as an adequate representation of the 
resistance of an ideal metal with a pure Debye lattice spectrum. 

The work of Born and Blackman (1935) on specific heats has shown 
that the Debye spectrum can only be regarded in general as a rather 
rough approximation to the actual vibration spectrum of specific 
crystals. Thus several detailed calculations have been made, based on 
crystal dynamics, of the actual frequency spectrum (e.g. Helen Smith 
(1948), Kellerman (1940)), and these reveal that the frequency spectrum 
has several peaks. This has the consequence that experimental specific 
heat data tend in general at low temperatures to give an apparent fall 
in the Debye temperature. 

It is therefore somewhat surprising that, as far as we are aware, no 
attempt has yet been made to extend the theory of metallic conduction 
to allow for a vibrational spectrum of more arbitrary form. This is 
probably due to the fact that a direct generalization is not immediately 
obvious as in the case of internal energy and therefore specific heats, and 
also possibly because it might be expected that the electrical resistance 


* Communicated by the Authors. 


On the Influence of Deviations from the Debye Spectrum 1407 


would depend less critically on the particular vibration spectrum. Thus, 
for example perhaps, one might look for a crude dependence on lattice 
entropy (¢. g. cf. Simon 1924, Mott 1934). 

From the fundamental theory of metallic resistance (e.g. Wilson, 
loc. cit.) it is quite straightforward to derive the following expression for 
the resistance as the corresponding generalization of (1) above : 


ee eee (2) 
PT J 9 (exp [{kT]—l)(1—exp[—mjkT)? °C 


where q is the lattice wave number, v is the frequency and y=ugq/2z, 
u being the wave velocity. ¢ is the upper limit of wave number 
determined by the total number of normal modes of the lattice. The 
Debye temperature is then defined in terms of g by KO=hdd. 


1. If then w is set =z (constant for all q), corresponding to the Debye 
spectrum, then (2) above reduces immediately to (1) as expected. 


2. If vv» (constant for all q) corresponding to an “ Einstein Model ”’ 
then we have 


el 1 
PT ~* (exp [Ay /KT]—)(1—exp [—/, kT) 


or, setting hyy=kO,, 


1 Ul 


P~ * (exp[OuT]—D—exp [Onl ‘ 


reducing to p~T/@?, for T>@, and p~1/T exp (—@,/T) for T<@,. 
Since, as Born has pointed out, the highly dispersive optical branches of 
the lattice spectrum may be approximated closely by Einstein oscillators 
it is clear from the latter formula that for T appreciably less than 0, 
the contribution to the resistance from these branches is negligible. 

We have turned our attention initially to the case of the monovalent 
alkali metal Lithium which would be expected a priori to conform most 
closely to the ideal free electron model on which the Bloch theory is 
based. However it is evident from the specific heat measurements of 
Simon and Swain (1935) that there is considerable deviation from the 
simple Debye spectrum, and the experimental @ value drops from 
~380° K. at 60° K. to ~330° K. at 15° K. At the same time electrical 
resistance measurements (cf. Meissner and Voigt 1980, MacDonald and 
Mendelssohn 1950) have shown that down to ~70° K. equation (1) fits 
the experimental results quite well with 6=360° K., but that at low 
temperatures there is very strong divergence. 

In order to apply equation (2) to a metal with a spectrum departing to 
some extent from the simple Debye model, we have considered the v—q 
function as made up of linear portions. The three such basic types 
involving a single discontinuity are sketched in the accompanying figure 
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together with the resulting N—v spectra. In our numerical work thus 


far we have limited ourselves to Case (a), i. e. explicitly : 


0<v<v, O<q<%1; 


y=Uy 


= 
ee 
( - 
vy MeV p= 2 SIM 


v) 


N(v) 


(b) 


(¢) 


ih v 
v—q and N—y» curves for lattices deviating from simple Debye model. 


| Note.—The direct complement to model (c) has been here omitted because a 
retrogression of wave-number, q, is involved leading to a negative value 


of N(v) at the discontinuity. | 


This leads to the following expression for the resistance : 


~[z0(8)+6@)-Ge@)) - + 6 
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ead G(m)= £ if > ae elle 
T/ \O) J, (—nd—e%’ 
i.e. G is the function tabulated by Griineisen (1933). Two sets of («, B) 


values have been used and they lead to the results (normalized at 273° K.) 
tabulated below. There is clearly improved agreement with experiment 


Table of Comparative Results. 


a 


ioe aay 374-5 | 273 80-1 20-4 11-4 6-0 


Observed resistance* | 1-443 | 1-000 | 0-1252 | 0-0013 | 0-00009_ | 0-000007 

Griineisen Eq. (1) 1-433 | 1-000 | 0-125 | 0-00042 | 0-000022 | 0-0000009 
(O=360° K. 

New model : 


1-434 | 1-000 | 0-144 | 0-0015, | 0-00009, | 0-000004 
1-435 | 1-000 | 0-129 | 0-00118 | 0-00010, | 0-000004, 


* With residual resistance subtracted, and normalized to unity at 273° K. 


over the Griineisen formula in the temperature range T<70° K. This 
does suggest that a major part of the apparent discrepancy with theory 
may be attributed to the departure of the lattice spectrum, or rather 
the v—gq curve, from the simple Debye model. It is, however, quite 
probable that part is still due to deviation of the electron characteristics 
from those of an ideal electron gas. . 

The abrupt gap in the frequency spectrum in the model we have 
chosen may appear somewhat surprising at first sight. Physically, 
however, one might rather expect to have v—q and N—yv curves somewhat 
as shown by the broken curves in fig. 1 (a). 

The specific heat for the model chosen is given by 


3 
Ci — (5) (Dies) aleve ices ah caaen 18 (0) 
where («, 8) are as before in (5) and D,=D(@,/T), where D is the specific 
heat for the simple Debye model. The apparent © values for a specific 
heat based on this model with «—1-3, B=3 are: at 60° K., 0=354° K., 
and at 15° K., @=286-5° K. This then compares quite favourably with 
the relative variation of @ observed by Simon and Swain, although the 
absolute values are not in agreement. In any case, of course, as Blackman 
(1951) has emphasized again, the effective spectrum will be different for 
the specific heat and the electrical resistivity, since in an ideal metal 
only the longitudinal vibrations produce significant electron scattering. 
We intend to extend this work to other metals and to consider also 
the electronic thermal conductivity and we hope to publish in more 


detail later. 
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CXLIT. Effects of the Recoil on Allowed B-Transitions. 


By O. Kororp-Hansen, 
Institute for Theoretical Physics, University of Copenhagen*. 


[Received July 9, 1951.] 


SUMMARY. 


The recoil corrections to the shape of the Fermi distribution for the 
-spectrum are calculated, and it is concluded that they are of the same 
order of magnitude as other effects, which are usually neglected. Further- 
more it is shown that the observed half-life and shape of the B-spectrum 
of the neutron excludes any coupling in B-decay depending on nuclear 
v./c. terms only. 


§1. INTRODUCTION. 


WHEN the shape of the f-spectrum is calculated one considers the 
f-decaying nucleus as infinitely heavy, and one neglects the kinetic 
energy of the recoil in the conservation of energy since it is small compared 
with the energy of the B-particle and the neutrino. 

In that approximation, where the neutrino rest mass is neglected, the 
maximum value of the kinetic energy of the recoil, TR?*, occurs when the 
f-particle is emitted with its maximum energy. Tj} is given by 


ee Teeeeatenoe Te Catster: bok (1) 


where T}** and Tj“ are given in MeV. A is the mass number of the 
isobars in question. For the f-decay of the neutron one finds TR** = 0-1 per 
cent of T;** and for He* and BY? TR**= 0-05 per cent of T;™*. 

The effect on the shape of the 8-spectrum is somewhat larger as will be 
shown in a following section. Although it is not necessary for this 
calculation to apply so general methods we shall first give an exact 
relativistic treatment of the phase space calculation for a general three 
particle process in the next section, and in the following section we shall 
apply the results of this calculation to the special case of nuclear B-decay. 

The effects from the recoil will partly enter in the phase space calculation 
and partly in the determination of the square of the matrix element for 
the transition. In the special case of the B-decay of the neutron some 
particular effects from the recoil appear in the matrix element and give 
rise to the exclusions of the pseudoscalar coupling case as responsible for 
the B-decay of the neutron. This special case will be treated in the last 


section of this paper. 


* Communicated by Professor C. Moller. 


1412 O. Kofoed-Hansen on the 


§2. Exact RELATIVISTIC PHASE SPACE CALCULATION. 


The probability of a process where a particle at rest disintegrates into 
three individual particles is proportional to the integral 


P=(47)°f| MP dp, dps dps 5(p,+ Po+ Ps) (M—E,—E,—E;) 

(2) 
where M is the matrix element for the process, p,; is the momentum of 
the ith particle, and E, is the corresponding energy including the rest 
mass m;. The abbreviation dp means dp, dp, dp, and the function 6(x) 
is the product of the three delta functions in each of the ocordinates of 
the vector x. M is the mass of the original particle, and c=1. The two 
delta functions express the conservation of momentum and of energy. 

If we assume a random orientation in the space of the original particles 
the process is completely described by two parameters in the momentum 
triangle, and we can average (2) over all parameters besides such two 
variables. 

Three of the integrations in (2) take away the momentum delta function. 
This means, however, that we have to take into account the momentum 
conservation in the form 


Ps =PP+Pe+2p;P2 008A. . . . - + - (3) 
in the delta function for the energy. Here @,, is the angle between the 
vectors p, and p,. By a well-known transformation (see Bethe 1933) we 
get . 

P= 4/E, dE, E, dE, E, dE,| M |? (M—E,—E,—E;). . . (4) 
One further integration leads to | 
P=f P(E,E,) dE, dE, 
=f]. |? (M—E,—E,)E, dE, Hodis.) ee 


P(E,E,) is the distribution function for the two variables E, and E,. 
These two variables determine the process completely. It must be 
remembered that for a fixed value of E, the variable E, cannot assume 
any value between m, and E}**. It follows from the laws of conservation 
of energy and momentum that for a fixed value of E, the variable E, 
must lie in the interval between the two limits 

a =[(M—E,)(M?— 2ME, +m,?+m,?—m;,?) 


Ep; {(M?—2ME, +m,?—m,?—m,")?—4m,?m," }1)/[2(M?— 2ME,+m,?)]. 
(6) 


This expression can be verified e.g. by means of the geometrical picture 
of the conservation laws of energy and momentum given by Blaton (1950). 
Fig. 1 shows as an illustration to (6) the area in an E,E, plane inside 


which these parameters may vary in the special case m,=m,=0-3M and 
M=0-2M. 
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Any one of the particles obtains its maximum energy when the-two 
others travel in the same direction and with the same velocity. This is 


easily seen by a simple extremum consideration. From this result we 
find 


BP*=[M?+m,?—(m +ms3)?]/(2M). a eae Sole ok Gg 


Let us consider a certain process in which |. |? is substantially 
constant. Then we find the probability distribution of one of the energies 
simply by an integration of (5) where |. |? is taken outside the integral. 
The integration must be carried out inside the limits given by (6). We 
find 
Ppy=37P1E, ((M?—2ME, +m,?—m,2—m,2)?— 4m,2m,2! {3(M—E,)? 

- [(M?—2ME, +m,?)?—(m,?—m,”)?]— (E,?—m,”)[(M?—2ME,+m,? 
—m,2—m,2)2—4m,2m,"]}/(M2—2ME,+m,2)> . . . . . (8) 


Fig. 1. 


The hatched area is the area in the E,E, plane inside which EK, and EK, may 
, vary in. the special case m,;=m,—0-3M and m,=0-2M. 


If one wants to determine the total probability for a disintegration 
one has to carry out one further integration over E,, where the integral 
is taken between the limits m, and EP“. This integral is the complete 
relativistic analogue to the Fermi integral in B-decay. The integral can 
be expressed in terms of complete elliptic integrals of the first, second and 
third kind in the general case. Only when all three particles may be 
treated non-relativistically or when one of the particles has zero rest 
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mass can the integral be expressed in terms of simple functions. Fig. 2 
shows the value of the integral of (8) in the special case m,=m.=m3. 
The figure gives this integral divided by 7M*/3840 and as a function of 
m|M where m is the common mass for the three particles. The dotted 


curve is the non-relativistic approximation. The curve illustrates the 
validity of the non-relativistic approximation. 


$3. Recor Errects oN THE PHasE Space CaLcuLaATION IN f-DECAY. 


In the f-deacy the square of the matrix element for the transition is 
given by an expression of the type (Hamilton 1947) 


| M |? =|, |? (1+-a(p,/Hz) cos 9) - . - . - (9) 


Fig. 2. 
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The figure shows the generalized Fermi integral divided by 7M5/3840 as a function 


of m in the special case |. | =1 and m,=m,=m,=m. The dotted curve is 
the non-relativistic approximation. 


where we have neglected the neutrino rest mass. |M, |? is the square 
of the nuclear matrix element, which we shall consider as a constant in 
this section. In the derivation of (9) one has neglected all nuclear v./c. 
terms and all higher order terms in the expansion of the lepton wave- 
functions over the finite extension of the nucleus 7. e. one has neglected 
terms of the order R/A or of higher order where R is the nuclear radius, 
and where A is the wavelength of the leptons. . 
Integrations of the type leading from (5) to (8) are easily performed 
also when matrix elements such as (9) are taken into account. The 
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resulting expression is in many respects similar to (8) but of course it is 
somewhat more involved. If the result is expanded in a series in E,/M 
one finds 


P(E;)=psE,(E3* — E,)?[1 + 3E,/M—Ep**/M + a(m,?/(ME;)—3E,/M)]. 
(10) 


This expression is identical with the usual Fermi 8-spectrum apart from 
the extra factor contained in the square bracket. 

The deviation is small since the extra factor is very close to unity. 
The maximum correction amounts to 0-2 per cent for He®, 0-5 per cent 
for the neutron and 0-6 per cent for B!2 with the choice a=—1. The 
effect is larger for Li’, but due to the complicated structure of the Be’ 
daughter nucleus this B-spectrum is not directly comparable with the 
simple Fermi spectrum. In any case higher order effects may influence 
the shape of the B-spectrum to the same order of magnitude. We shall 
mention that v./c. terms and R/A terms may change the allowed shape to 
the same order of magnitude, namely about 0-1 per cent. 

Thus we may conclude that the effects from the neglect of T, in the 
energy conservation is of the same order of magnitude as other effects 
which are usually neglected in the discussion of allowed B-decay. 


§4. B-DEcAY OF THE NEUTRON. 


Usually the nuclear matrix element -M, cannot be calculated since 
so far no satisfactory theory of nuclei exist. In the case of the 6-decay 
of the neutron, however, one may assume that a satisfactory treatment 
of the behaviour of the neutron and the recoil proton is obtained when 
these particles are described by a Dirac equation. In this case it is a 
simple problem to determine Mx. 

In the scalar coupling case one finds 


| Mx 2—4(EpM-+Mmyp)/(E,M) 
See en oy eS Oe ec Fer) ws (1) 


where the index P refers to the recoil proton. In the pseudoscalar 
coupling case one finds 


IM « ?=3(EpM—Mmp)/(EpM) 
=p". ° 2 ° . e ° ° . . . e e (12) 
We remember that we have put c=1. : 

In case of B-decay of other nuclei than the neutron one would still find 
in the pseudoscalar coupling case that |M, |’ is of the order of vp but 
now vp means the velocity of the nucleon in a bound state, and thus 
| M@|? would be ~104 as large as for the neutron (see also Marshak 
1949) in contradiction to experimental results (Robson 1950, 1951, 
Snell, Pleasonton and McCord 1950). Thus we may conclude that the 
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experimental detection of the B-decay of the neutron has ruled out the 
possibility that the pseudoscalar coupling case can be responsible for the 
B-decay of the neutron. 

Also in case of a pseudoscalar coupling the £-spectrum of the notte 
would not be the same as for normal allowed f-spectra, since the variation 
of (12) with EK, is proportional to 


[pat p, [PoP Se re 


In other words the shape of the f-spectrum of the neutron would be 
of the same type as for first forbidden spectra. This result is also in 
contradiction with experimental results (Snell 1951). 


I wish to thank Professor N. Bohr and Professor C. Meller for their 
interest in this work. 
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SUMMARY. 


' The modifications produced in the Rayleigh-Jeffreys theory of slow 
thermal convection by magneto-hydrodynamic effects in a conducting 
fluid placed in a magnetic field are examined. 

Even for a non-viscous fluid, a critical temperature gradient 8, must be 
exceeded in order that convection occur. In this, the place of viscosity 7 
occurring in the Rayleigh-Jeffreys formula is taken by a quantity yy 
depending on the conductivity of the fluid, the magnetic field strength H 
and depth of fluid d, thus.» = (4/27)(d?/7?)(ou2H2/c2). If the fluid is viscous 
a multiple of the normal viscosity depending on H must be added to ny. 

An estimate of the critical gradient is made for somewhat artificial 
boundary conditions and it is found large enough to be experimentally 
detectable. 

The applicability of Jeffreys’ method of marginal stability is discussed 
and the nature of possible oscillations investigated. 


INTRODUCTION. 


THE interaction between the flow of electrically conducting fluids and 
magnetic fields has been considered by Alfven (1950) and Walen (1944) 
in discussing processes in the sun, by Fermi (1949), Richtmeyer and 
Teller (1949) in a theory of the origin of cosmic rays, and by Elsasser 
~ (1946) and Bullard (1949) in an explanation of terrestrial magnetism. 
Bullard (1949) has proposed a reasonably simple scheme of thermal 
convection currents and magnetic fields within the core of the earth 
which might be able to maintain the observed field at the surface by the 
self-inductive mechanism originally proposed by Larmor. He also showed 
that the magnetic field is of prime importance in determining the flow. 
As a first step toward investigating the hydrodynamical aspects of 
this process a study has been made of the slow thermal convection currents 
set up in a plane layer of conducting fluid placed in a homogeneous 
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magnetic field and heated from below. The motion occurring when no 
field is present has been studied experimentally by Bénard (1901) and 
by Sutton (1951) and examined theoretically by Rayleigh (1917) and 
Jeffreys (1926). 

The stream lines form a system of horizontal vortex rings dividing the 
plane into cells with fluid flowing up in the centre and down at the edges. | 
Rayleigh’s theory described this system, predicted the temperature 
gradient needed to produce it, and made statements about the shape of 
the cells. Here the analysis of Rayleigh and Jeffreys will be modified by 
including the ponderomotive effects of the magnetic field, which involves 
adding to the hydrodynamic and heat flow equations Maxwell’s equations. 
for the electromagnetic fields. 


$1. MarHemMatTicaL FORMULATION OF THE PROBLEM. 


We shall consider an infinite horizontal layer of depth d of an incom- 
pressible fluid with electrical conductivity o placed in a homogeneous | 
maznetic field H, and heated from below. At first viscosity will be 
neglected and Euler’s equation for the velocity V=(u, v, w) written as 


oV : 
of e+ (V- WV] =—vetos tH ixm ©. . (Ly) 


where p is the density of the fluid, p the pressure, g the gravitational 
acceleration, and jx H the ponderomotive force due to the magnetic 
field H, j being the electric current density. 

The motion arises because of the temperature dependence of the 
density, p. This dependence may be taken in a form appropriate to an 
incompressible fluid with a coefficient of thermal expansion «, 


p= poll) ca «oe ee 


The temperature distribution is determined by Fourier’s equation 
modified by the motion to read 


or cere | 
at V.VEISKVAL, eo ee 


where K is the thermal diffusivity. 
In addition, the motion is subject to the equation of continuity 


diy pV—an in. see ee ee 
The current j is given in Gaussian units by 


; V 
Joo(E +p XB) ior ee ee ee 
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where co is the electrical conductivity, and the fields E and H satisfy 
Maxwell’s equations with a dielectric constant e, and permeability ju. 


curl H= <7 j 08 

C C Ob feel 
analy aon LE 

ot 

div H=0, > Piss ee eee me aa.) 
div E47, 2. 

€ 
div g=— 2. 


Here ¢ is the electric charge density and c is the velocity of light in vacuum. 

To render the problem tractable a number of assumptions must be 
-made. First, the equations must be linearized by considering departures 
from equilibrium so small that their squares can be neglected. By 
expanding all functions in powers of a constant A=#/v) where v represents 
the mean velocity and v, is defined by 4pv)2=/(87)H,? one can show that 
the neglect of terms of higher order than the first leads to small errors 
in the solution provided that A is small. That is, the linearization can 
be justified if the kinetic energy per unit volume of fluid is small compared. 
to the magnetic energy density. 

A second assumption is that the fractional variation in the density is 
small and that the continuity equation (1.4) may be replaced by (Hales 
1937) 


divcV= 0: eee ieee A C7) 


The third assumption is that a physically meaningful result, can be 
obtained by solving the time independent problem. Jeffreys (1926) gives 
an argument justifying this procedure derived from the principle of 
marginal stability. ,In a later section the modifications produced in this 
argument by consideration of the magnetic field are discussed and it is 
shown that Jeffreys’ argument remains valid provided c#/47yKo<1. 

With these assumptions the zero’th approximation yields the condition 
of static equilibrium. If a Cartesian coordinate system with the origin 
in the lower surface of the fluid and the positive 0z axis directed vertically, 
upwards is introduced, the relevant quantities can be written : 


Po= Pog 2—4)+ poxgh(2*—@*), . « « « « (1.8) 
ine Oe ee are eee sn eres) cee t= 220. (19) 


where f is the negative temperature gradient, and po is the density at’ 
the lower surface of the fluid, where the temperature is considered as 
zero. 

502 
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If the static terms are removed, the equations determining the first 
approximation to the velocity, and the first order corrections p, and T’ 
to the pressure and temperature are 


0=—yp—porgl’ +5 (jx Hy) - 24% Sa 


~ Aw=K V2" 9. hag: ae 


to which must be added Maxwell’s equations for the approximations to 
the field, and the equation of continuity for velocities. 


§ 2. DEDUCTION OF THE FUNDAMENTAL DIFFERENTIAL EQUATION. 


From (1.10) and (1.11), the static form of Maxwell’s equations, and the 
continuity equation it is possible to deduce a differential equation for 
T’. This is done by first applying the operator curl to (1.10) obtaining 


0=—curl (pyxgT’) +£ curl (j X.Hy)5 9-85 | 
The second term may be developed as 
(H, . V)j—H,(V - Jj) 
which, in view of (1.5) and (1.6) becomes 
(H, . V)o(E + Vx H,). 


When this replacement has been made the curl operator may be applied 
to (2.1) yielding 
OL 


0=—curl curl (pyxgT’) + ai 


H,.V | cw E +& curl (Vx Hy) | eae a (7-3 


From this equation curl E vanishes by (1.66), curl (Vx H,) reduces to 
(H, . V)V from (1.7) and curl curl (p)xgT’) can be expanded to leave 


2 
0=V*(pg2gT’)—grad (pyzg. VT’) += (Hy- VV. . . (2.3) 


Since p)xg is a constant vector directed along the negative 0z axis, 
the 0z component of (2.3) may be written 


C2 , 2 
— poxg (v-s)T + =5 (Hy. V)w=0. sje en tay 


Finally, w can be eliminated with the help of (1.11) leaving the desired 
equation 


o2 < ane 
Po%Y (v" — zs) fh +. Sy Sane (H, . V4 0), . at ee (2.5) 
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It proves convenient to write this in dimensionless form introducing as. 
a unit of length the depth d of the fluid, and writing the inducing field 
H,=HU where U is a unit vector. Then (2.5) becomes 


2 
| w -ViPV2+A (v* —_ =) | ia Ose her a ete ee ON 
where W==oe oe int A ye nso MOA 


This differential equation with its boundary conditions forms an 
eigenvalue problem for A and its solution determines the minimum 
temperature gradient required for the onset of any convective mode, 
and the form of the temperature perturbation produced. From this, 
the remaining variables can be determined. 


§ 3. SOLUTION TO THE FUNDAMENTAL DIFFERENTIAL EQUATION. 


The solution ot (2.6) is obtained by the method used by Rayleigh for 
the analogous problem in the absence of a field. This involves considering, 
not the close packed hexagonal cells observed experimentally but 
hypothetical rectangular cells. If the plane layer is infinite we write 


T’=T'(z) exp (sdx+imy)+complex conjugate . . . (3.1) 


where zd/l, wd/m are the dimensions of the cell considered and x and y 
are coordinates measured horizontally. 

Suitable boundary conditions are those suggested by Rayleigh (1917), 
the vanishing of the temperature perturbation T’, and of the vertical 
component of velocity w, corresponding to rigid conducting walls. 

(1.11) implies that these boundary conditions are 

a7 (0). d?T’ (1) 


ee Sor) 


FT dat dz? 
Since (2.6) is homogeneous with constant coefficients we may write 
T’(z)~exp isz 
where s is a root of the quartic, 
(s?--/?-+-m?)(s sin 6+1 cos ¢)?=A(l?+m?). . . - (3.2) 


Here, the unit vector U=(cos ¢, 0, sin 4) making ¢ the inclination of 
the magnetic field H, to the horizontal. For ¢=7/2 the field is directed 
vertically upwards and for ¢=0 along the positive OX axis. 

Any solution is a linear combination of four fundamental solutions 
involving the roots sj, 52, $3, and 84, of (3.2), 


4 
T'(z)= & A, exp 18,2. 


j= 
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The boundary conditions put four homogeneous conditions on the 


constants A., whose consistency restricts the roots of (3.2) to those 


yp? 


satisfying 
\ I ] 1 
Fie fe ooo OE 
Atel one : ae “3% =0. . (3.3) 
exp 78, eXP 185 EXP 283 eXP 18,4 


82 exp tS, S_2eXPtS_ Ss2eXPts, 8,7 EXP 18, 


Equations (3.3) and (3.2) define A as a function of the parameters /, m. 
The minimum value of A(l,m) determines the temperature gradient 
needed to start convection, and the corresponding values of / and m 
determine the nature of the initial motion. The good agreement between 
Rayleigh’s and Jeffreys’ results and the measurements of Bénard and 
Chandra (1938) show that this cell size is preserved for finite values of 
the velocities. 


§ 4. ExIsTENCE OF A MINIMUM TEMPERATURE GRADIENT. 


It is possible to show that unless the field is horizontal (sin d6=0) the 
negative temperature gradient must exceed a certain critical value before 
Bénard—Rayleigh convection can occur. With A=0 (3.2) becomes 


(s?-+- 2+ m?)(s sin 6-+1 cos $)2=0 Sg, a aye ae 
which for sin 60 has roots, 
8y=14/(?-+m?), Sy=—1/(?+m?), S3=S,=—I cot ¢. 


Since the characteristic equation has a repeated root the solution can 
be written 


T’(z)=A, exp ts,z+ A, exp ts.z+(A3;+A,4z) exp issz 
and (3.3) must be replaced by 


1 1 7) i) 

8," Sa" ORS — 2is, 
A= . =0. (4.2) 

Xp 18, EXP US» EXP 283 EXP 18s 


Sy7expis, sp2expisy 8,2 expis, (S32—2is3) exp is, 
or (28?y?--y4-+-s4) sinh y=(y?+-s?)? sinhy=0 . . . (4.3) 
where y= V (?+-m?) and s=1 cot ¢ are both real. Hence the only possible 
solution is /=m=0; and unless sinfd=0, A and consequently 8 must 
attain some positive value before convection begins. Physically this 


means that convection cannot occur until the temperature gradient 
attains some negative value, as would be expected. However, it is 
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possible to prove a stronger statement than this, namely: for a non- 
horizontal field (640) solutions do not exist unless A exceeds some 
definite positive quantity, and hence that convective motions only occur 
when the negative temperature gradient exceeds a definite critical 
value Bp. 

In order to do this we must show that not only is there no solution for 
A=0 but that it is not possible to obtain solutions for arbitrarily small 
positive values of A. To inspect solutions for small values of A it is easiest 
to expand the roots s in a series of powers of A and consider only the 
first terms. For convenience (3.2) is written 


(ete O")(e=a)*=67b2a et oe wee oe (44) 
where a=!/ cot ¢, b?=/?+m? and c?=A cosec? 6. For c=0 the roots are 
$,=1b), 8,=—ib, s,=s,=—a. 

We expand s as a Taylor’s series in c, thus 
Os 
s=s(0)-+e(5:) +... 


c=0 
and from (4.4) 


Os | cb? 
dc (s+a)[s(s+a)+s?-+b?]’ 
and for s;=—s,=ib, c=0; ds/éc=O0, while for s=—a 


Cs. eo te C 
dc ab? ee sta)’ 


Os ea b 
e = PEP)’ 


To first powers of c the four roots may be written 


: be be 
A aia om —(a+ en s=—(0— Ta) 
end) 
The determinant 4 which must vanish if the boundary conditions are 
to be satisfied may be expanded to read, 


22 
== | b+ (« —_ =m) + b2 («? + aa) | sinh 6 sin eres A (4.6) 
The only zeroes of these yielding non-trivial solutions are 
be 
v (a?-+6") 


. 
or cmun, [(1 + 3) : 


and since a and b are real, c>z7 and A>7? sin? ¢. 

This establishes that solutions do not exist for arbitrary small values 
of A, but, in view of the approximations used cannot be considered as 
establishing the actual minimum value for 4; however it will be shown 
below that this value is in fact a good approximation. 


= ir ; 
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§ 5. SOLUTION FOR SPECIAL CASES. 


For certain orientations of the inducing field Hy it is possible to write 
down exact solutions to the problem. 


(a) Vertical Field. 
(3.2) becomes 


(s?--[2--m?)s?=A(l?-+-m?) . . . . . « (8.1) 


pee 4 1/2 
with roots $= 4° a {/ 1+ a) -- i} ] 
2+-m?* 4A 12 
=+:i[ ee Witee = hala oleae 


and (3.3) may be written 


_, Pe+m? 4A 12 
ont et es) 
2 2 1/2 


Therefore for a non-trivial solution 


A= (1+ ata). - oo ae be aoe oes 


This has its minimum value 7? at 1/(/2-+-m?)=0 which implies that in 
at least one direction the convection cells must be very small; the size 
in the other direction is indeterminate. 

The critical gradient needed for convection is from (2.7) 
a KO on*H? 


Pu B pg  ¢ 


. (5.5) 


For thermal convection in a viscous fluid with the boundary conditions 
imposed here, Rayleigh (1917) found the critical gradient 


Pat rated 0% . 
= ee a pong n (5.6) 
where »=coefficient of viscosity. 
This suggests introducing a ‘‘ magnetic-viscosity ” 
4 d*® op*H? 
ta = oe a ae (5.7) 


which depends on the scale of the system. This need not be small; in a 
one-centimetre layer of mercury a field of 300 gauss produces one centipoise 
magnetic viscosity. Hence, for conducting fluids magnetic forces should 
be of considerable importance in the process of thermal convection. 
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(b) Horizontal Field, Convection Rolls. 
If the field is horizontal (3.2) becomes 


(sealer (L2 on? )et ae ane, ee ea (58) 


or s= + (4 (e-) a (Pm) : 


The boundary conditions are satisfied if s=-z, ice. if 


A= (oe ati), pee eee TY, 


which is zero for J=0; hence in a horizontal field there is a stable mode 
of convection for arbitrarily small gradients, in the form of rolls with their 
axes along the field direction. This type of motion has been observed in 
the absence of a magnetic field by Chandra (1938) when convection 
occurs in a gas in non-uniform horizontal motion. 


(c) Convection Rolls in an Arbitrary Field. 

There is an exact solution to the fundamental equation corresponding 
to Chandra rolls with their axes along the direction of the horizontal 
component of the field. This is obtained by setting /=0 in (3.2) which 
becomes 

SIN ¢S2(82 =e A ee = (5.10) 
with roots 


Sie MN i 
is LJ (+ Font a waa) | % 


As in (a) the boundary conditions cane 
tf 2 
A=n? sin? (1 +3) aerate ne oan 5c 11) 


which has its minimum A=7? sin? ¢ at infinite m, yielding two finite 
values for s. This value of A is not, of course, necessarily the minimum 
possible one as / has been arbitrarily set equal to zero. 


§ 6. APPROXIMATE SOLUTIONS FOR ARBITRARY FIELD. 


Although the special solutions described above were easily obtained 
only an approximate solution has been secured for an arbitrary field. 
One obvious approximation procedure is to use the solution for a vertical 
field as starting points for an expansion in powers of / cot ¢. 

For sin 640, (3.2) can be written again as 


(8? =202)(S1-@)" = CU ees es ws te (6.1) 
This equation has at least two real roots, the intersections of the curves 
cb? 


y=s?+b? and C— (s-a)2” 
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The four roots can be written s,, 8, e+ty, x+iy, where all quantities 
are real with the exception of y which may be purely imaginary. 
The determinant 4 becomes 


A= 2xy(s,2—s2)| cos (a+s,+8,)—cosh Yy COS a 


+ [laP-+y2)2+9,26,2—(6,2-+5,9)(22—y?)] sinh y sin S29) (6,2) 


From (6.1) and (6.2) roots s can be eliminated leaving ¢ as a function 
of a and 6, which can then be minimized to determine the critical gradient 
and the conditions at the onset of convection. In fact, however, algebraic 
complexity prohibits this programme, and an expansion in powers of 
a about a=0 must be used, by writing : 

a? ds” 


Pe sh aaah 
s=s(0)+a = (0) +5 55 (0+... 


2 922 
c2—=c2(0) +555 (04... fo 


0d ; 
A=A(0)+a>—(0)+.... 


Since 4 is not an explicit function of a we must use 


04 ‘=*404 Os; 
065 a 08.00. 
It will be noticed that dc?/da~dA/da has been equated to zero, the 
condition for a maximum or minimum of A with respect to a. 
The values of the partial derivatives of 4 are calculated in terms of 
6 and c with aid of the known values of s;,(0). Requiring that successive 
terms in the expansion for 4 vanish yields a system of equations for c? 


and its derivatives, and hence to an expression for A(a,b) which to 
second order in a@ is 


A ti 1 m?+b? \2 m+ 5? 
saya fee Vahl sible sa 
ee a ¢ = 7) +B al 4 (Se) sh 3(So) | Pree 5, 
Since a and 6 are necessarily real, this has its least value at a=0, that 
is, at [=0, corresponding to the convection rolls described in § 5. A slight 
inclination of the field therefore removes the ambiguity encountered in 
§ 5 (a). 
The corresponding temperature gradient is 
ao er Oia 
By = Wan a 2 
PoX7, © 


This suggests that only the vertical component of the field should be 
included in the definition of the “ magnetic viscosity ”’. 


Sin? Jy vev ko eee ai OE) 
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§ 7. THe Errect oF Viscosity. 


It will be noted that in all the solutions discussed the minimum eigen- 
value A occurs when one of the cell parameters becomes infinite, corres- 
ponding to a cell of vanishing dimensions. ° This obviously arises from 
neglect of viscosity, which leaves the possible velocity shears unlimited, 
and in order to obtain physically meaningful flow patterns it is necessary 
to consider viscous effects. The viscous stress adds to the hydrodynamic 
equation of motion a force Fy=—» curl curl V which for an incompressible 
fluid reduces to 

Fn=nV’?V. eR Sh et (tome eo A 


This addition calls for no change in the procedure used to get the 
fundamental equation (2.6) which is now modified to read 
Pepe 
4(v*— 53) TN Ue wav jo? beady ee Orn de fess (722) 


2 7 
c 27 : 
where il 21 sin? Q. 


ae oH 4” ny 

This equation may be solved as was (2.6), and if the same boundary 
conditions are used (although these have admittedly become somewhat 
artificial) the solution presents no great difficulties for the special cases 
studied in § 5. In particular, for convection rolls of width wd/m 


A qr art m\3 
aati i aan ro(+3) ee eS ar (ase) 


Minimizing this yields an equation determining m which is most neatly 
written by introducing a variable y such that 
72+ m= 3/27?y 


whereupon the condition that A be a minimum becomes 


Nile eLearn mecmne) 0 fo 2i-4(°7 14) 
7 ; 
The corresponding critical temperature gradient is 
27. VAS are 
= i (* Awe ga ed Ie aS Wess) 
errern) pcre u ttl x°B, + Bax (7.5) 


For 7q=0, x= 1 and f reduces to the value given by Rayleigh. Identical 
results holdin a vertical field with m? replaced by b?=(?-+-m?. 
. The width f of the roll is given by 


pee ee ive) 


For square cells occurring in a vertical field this is multiplied by 1/2. 
The formula determining the width’f of the cell renders (7.5) intelligible. 
The effect of the magnetic field is to establish a resistance to flow across 
the lines of force, hence the cell width is reduced, thus diminishing the 
excursion of any element of fluid across the field lines. This squeezing up 
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of the cells increases the velocity gradients, hence increasing the viscous 
stresses. The effect of a field then, is to add a term By due to the magnetic 
viscosity, and a multiple y2>1 of the normal viscous term f,,. 

The effects here described are of sufficient magnitude to be easily 
detected. For example, for a layer of mercury one centimetre in depth 
placed in a vertical field of 1000 gauss the relevant quantities are 


n=0-014 poise, 8, =0-162° C./cm., 


y= 0-152 poise, B=2-86° C./em., 
y=2°6, f=0-69 cm., 
f,=1-41 cm. 


§ 8. THe NATURE OF THE SOLUTIONS. 

Once the cell parameters have been fixed the remaining dynamical 
quantities can be determined in terms of the temperature perturbation T’. 
Unfortunately, the theory is linear and homogeneous in the dependent 
variable, and the solution is only obtained to within a multiplicative 
factor. This means that in common with existing analytical theories of 
thermal convection it does not predict the heat transported. 

If now, we remove the convention that displacement be measured in 
units of d, the depth of the layer, the temperature perturbation for the 
two-dimensional roll solution becomes 


ee my . 8% 
T’=A cos 7 *@q- (8.1) 
From (1.11) the 0z component of velocity can be written 
K K my SZ 
Bye as ea ee abe Dace, BIS 2 
w BY Ba he Og are «aoe Bee 


The second velocity component v is determined from the equation of 
continuity as 
7 7 C8 s- (8.3) 
The electric charge density ¢ is obtained from the continuity equation 
for the current 


K m: Zz 
=— 5, Bd? (m2-+82)A sin —” cos 


=v j =o (aiv E--div - Vx H) : 
whence 
SAE d eo SE : a Eu Ow Ov 
{[=7- div E= vars H. curl V=— aie H cos ¢ & = =) oF Le eee) 
The electric field can be written 


Hw Hv 
E= (0 i cos ¢, — cos $) <= sag ac ee 


r . 4 . . bs . . 

This satisfies E,,,—0 at the boundaries which are assumed perfectly 
conducting, although the normal component does not vanish but gives 
a surface charge Y 

ew K m'+s8" 8 


Te ayaa Heos¢.Asin my. ©... > (8:6) 
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The current density j is given when E and V are specifid as 


: le ye 
j=(% sing.) 0,0). PO geo S77) 


The induced magnetic field can now be determined since the vector 
potential within the fluid, a solution of \A=(47/c)j can be written 


Samy SZ 
SN (2S Veen ree Be Be 
A ( aa a ee 7 cos 7, 0,0), 


and the field h within the fluid 


4nuoH Ks? | |. my . sz 4apoH K my 82 : 
=| 0, ———- = — — —-, —-_ > Ss — cos — 8 
h (0, Ol Bi daa od F&A 008 “H cos =) (8.8) 
Outside the fluid h is most easily determined by noting that the 
boundary conditions at the surface are satisfied if h—=—\/y where wp is 
the magnetostatic potential due to a surface pole distribution 
po HK my 
oe d B sA cos % 


on the surface z=0 and a similar distribution with opposite sign on z=d. 
The field then has the form 


n= (0, dais aN sin —? exp Bad pee a cos xp 7) 


CBR d iit eenicied ee Cdwutt 
(8.9) 
for z<0 and a corresponding quantity varying as exp m ( 1 — 5) for z>d. 


Finally, the pressure perturbation is given from the hydrodynamical 
equations 


g 
sin? ¢ . v, ene: —pyxgT’. 


Oz 


The integrebility condition is seen to be equivalent to the equation for 
T’ and the pressure perturbation may be written 


d m Sz 
p'=poxg — . A cos * cos >... ee eee ee S10) 


§ 9. StaBiLiry CONDITIONS. 


An important assumption has yet to be considered, the condition for 
the existence of stable convection currents. The argument used by 
Jeffreys to justify treating the time independent problem was as follows. 
Consider a solution with time dependence in the form exp yt, where y is 
a function of 8 the temperature gradient. As f is increased y moves from 
negative to positive values, negative values corresponding to static 
equilibrium and positive ones to instability and the establishing of 
convection currents. The condition for marginal stability, the case 
separating these two regions, is that y vanish. It is however, necessary 
that y be real, and Jeffreys appeals to experiment to show that undamped 
convective oscillations do not occur as they would if y had a pure imaginary 
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value at marginal stability. Such an appeal cannot be made for the 
magneto-hydrodynamic problem, and it has been suggested (Walen 1944) 
that oscillations do occur, therefore to determine the significance of the 
solutions for marginal stability it is necessary to examine the time 
dependent problem. We shall seek a solution to the time dependent 
problem in the form f(z, y, z) exp yt where y will be found to satisfy a 
cubic equation with coefficients depending on the constants of the material, 
and flow pattern, and the negative temperature gradient 8. For low 
enough values of B, y has a negative real part indicating that small 
motions are damped out, and the conduction of heat through the resting 
fluid is a stable process. If f is larger the real part of y is positive and any 
motion increases with time until limited by non-linear effects. We will 
find that y is real if this latter condition holds unless.the constants of the 
material satisfy very stringent conditions, and therefore that steady 
currents and not oscillations are to be expected. 

Again, the linearizing procedure is carried out, but the time independent 
assumption is not made. The first order equations can be written to 
determine the significance of the solutions for marginal stability, 


OV) tent ; 
E75, VEO g + IX My Re Si 
aT’ 
a 7 WP= KIT, . ot ae hele SN oy fk ina 
div, V=0,"5, ya Ae ne re 


Again, it has been assumed that variations in p are small, and moreover 
that the time variation is slow. Then, the displacement current can be 
omitted from Maxwell’s equations, 


A 
curl h=— j, 0 0 BU. e tine flies an aie 
0h 
curlE=—" ere er ee 
div E=.4rq, *. 2 GA, yee a 
div Rh ==.0) tM nn ye en OY 
ee EN) 
divj= 3. _ a i enact eM 
j=o[E+% vxu]. LE (ora 
From (9.4), (9.5) and (9.9) . 
0 c : 
Ot = 4c Vi ] h=(H4 V)V. . . . . . (9.10) 
And from (9.1) 
0 ' 
ai (VV) + 2 (Hy - V)V?A=—curl curl («T’g). . (9.11) 
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Using (9.10) to eliminate A, and (9.2) to eliminate the 0z component 
of velocity leaves the fundamental equations for T’ in the form 


G-ie 9 Ke c. aq ae ry o(2 oe) cep 
ot Atrpuo Ot ap rr ee ni2¥) Casa 


rs) 2 2 
— 08 (5 : v*)(v?— 53) T'=0. ewe (Oa 2) 


ce Amps 
If we consider a vertical field and use the Rayleigh boundary conditions 
the appropriate solution can be written as 
T’~exp (ila+imy-+yt) sin sz 
where the parameters satisfy the characteristic equation 


2 Kerk 
3 ke Kis) Je2y2 [ ie 2o2]2__ poe 6D 
say (— - oa Air poy ee ae Atruo 
ie 2K c2p4 __ c say ee 
+ gio HK sti — 7 agh(kt—st)k*=0, . . . (9.13) 


where k?=/?+-m?-+s?. This equation can be written 
y+ Ay? +By+C=0, 
and in order. that steady convection currents appear C<0. 


If the equation is to have a complex root with positive real part 2, 
a condition for instability to oscillations, must satisfy 


x3+ 2Ax?+ (A?+B)x+(AB—C)=0 
hence the condition that oscillations, not steady motion should occur is 
C>0>AB—C. 
This in turn, is a condition on the negative temperature gradient 
co? Hstk2  k6e2 > ce? 
a —+tkK , 14 
~ Tént Poko —§ 4mpuc E= er ] ; Ole 


and the condition for oscillation 


2 
eat KH2k2s2 >agB(k2 —s?) 


ct c2 \2 CK 7] poe? Te 
—— aaa Kors Jee (9.16 
Seam 167? Ag | (za) v 4c _| H? ar ee) 
Hence a necessary condition on the material involved is 
ro 
— =< l. Pe aS cos a (OGL 
4auKo i ( ) 


The values of , K, o for solid copper give 6=0-7, hence the condition 
is not easily satisfied in the laboratory. 
A necessary condition on @ for oscillation is obtained by substitution 


in (9.15) 
eS eae 

re PAG Cal ees rey . 5 . 5 : 6 ( . 7) 
Tf such oscillations do occur it is possible also to determine the gradient 

required for the onset of convection, and the frequency of the oscillations 


at marginal stability, 


Cre eee 8 ec? | gos 4 
8 eee posi ER eet be Riel eae ad BL be 
#9Po 167? pKo ne 4 myo |_4npo 4 ( ) 
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which with (9.14) gives for the cell parameter k?=3/2ys? and x a root of 
So rl ase 


OX) or Bape GT (9.19) 


The corresponding frequency obtained from (9.13) is 


tu? 1-9 8my?s? 5 


Introducing the magnetic energy per unit mass &=(jH}?)/(87p) and 


s=7/d this is 
7 1—6é , 1 ; 
v= /(2- 1) 61 (ae 6 —Kx” =) cate C (D2) 


§ 10. CONCLUSION. 


We have examined the modifications produced in the Rayleigh-Jeffreys 
theory of slow thermal convection by the addition of magnetic forces 
acting on induced currents. It has been shown that for most materials 
the flow is steady, although oscillations seem possible in principle. An 
approximate value for the critical gradient has been found, but as was 
the case for Rayleigh’s analysis this is probably somewhat in error due 
to improper choice of boundary conditions. A numerical analysis similar 
to that given by Jeffreys (1926) for ordinary convection would improve 
the estimate of A. On the other hand, the present analysis is a great deal 
simplified and the proper dependence of on depth and field strength has 
probably been deduced. This predicted dependence could be tested 
experimentally by determining the temperature gradient for which the 
thermal conductivity law ceases to hold for a thin layer of fluid heated 
from below as was done by Schmidt and Saunders (1938) for ordinary 
convection in water. 

Finally, the author wishes to express his thanks to Dr. E. C. Bullard 
who suggested the problem and gave valuable advice and encouragement. 
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CXLIV. The Boundary Layer in Three Dimensional Flow.—Part II. 
The Flow near a Stagnation Point*. 


By L. Howarrna, F.B.S., 
Department of Mathematics, University of Bristol +. 


[Received August 16, 1951.] 


SUMMARY. 


The equations of boundary layer flow in the vicinity of a stagnation 
point on a general (three dimensional) surface are discussed and shown 
to be reducible to a pair of simultaneous ordinary third-order differential 
equations containing a single parameter c which is determined by the 
mainstream flow. The variation of c can be effectively limited to the 
range from 0 (corresponding to two dimensional flow) to 1 (corresponding 
to the axial flow past a body of revolution), and solutions have been 
computed for the cases c=0-25, 0-50, 0-75 and are tabulated below. 
A series expansion useful for small c is also given. 


INTRODUCTION. 


In Part I. (Howarth 1951) the equations of flow in the boundary layer 
along a general curved surface were derived using an orthogonal coordinate 
system £, 7, ¢ comprising respectively the surfaces generated by the 
normals along the lines of curvature and the parallel surfaces. The 
equations for steady flow were found to be 


Pcie EE Gl 
Bee ic —K,w=— eee (2) 
3" neo ® Eee (ate a Pret BF) 2°73 

where h,, hz, hz are the length elements and 
ee 


The yuanet’ of this paper is to discuss the flow in the neighbourhood 
of a stagnation point S on the surface. We shall suppose that S is the 
oint €=0, 7=0, C=0, that the surface is regular at S and in particular 


P 
that (2), (= 3) ea (32), (=) (Ze : are bounded. We 


* This paper forms part of one of the essays for which the Adams Prize 
was awarded in 1951. 
+ Communicated by the Author. 
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may, without loss of generality, suppose that é and 7 are so defined in 
magnitude that at S hj =h.=1. Then, in the vicinity of 8, 


n=1+ ($2), e+ (F), nt (52), 5 eae 
mu BeBe Be 
approximately. 


The Nature of the Solution. 

We shall suppose that the flow in the mainstream outside the boundary 
layer is irrotational, so that, if we denote by U, V its component velocities 
in the directions €, 7 increasing respectively, 

ov 0U oh, oh 

howe rue, Ae ya 

At S,. by definition, U and V vanish and since § is a regular point of the 
surface, we shall have in its vicinity 

Us At Baye ae ct Rie ae (8) 

V=DE“Cam Ss. 6 Re | 

where A, B, C, D are constants. Hence we see from the leading terms 

n (7) that 


=0*.;_ «i, ee 


D=B: 


Since the velocities in the boundary layer will be dominated by U, V, 
it is evident that the curvature terms in equations (1), (2), (3) may be 
neglected and that h, and h, may each be replaced by unity in determining 
the leading terms near Sj. As we might have anticipated on physical 
grounds, we see that the equations of flow and the condition of 
irrotationality of the main flow are to this order exactly the same as 
if the surface were plane, ¢ measured normal to it and € and 7 along two 
axes at right angles in the surface. We are at this stage, therefore, free to 
rotate our €, 7 axes about Sf to any convenient position whilst at the 
same time maintaining the form of our equations. 

Now, by the irrotational condition, the tensor 


|e), (Fe) 
(i), (er), 


is symmetric and by suitable rotation of the axes Sé, Sy to Sa, Sy, say, 


we may express it in its principal form. In other words we may take 
the mainflow outside the boundary layer to have components (ax, by) 


* It may perhaps be mentioned here that, just as in the corresponding 
two dimensional problem, we do not suppose that U, V in themselves satisfy 
the equation of continuity, but only after due account is taken of the outflow 
from the boundary layer. 

+ These assertions may be justified a posteriori. 
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and proceed with the solution as an immediate extension of the standard 
two dimensional solution. We may assume, without loss of generality. 
that @ is positive. 

We find, now using w, v, w for the velocity components referred to our 
new axes S(a, y, C), that 


U—=ON (2) 
ee RB ee tee eee a (LO) 
UO va lafle )+-bg(z) 
ee qu2g 
where Ce ie ee ag or ee (11) 
and Fe ete eee ee ees (19) 
1929 "” 1 ” 1 Uda 
g*—gq eens Tee ee A ks) Pip ears er ee ag 0 esd) 
and EAU, OU ae te ae ee Seer ag eS 


The boundary conditions for (12) and (13) are 
f=g=f’=g'=0 at 2=0, 
Yh g >1l as z>0. 


This slightly asymmetric form of the equations is advantageous for 
discussing conditions when c is’small. 
The skin friction components 7;,, T,, are given by 


cop = pv 2a? *xf"(0) 
t= = pvl2q3/2cyq" (0 , 
The following special cases aie previously been discussed :— 
(i.) c=0, which corresponds to two dimensional flow, 


(ii.) c=1, which corresponds to flow round a body of revolution 
placed symmetrically in a stream. Here we have obviously 
g=f, so that 


Fee feet egress ee) (15) 


and this is Homann’s (1936) equation. 
We may confine attention to the range c=0 to c=1 since it can readily 
be verified that 


EEE Py f(z, —=f'(, ¢), (16) 


J\e; —c)=—g(2, c); g' (2, —c)=—9'(z, c) 


(2.4) =et%9 (ste ¢): r(a2)=0 Ge): | . (17) 
a(« 2) Sal Gee c), a (22) =F (Ga) 


In any event it is simply a matter of notation to ensure that in any given 
case 0O<c<l. 


and 


5D2 


1436 L. Howarth on the 


Numerical solution. 

The solutions of equations (12) and (13) may be obtained by standard 
numerical methods, but before embarking on them it is convenient to 
have expansions when c is small, for it is here that conditions change most 
rapidly with c. To obtain these expansions we write 


f=fotefitefet 
J=JotCJotegr+.--5 j 


substitute in equations (12) and (13), equate coefficient of powers of c, and 
obtain the following equations :— 


(18) 


foe thfo— Hercule tee tS Sao ete ey ele 
1 +hofi— 2fofittoii=—Golos = ee 
fo t+fhofe—2fofetto fe=fP —f Sit Osi too eee 
CS Fi I re See. MP Re Qs 
1 +honii= ft gage ; (23) 
92 +ho92= fo —fiGi + 29091 -IJoI iJon» + + (24) 


The boundary conditions are f,=9,=f;=g,;=0 at z=0 (r=0, 1, 2,...) and 
foHgoeal. ~_fiSe 0 (ra es 


Equation (19) is the standard two dimensional equation (see Howarth 
1934 for the solution) and equation (22) has previously been integrated by 
Prandtl (1945) and by Sears (1948). Equations (20), (21), (23), (24) have 
been integrated numerically and the results for the first derivatives are given 
in Table I. They enable velocity distributions to be evaluated for small 
values of c. For completeness, values of fj and gj have been included in 
the table. The corresponding results for the skin friction components 
are determined by 


”(0)=1-2326, —f%(0)=0-0478, —f(0)= 0-059, 
gi(0)=0:5705, g%(0)= 10535, —-gf(0) = — 0-537. 


These results, although insufficient for accurate computation as far as 
c=1, were found to give, in conjunction with the known solution when 
c=1, useful starting points for ‘the detailed numerical solution of 
equations (12) and (13), which has been undertaken for c=0-25, 0-50, 0°75 
by the Adams method. The results for the first derivatives are shown in 
‘Table II. and include the results for c=1 computed in the course of an 
interesting investigation by Miss Hannah (1947) from Homann’s equation. 
(The values of fj, gj in Table I. are, of course, the results for c=0.) The 
skin friction is determined by the values in Table ITI., and a few values 
of f and g are contained in Table IV. 

It may perhaps be remarked that, throughout the range c=0 to 1, 
jf varies very slowly—in fact f’(0) only changes from 1-233 to 1-312— 
whereas in g the changes are much more rapid, the corresponding 


— ” 
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TABLE I. 


0-0 0-000 0-000 
0-1 0-118 0-005 
0-2 0-227 0-010 
0:3 0-325 0-014. 
0-4 0-415 0-019 
0:5 0-495 0-023 
0-6 0-566 0-027 
0-7 0-630 0-030 
0-8 0-686 0-033 
0-9 0-735 0-035 
1-0 0-778 0-036 
1-1 0-815 0-037 
ey 0-847 0-037 
he33 0-874 0-036 
1-4 0-897 0-035 
1-5 0-916 0-033 
1-6 0-932 0-031 
ley 0-946 0-029 
1-8 | 0-957 0-026 
1-9 0-966 0-024 
2-0 0-973 0-021 
2-1 0-979 0-019 
2-2 0-984. 0-016 
2:3 0-988 0-014 
2-4 0-991 0-012 
2-5 0-993 0-010 
2-6 0-995 0-008 
2-7 0-996 0-007 
2-8 0-997 0-005 
2-9 0-998 0-004 
3-0 0-998 0-003 
3-1 0-999 0-003 
3°2 0-999 0-002 
3°3 0-999 0-002 
3-4 1-000 0-001 
3°5 0-001 
3:6 0-001 
3-7 0-000 
3:8 

3-9 

4-0 


0-000 0-000 
0-006 0-057 
0-012 0-114 
0-018 0-171 


0-031 0-339 
0-035 0-393 
0-037 0-446 
0-038 0-498 
0-039 0-547 
0-038 0-594 
0-036 0-639 
0-034 0-681 
0-032 0-720 
0-028 0-756 
0-025 0-789 
0-022 0-819 
0-018 0-846 
0-015 0-870 
0-012 0-891 
0-009 0-910 
0-007 0-926 
0-005 0-939 
0-003 0-951 
0-002 0-961 
0-001 0-969 
0-001 0-975 
0-000 0-981 
0-985 
0-988 
0-991 
0-993 
0-995 
0-996 
0-99 
0-998 
0-999 
0-999 
0-999 
1-000 


0-488 
0-517 
0-536 
0-546 
0-547 
0-540 
0-527 
0-507 
0-481 
0-452 
0-420 
0-385 
0-350 
0-314 
0-279 
0-246 
0-214 
0-184 
0-157 
0-132 
0-111 
0-092 
0-075 
0-061 
0-049 
0-039 
0-031 
0-024 
0-018 
0-014 
0-011 
0-008 
0-006 
0-004 
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variation being from 0-570 to 1-312. These changes are reflected in the 
series solution in the much larger coefficients, after the first, in the g series 


as compared with the f. 


CONCLUSIONS. 


The main conclusions are embodied in Tables I. to IV., but some of the 


general trends may, however, be worthy of special note. 
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TABLE II. 

c=0-25 c=0-50 c—0:75 e=1-0 

2 t: , g’ i g ft g’ Ts = 

0-0 0-000 0-000 0-000 0-000 0-000 0-000 0-000 
0-1 0-120 0-079 0-122 0-097 0-124 0-113 0-126 
0-2 0-230 0-156 0-233 0-190 0-238 0-218 0-242 
0:3 0-330 0-230 0-335 0-277 0-342 0-315 0-349 


0-420 0-301 0-428 0-359 0-436 0-405 0-445 
0-502 0-369 0-511 0-435 0-521 0-488 0-532 
0-574 0-434 0-585 0-506 0-596 0-562 0-609 
0-639 0-496 0-650 0-572 0-663 0-630 0-677 
0-696 0-554 0-708 0-631 0-722 0-690 0-736 
0-746 0-608 0-758 0-686 0-772 0-742 0-787 
0-789 0-658 0-802 0-734 0-816 0-789 0-830 
0-826 0-704 0-839 0-777 0-852 0-828 0-865 
0-858 0-746 0-870 0-815 0-883 0-862 0-896 
0-884 0-784 0-896 0-848 0-909 0-891 0-920 
0-907 0-817 0-918 0-877 0-929 0-914 0-940 
0-926 0-847 0-936 0-901 0-946 0-934 0-955 
0-941 0-873 0-951 0-921 0-959 0-949 0-967 
0-9F4 0-896 0-962 0-938 0-970 0-962 0-975 - 
0-964 0-915 0-971 0-952 0-978 0-971 0-983 
0:972 0-932 0-979 0-963 0-984 0-979 0-988 
0-979 0-946 0-984 0-972 0-988 0-984 0-992 
0-984 0-957 0-988 0-979 0-992 0-989 0-995 
0-988 0-966 0-992 0-984 0-994 0-992 0-997 
0-991 0-974 0-994 0-988 0-996 0-994 0-998 
0-993 0-980 0-996 0-991 0-997 0-996 0-999 
0-995 0-985 0-997 0-994 0-998 0-997 0-999 
0-996 0-988 0-998 0-996 0-998 0-998 

0-997 0-991 0-998 0-997 0-999 0-998 

0-998 0-993 0-999 0-998 0-999 0-999 

0-998 0-995 0-999 0-998 0-999 0-999 

0-999 0-997 0-999 0-999 1-000 0-999 


WNNNNNNNNNN YE SE Pe Eee ooooo°o: 
COWMDAADTNLPWNHCHOWDAAGATIEWNH CHO WMDATGIS 


Taste ITI. 


0-00 0-25 0-50 0-75 


1-233 1-247 1-267 1-288 
0-570 0-805 0-998 1-164 


As a matter of notation, since we have restricted attention effectively 
to 0<c<1, we may call the axes Sx, Sy respectively, the major and minor 
principal axes of the mainflow. We may remark first of all that the 
changes with c in this range are then relatively small for the velocity 
component in the major direction, though marked in the minor direction. 


——— 
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Secondly, the direction of the velocity vector at a particular station 
(that is for given x, y) may undergo appreciable changes as we move 
through the layer. The velocity is zero at the surface, but there is a 
limiting direction of flow there (which is also the direction of the resultant 
skin friction), and this direction is inclined to the mainflow direction at 


an angle 
1 (9 (9)\ cy cy 
tan! ee ee 
an ( 70)) tan : 


it appears from an examination of Table II. that, as we pass from the 
mainstream towards the surface, the angular deviation of the velocity 
vector steadily changes towards this value, the sense of the change being 
from the mainflow towards the major principal direction. For c=0-00, 
0-25, 0:50, 0-75, 1-00, g"(0)/f’(0) takes the values 0-463, 0-645, 0-788, 
0-903, 1:0, and the maximum values of the corresponding changes in 
direction are then 21°*, 14°, 6°, 3°, 0° respectively; the first four occur 
respectively when cy/x takes the values 1-47*, 1-25, 1:13, 1-05. 


TABLE IV. 


c=0-50 C=00 c=1-00 


fi g A g f=9 
0-000 0-000] 0-000. 0-000} 0-000 
0-138 0-114]0-140 0-130] 0-143 
0-474 0-412| 0-482 0-456] 0-492 
0-913 0-826] 0-928 0-892] 0-944 
1-396 1-207] 1-413 1-374] 1-433 
1-892 1-789|1-912 1-870} 1-932 
2-392 2-289| 2-412 2-370] 2-432 


0-0 
0-5 
1-0 
1:5 
2-0 
2°5 
3-0 
3:5 
4-0 


Thirdly, we may notice that, apart from c= 1, there is a general tendency 
for g to approach the mainstream value less rapidly than /, this tendency 
being most marked at c=0. In fact we might artificially speak of two 
boundary layer thicknesses, one associated with each principal direction, 
though it must be emphasized that they do not differ vastly. There is 
a general tendency for both thicknesses to decrease as c increases towards 1, 
as can be seen from Table II., or from the values of the integrals 


| (=) dz, \ (1- iy) & 
J0 On 0 by 
which are 0:648, 1:026 when c=0; 0-628, 0-829 when c=0-25; 0-608, 


0-711 when c=0:50; 0-588, 0-630 when c=0-75; and 0-568, 0-568 when 
c=1:0. 


* This is the limiting value as c + 0, of course. 
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Finally, it may be remarked that, just as in the corresponding two 
dimensional solution, the solution given here is in fact a solution of the 
full equations of viscous motion in Cartesian coordinates. 
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CXLV. CORRESPONDENCE. 


On the Calculation of Characteristic Temperatures from the Elastic 
Constants. 


By M. Brackman, 
Physics Department, Imperial College, London *. 


[Received September 18, 1951.] 


Rucent work (Leigh 1951) on the free energy of alloys has drawn attention 
to the calculation of 6, values in cases where the solids are elastically 
anisotropic. This anisotropy is due to the smallness of the (C,,—C,,) 
value, where the elastic constants of the (cubic) solid are C,,, C,, and C4, 
respectively. In such cases no simple method exists for calculating 0p. 
One has either to calculate the mean velocity directly, which is a lengthy 
procedure, or to use the Hopf—Lechner (1914) method, which in these 
cases is equally lengthy and tedious. While this is immaterial in single 
cases, the time spent in calculation can become disproportionately large, 
where a series of alloys has to be considered. 


TABLE I. 


Substance Li 


The purpose of this note is to communicate a semi-theoretical formula 
for 6, designed to cover cases where (C,,—(C,,) is relatively small. The 
method used is to consider the shape of the surfaces of constant frequency 
(cf. Blackman 1935). While this does not lead to a definite formula, it 
does suggest a reasonable form ; a constant which then remains is fixed by 
fitting close to the value for Lithium calculated by the standard methods. 

The resultant form is 

| 215 ol 
vps. io . ps2A 


(Cy —Cy9)29(C py + Cy + 2Cg4)2?(Caa)?®, omen) 
where 4=volume of the unit cell, s=number of particles per cell, 
p=density, and v=Debye frequency. The Debye 6, value is equal to 
hyy/k (h=Planck’s constant, k=Boltzmann's constant). 

Table I. shows the agreement between (1), and the values of vy calculated 
by standard methods. These, it should be remembered, are not exact 
calculations and may be several per cent in error. The @ values so 
obtained are labelled 65. 


* Communicated by the Author. 
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The validity of (1) is surprisingly wide, but it should be emphasized that 
it is nevertheless limited. When the above formula is used, the elastic 
constants should be such that the C,,/C,,, Cy2/C,, ratios lie inside the 
region covered by Table I, (cf. Blackman 1938, p. 76). 

As an example, we take the Ag—Au alloys, the elastic constants for 
which have been given by Roehl (1933). The @p values were calculated by 
Roehl not from these values, but from averaged elastic constants. The 
result of direct calculation using (1) is given below. 


Tasce II. 


25 At. 9 5 yh 5 At. % 
Substance Ag 5 Re Yo D0 i He 75 Yo 


rr Au 


Op 203 195 185 172 160 


The conclusion to be drawn from this table is the same as that drawn by 
‘Roehl, namely that 9 is not a linear function of the percentage of gold 
added to silver. 
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Inelastic Scattering of Deuterons. 


By R. Husy and H. C. Nrewns, 
Department of Theoretical Physics, The University of Liverpool *. 


[Received September 19, 1951.] 


THE inelastic scattering of protons has been much studied in experiments, 
the results of which seem understandable in terms of the formation and 
decay of a compound nucleus. 

Recently, evidence has been obtained showing that deuterons also can 
be inelastically scattered by nuclei (Greenlees, Kempton and Rhoderick 
(1949) ; Holt and Young (1949)). In some cases, e. g. the excitation of the 
1-38 MeV. level in Mg™, the total cross section seems much larger, at about 
10-°° cm.’, than would be expected for this process on the theory of 
compound nucleus formation. Further, the angular distribution, which 
has been measured in the above-mentioned case by Holt and Young, is of a 
complicated form, with a strong concentration at rather small angles of 


————— 


* Communicated by the Authors. 
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scattering, such as would seem unlikely on that theory. The present 
authors therefore attempted an explanation through electrical interaction 
of the deuteron with the nucleus, but this was found to yield much too 
small a cross section to account for the present experiments (Huby and 
Newns (1951)). Ina similar treatment, Mullin and Guth (1951) were able 
to obtain an angular distribution of approximately the required shape, 
but to account for the measured total cross section, they would require a 
dipole moment of the order of three times the nuclear radius, which seems 
too large to be physically plausible. 

The purpose of this letter is to propose a very simple mechanism to 
account for the main features of the experimental results. ‘The mechanism 
is suggested by the similarity of the observed angular distributions to those 
obtained in various (d—p) and (d-n) reactions, and it has something in 
common with the stripping theory for the latter. We assume that only 
one of the constituents of the deuteron, say the neutron, interacts with 
the nucleus, the proton being outside the range of the nuclear forces. The 
neutron interacts at the surface of the nucleus, exciting it and then 
“ bouncing off ”’, the deuteron emerging as a whole having been inelasti- 
cally scattered*. A calculation can be made, based on the Born approxi- 
mation, and following closely a treatment of the (d—p) and (d—n) reactions 
of Bhatia and Huang (reported by Huby (1950)). Despite its crudity, the 
latter had considerable success in explaining the observed (d—p) and (d—n) 
angular distributions, and this seems an encouragement for the present 
calculation for (d—d*) reactions. 

The differential cross section can be written 


9 > 2 
o(O= | | UitxtxbWdixixn: do AR, aR, | 


where %, and , are the initial and final deuteron translational wave 
functions, y,; and x; the initial and final nuclear wave functions, and x; 
and yp, the initial and final internal deuteron wave functions. V repre- 
sents the interaction between the neuteron and the nucleus. This can be 
evaluated if we assume (i.) the interaction V operates only on a sphere of 
“‘ radius of interaction ”’ a, (ii.) the incident and scattered deuterons waves 
are plane waves of wave number k, and k, respectively, (iii.) the ground 
‘state of the deuteron is a pure triplet-S state, with radial wave function of 
the form exp {—«|R,—R, |}/| R,,—R, |- 
The result obtained is 


4, k 1/2 2 
o(0)=5| AUP E tant > SS J elk) | 
1 


where k?—=(k,—k,)?+4k,k, sin? $0, and the A;’s are unknown constants. 
The formula for o gives scattering predominantly at small angles. The 
sum is taken over integers / which are related to the angular momentum 


* A rather similar model was in fact proposed by E. Guth (Phys. Rev., 68, 
280 (1945)) in predicting inelastic scattering of deuterons. 
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change of the nucleus. There are selection rules requiring that if J; and 
J, are the angular momentum of the initial and final states of the nucleus, 
then J,=J; Site (but when J;=0, and /=0 then J,=0 or 1), and / is odd 
or even according as there is a change of parity or not. The values of / 
can be determined from the angular distribution. If more than one value 
of 1 can contribute, the lowest one usually predominates. Thus the 
angular distribution of inelastically scattered deuterons can lead to 
information about the spins and parities of nuclear states. 


8 


Angular distribution of inelastically 
scattered deuterons from the 
1-38 MeV. excited state of Mg”. 

Experimental curve 

Theoretical curve — — — — 


o(@) arb. units. 


A particular example is the case of deuterons inelastically scattered from 
the 1-38 MeV. excited state of Mg*4. A theoretical curve for /=2 is shown 
in the figure together with the experimental curve (Holt and Young (1949)). 
Since the ground state of Mg®4 has spin zero and even parity, the first 
excited state has spin 1, 2 or 3 and even parity. 

The parameter a is not well defined. One might expect it to be given 
by the Gamow formula (Gamow and Critchfield (1949)) for the nuclear 

radius, but as in the (d—p) reaction, the values of a needed are always 
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larger than this. In the magnesium case discussed above, we used 
a=7X10-%cm. To make curves for different values of / fit, values of a 
which seem unreasonable would be required. A spin of 2 for the excited 
state would be in agreement with the experiments on the angular corre- 
lation of successive y-ray emissions from the nucleus (Brady and Deutsch 
(1948)). 

The angular distribution of deuterons inelastically scattered from the 
1-66 MeV. excited state of Ne?° has been measured by Middleton and Tai 
(unpublished). This can be fitted by the theoretical curve for J=1. 
Since the ground state of Ne®® has spin zero and even parity, the 1-66 MeV. 
excited state must have spin 0, 1 or 2 and odd parity. 
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Experiments on the Unsaturated Helium II Film. 


By R. Bowers, D. F. Brewer and K. MENDELSSOHN, F.R.S., 
Clarendon Laboratory, Oxford *. 


[Received September 17, 1951.] 


Tue transfer of the helium IT film at pressures lower than the vapour 
pressure of the liquid has been investigated by a method which consists 
essentially in measuring the heat conduction of a tube containing He 
gas under variable percentage saturation. 

The principle of the apparatus is shown in fig. 1. The bottom of the 
(thin-walled german silver) tube is connected to a constant-temperature 
heat sink, while the top contains a heater and a thermometer. When, 
under fixed conditions of temperature of the sink and of pressure in the 
tube, heat is supplied at increasing rates Q, the thermometer remains 
constant at its original reading until a certain critical heat input Qe is 
reached, when it rises sharply. The constancy of temperature is main- 
tained by flow of the film up the inside wall of the tube, the film evapor- 
ating at the top, and the helium returning as gas to the bottom of the tube. 
This circulation provides a powerful means of heat transport, in which the 
latent heat of the helium is taken up by evaporation at the top, and given 


* Communicated by the Authors. 
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off by condensation at the bottom to the heat sink. Incidentally, the heat 
of evaporation from the film is slightly larger than that of the bulk liquid 
(L,,), since it also contains the energy of excitation of the superfluid into 
the normal state (7148). The circulation process must break down when 
a film flow rate larger than the critical rate 7, is required to maintain it. 
The critical heat input Q, thus provides an accurate means of determining 
the critical rate of mass transport of the He IT film since 
fea doer 
° (Lpgt+TAs) | 


Higa: 


HEATER 


THERMOMETER 


HEAT SINK 


Such factors as might influence the validity of this equation, e. g. viscous 
drag in the returning gas, heat flow by conduction along the metal tube 
and through the gas, and heat resistance of the partition between the tube 
and the heat sink, were negligibly small in our case. A check on the 
validity of the method was provided by plotting the critical flow rates 
obtained at various temperatures under saturation conditions (7,) against 
the absolute temperature. The values derived are in satisfactory agree- 
ment with the shape of the curve of. transfer rate out of a beaker 
(Mendelssohn and White 1950). 

Transfer rates under various percentage saturations have been obtained 
for four different temperatures, and the results are given in fig. 2, where 
the reduced transfer rates7,/r, are plotted against the percentage saturation. 


— 


— 
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The method is sensitive enough to detect mass flow of the order of 0-1 per 
cent of the saturation flow rates (r,/r,~0-001). 

The results show clearly that any appreciable film transfer will only 
occur hear saturation. In fact, even for the lowest temperature (1-37° K.), 
the flow rate at 90 per cent saturation is only 5 per cent of the full rate ; 
while nearer the lambda point (1-84°K.), this fraction of the full rate is only 
reached at 98 per cent saturation. The curves also indicate that for any 
particular percentage saturation, there exists a temperature above which 
no film transfer occurs. This “ critical ’’ temperature decreases with 
decreasing percentage saturation. 


0:10 


0-05 


60 8 100 % 
‘ SATURATION 


Our values give mass flow rates which can be expressed as a product of 
the film thickness d, and an average velocity v. The only measurements of 
d close to 100 per cent saturation (Long and Meyer 1949) indeed indicate 
that this quantity is rising rapidly above 90 per cent saturation. At 
80 per cent saturation, where according to Frederikse and Gorter (1950) 
the film thickness is about 10 per cent of that of the saturated film 
(Daunt and Mendelssohn 1939, Burge and Jackson 1951), d hardly varies 
with temperature. On the other hand, it is evident from our measure- 
ments that the mobility of the film increases rapidly as the temperature is 
lowered. It is hoped to obtain better determinations of v from improved 
measurements of the adsorption isotherms near 100 per cent saturation, 
and from more detailed determinations of the transfer rates which are now 
in progress. So far, our observations do not suggest that v xd is constant. 
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The Upper Limit for the Neutrino Rest Mass. 


By O. Kororp-HANSEN, 
Institute for Theoretical Physics, University of Copenhagen*. 


{Received July 9, 1951.] 


Since the B-spectrum of H? has a very low maximum energy it provides 
us at present with the best possibility for the determination of the rest 
mass of the neutrino. The following comments on the problem may be 
of relevance in view of the great experimental progress recently achieved 
(Curran, Angus and Cockroft 1949, Hanna and Pontecorvo 1949). 
For a finite neutrino mass the shape of the 6-spectrum is given by 
P(E)=p.E. (E™*+m,—E) . {(E™*+m,—E)?—m=}# 

x [1+bm,m,/(E(E™*+-m,—E))]. F(Z, E) . . . (1) 
where E is the energy of the B-particle including the rest mass m,, p is 
the corresponding momentum, E™* is the maximum energy of the 
B-particle, m, is the rest mass of the neutrino in energy units (we put 
C= 1)e NZ: E) is the Coulomb factor and 6 is a constant which we shall 
discuss in the following. 

It has been shown previously (1947) that even if the neutrino mass is 
not zero the extrapolated maximum f-energy, if sufficiently large, is the 
same as for vanishing neutrino mass, 7. é., it equals the mass difference 4 
between the two nuclei involved in the process. Consequently, in this 
case the maximum f-energy and the mass of the neutrino can only be 
determined by the original Fermi method, 7. e. from the shape of the 
f-spectra in the immediate vicinity of the upper limit. 

Pruett (1948) has pointed out that the term in (1) containing 6 is 
important for the discussion of the influence of a finite neutrino mass 
at low energy f-spectra. Pruett adopted the value —1 for 6. However, 
6 may also assume the value +1, and if linear combinations of the five 
coupling cases are used also values between --1 may occur. Furthermore, 
no definite sign of 6 can be attached to a definite coupling case (Yang 
and 'Tiomno 1950). 

Many authors assume that the result of an extrapolation of the observed 
B-spectrum when the term containing b in (1) is taken into account and 
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for low maximum energies leads to an extrapolated upper limit of the 
spectrum which deviates by m,/2 from the value obtained for vanishing 
neutrino mass. Such a result is found when the f-spectrum is expanded 
in a series in m,/(E™*+m,—E). Unfortunately, however, this series is 
not adequate for the low energy spectra since the condition for such a 
series development, viz. m,<E™*+m,—E, can not be assumed to hold 
for the major part of the spectrum. A correct procedure is to apply a 
numerical evaluation of (1). The result is illustrated in fig. 1 where three 
hypothetical Kurie (or Fermi) plots are given. In this kind of plot one 
draws 


(P(E)/(pE . F(Z, E))}! SUT a ote? We ee gore ia) 


Fig. 1. 


[—> keV 


figure shows hypothetical Kurie plots for (4—mg)=20 keV. and for 

ag Be hine eetinde mass (curve 1) and for m,=2-5 keV. and b=] 
(curve 2) and 6=—1 (curve 3). The extrapolated curves will cut the 
abscissa very close to 20 keV. for curve 2 giving no effect of m, and at 
~18 keV. for curve 3 giving an effect of ~4/5m,=m,. Tis the kinetic 
energy of the B-particle. 


as a function of E. If m,—0 the Kurie plot is a straight line cutting the 
abscissa at E=E™*. For m,40 effects of the neutrino mass appear. 
In the present low energy case the conclusions as regards the mass of the 
neutrino are that if b=1 one finds a curve which lies rather close to the 
curve for zero neutrino mass, and in the case b=—1 one finds a curve 
which deviates from the curve for zero neutrino mass by almost the full 
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amount corresponding to m,. This means that if ever a finite neutrino 
mass is observed then this procedure could be used for a determination 
of b. 

In the case of H® one does not know the precise value of the mass 
difference 4 between the H® and He® nuclei. The main uncertainty 
comes from the uncertainty in the n-p mass difference. Nuclear cycles 
involving B-decay processes are our best tools for the determination of 
the n-p mass difference (Li, 1951) but obviously they involve the 
assumption of zero neutrino mass so that we cannot apply the 4 value 
found in this way. Consequently we cannot compare E?* and J directly. 
However, the shape of the spectrum is known very precisely close to the 
upper limit of the energy, and thus from the shape of the spectrum one 
has the possibility of a determination of the neutrino rest mass. The 
experimental result is (Hanna and Pontecorvo 1949) that the Kurie 
plot deviates from a straight line by less than 1 keV. Since 6 is not 
known this deviation may be a measure of ~1/5 of the neutrino mass 
only. This represents the lowest curve in fig. 1 corresponding to b=—1. 
Thus we may only conclude that the neutrino mass is smaller than 5 keV. 
or that the neutrino mass is smaller than 1/100 of the electron mass. 
This statement is more conservative than the result given by Hanna 
and Pontecorvo (1949). 
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ERRATA 


A Theoretical Investigation of the Compression of a Ductile Material between 
Smooth Flat Dies, by A. P. GREEN, 1951, Phil. Mag., 42, 900. 


On page 916 the tranformation equations should read : 
u=y’, v=—2Z', ara’, d=¢' (0d), 


ae ‘oe —y’, a>B’, = —¢’ (c), . . . (5) 
UY", v=2', (o haoe ae € o= —¢’ (d). 


Electron and Gamma Ray Spectroscopy with Scintillation Detectors, by R. C. 
BANNERMAN, G. M. Lewis and 8. C. Curran, 1951, Phil. Mag., 42, 1097. 


On page 1102, line 6, the expression “15 per cent’ should read “ 1 per cent”. 


[The Kditors do not hold themselves responsible for the views 
expressed by their correspondents. | 


OBITUARY NOTICE 


Dr. ALLAN FERGUSON 


We regret to announce the death of Dr. ALLAN FERGUSON, for 
many years an Editor of the Philosophical Magazine, who passed 
away on Friday, November 9th, 1951. 

An appreciation and tribute will be published in a later issue of 
this publication. 
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